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ABSTRACT 


A comprehensive review of the general theory of thermal radiation from 
liquid and solid propellant rocket exhausts is presented. More than 600 refer- 
ences are discussed and almost 500 equations describe the details of gaseous 
radiation, radiation from carbon particles, and radiation from aluminum oxide 
(A1 2 0 3 ) particles. The equation of radiation transfer is derived in detail, and 
solutions include the effects of scattering (Mie and Rayleigh) , as well as emis- 
sion and absorption. The analyses of spectral line broadening and band models 
are presented for gaseous radiation, and formation and sizes of carbon and 
A1 2 0 3 particles are discussed for particle radiation. Numerous methods of 
predicting gas and particle emissivities and radiative heating are critically 
analyzed, and the most accurate methods currently available are discussed in 
detail. Although the application of the general theory is directed toward radi- 
ation from the exhausts of the Saturn H-i, F-i, J-2 and RL-10 liquid propellant 
engines (for base heating) and of the solid propellant ullage and retro motors (fo 
stage-separation heating) , the basic theory applies to radiation from the exhaust 
of any liquid propellant engine or solid propellant motor. 
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TECHNICAL MEMORANDUM X-53579 


REVIEW OF THERMAL RADIATION FROM LIQUID AND 
SOLID PROPELLANT ROCKET EXHAUSTS 

SUMMARY 


A comprehensive review of the general theory of thermal radiation from 
liquid and solid propellant rocket exhausts is presented. More than 600 refer- 
ences are discussed and almost 500 equations describe the details of gaseous 
radiation, radiation from carbon particles, and radiation from aluminum oxide 
(A1 2 0 3 ) particles. The equation of radiation transfer is derived in detail, and 
solutions include the effects of scattering (Mie and Rayleigh) , as well as emis- 
sion and absorption. The analyses of spectral line broadening and band models 
are presented for gaseous radiation, and formation and sizes of carbon and 
A1 2 O s particles are discussed for particle radiation. Numerous methods of 
predicting gas and particle emissivities and radiative heating are critically 
analyzed, and the most accurate methods currently available are discussed in 
detail. Although the application of the general theory is directed toward radi- 
ation from the exhausts of the Saturn H-l, F-l, J-2 and RL-10 liquid propellant 
engines (for base heating) and of the solid propellant ullage and retro motors (for 
stage -separation heating) , the basic theory applies to radiation from the exhaust 
of any liquid propellant engine or solid propellant motor. 


INTRODUCTION 


Among the many problems which designers of modern rockets and 
missiles face today are the accurate predictions of radiation and convective 
heating from the exhausts from the large main-stage engines and the smaller 
ullage and retrorockets during stage separation. Excessive heating from the 
exhausts of the first-stage engines in the base region destroyed several of the 
early Atlas, Jupiter and Titan missiles during lift-off. On the other hand, too 
much insulation for the heat shield in the base region is not desirable because 
of the weight penalty this imposes. Thus, it would be advantageous to obtain 
the best possible means of predicting this base heating. 

Heating in the base region from clustered engines is caused mainly by 
radiation at the low altitudes; however, at some intermediate altitude (25,000 
to 50,000 feet) reverse flow begins to occur, and convective heating becomes 
the dominant mode of heating. For the solid propellant ullage and retro motors, 
impingement of the exhaust directly on the vehicle causes an additional heating 
problem because of the thermal and kinetic energy transfer from the hot parti- 
cles to the surface. This report, however, is concerned only with the radiation 


heating from these rocket exhausts. This includes the detailed analysis of the 
selective emission of the gases (emitting in certain spectral bands only) , the 
near blackbody continuum emission of carbon particles found in exhausts of 
hydrocarbon-fueled engines, and the graybody continuum emission (plus scatter- 
ing) associated with A1 2 C>3 particles found in the exhausts of solid propellant 
motors. 


For the past several years the Thermal Environment Branch in the 
Aero-Astrodynamics Laboratory at NASA/Marshall Space Flight Center (MSFC) 
has been analyzing, both theoretically and experimentally, the effects of heat- 
ing to the base of the Saturn booster stages (S-I, S-IB, and S-IC) and upper 
stages (S-IV, S-II, and S-IVB) caused by the exhausts of the hydrocarbon- 
fueled H-l and F-i engines and the hydrogen-fueled RL-10 and J-2 engines. For 
the past two years this branch has been analyzing the effects of heating to various 
components and structures exposed to the exhausts of the Saturn solid propellant 
ullage (S-II and S-IVB) and retro motors (S-I, S-IB, S-IC, S-II, S-IVB, and 
Centaur) . In recent months various presentations have been made and publica- 
tions written by members of the Thermal Environment Branch and its associated 
contractors on the general subject of the thermal environment of the Saturn 
vehicles caused by the exhausts of both liquid and solid propellant rockets. 

At the American Institute of Aeronautics and Astronautics (ALAA) Second 
Propulsion Joint Specialist Conference held at Colorado Springs, Colorado, in 
June 1966, Farmer et al. [1], Chu, Neimann, and Powers [2], and D'Attore, 
Nowak, and Thommen [3] presented analyses of flow fields from plumes of 
liquid propellant engines; Huffaker[4] and Carlson [ 5] presented, respectively, 
analyses of radiation from liquid and solid propellant exhausts; Rochelle [6] 
discussed experimental measurements and the theoretical correlations of heat- 
ing (radiative, convective, and particle impingement) in Saturn solid rocket 
exhausts; and Hendershot [7] and Hopson and McAnelly [8] " discussed base 
region thermal environments from the exhausts of clustered nozzles. Earlier 
this year, Payne and Jones [9] had an article published in the AIAA Journal 
which summarized the Saturn I base thermal environment from the exhausts 
of the H-l engines to the S-I stage base region and from the exhausts of the 
RL-10 engines to the S-IV stage base region. Wilson [10, 11] late in 1966 
presented applications of the short-duration technique developed at Cornell 
Aeronautical Laboratory (CAL) , under contract from this branch for predicting 
heating to the Saturn vehicles from exhausts of scale models of the Saturn liquid 
and solid propellant rockets. 


* Propulsion and Vehicle Engineering Laboratory, MSFC. 
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This report presents a comprehensive review of the general theory of 
radiation from rocket exhausts and discusses the most current methods avail- 
able for predicting radiation heating as developed by this branch, its associated 
contractors, and other agencies. More than 600 documents have been reviewed, 
and nearly 500 equations are presented which pertain to the general theories of 
gaseous radiation, radiation from carbon particles, and radiation from A1 2 0 3 
particles. This is believed to be the most comprehensive survey presented on 
this subject, although Beheshti and Thibbodeaux [12] and Robertson and Usher 
[13] have recently published summaries of various methods of predicting radi- 
ation from rocket exhausts. Several years ago Carpenter, Foreman, and 
Goldstein [14] presented a detailed review of rocket exhaust radiation and pre- 
sented numerous references through the end of 1957. 

Many Thermal Environment Branch contractors have contributed to the 
radiation programs mentioned in this report. For the gaseous radiation analysis 
some of these contractors include General Dynamic s/C onv air at San Diego, 
California; Rocketdyne at Canoga Park, California; Warner-Swasey at Flushing, 
New York; Brown Engineering at Huntsville, Alabama; and Hayes International 
at Birmingham, Alabama. Rocketdyne and General Dynamics/Convair have 
contributed extensively in the analysis of radiation from carbon particles. 
Aeronutronic at Newport Beach, California, has been chiefly responsible for 
the A1 2 0 3 particle theoretical radiation programs. Heat Technology Laboratory 
at Huntsville, Alabama, and Hayes International (using instrumentation from 
Hayes and from Block Engineering at Cambridge, Massachusetts) have been 
responsible for radiation measurements obtained from exhausts of both solid 
and liquid propellant rockets. The stage contractors - Boeing at Huntsville, 
Alabama, and New Orleans, Louisiana (S-IC stage) ; Chrysler at New Orleans, 
Louisiana (S-IB stage) ; North American at Downey, California (S-II stage) ; 
and Douglas at Huntington Beach, California (S-IVB stage) , as well at MSFC's 
Propulsion and Vehicle Engineering (P&VE) Laboratory (and its contractors, 
Rocketdyne, United Technology Center at Sunnyvale, California, and the 
University of California at Berkeley, California) - have all made valuable 
contributions to the radiation heating problems associated with the Saturn 
vehicles. Lockheed at Huntsville, Alabama, has been responsible for the 
development of the flow-field (method-of-characteristics and thermochemical) 
programs used by MSFC and many contractors for the radiation analyses. 
Aeronutronic has combined the Lockheed single-phase program with its two- 
phase program to handle the radiation from the A1 2 0 3 particles. 

Experimental programs sponsored by MSFC involving radiation measure- 
ments have been performed in altitude cells at Arnold Engineering Development 
Center (AEDC) at Tullahoma, Tennessee, for the S-II and S-IVB ullage and 
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S-IVB retro motors and J-2 engines (all full-scale) , and scale models of H 1 
engines (on the S-I and S-IB stage) with external flow. Experiments on radiant 
base heating from the exhausts of model F-l engines (on the S-IC stage) have 
been performed at Lewis Research Center, Cleveland, Ohio, in transonic and 
supersonic wind tunnels (at altitude) . Radiation measurements have also been 
made in altitude cells at Cornell Aeronautical Laboratory (CAL) at Buffalo, 

New York, which involved scale models of H-l, F-l, J-2, and RL-10 engines 
and S-II and S-IVB ullage and S-I, S-IB, S-IC, and S-II retro motors. Measure- 
ments of plume radiation have also been obtained in altitude cells at the 
Ordnance Aerophysics Laboratory (OAL) at Daingerfield, Texas, for the S-II 
ullage motors and at MSFC's Test Laboratory for the Centaur retro motors. 

Sea level radiation measurements have been obtained at MSFC’s Test Labora- 
tory for the J-2, H-i, and F-l engines (including a static firing of the S-IC 
stage which has five F-l engines) , and total heating measurements were obtain- 
ed at Rocketdyne at McGregor, Texas, on objects placed in an S-II ullage motor 
plume 20 inches from the nozzle exit. 

This report is primarily concerned with the development of the theoreti- 
cal aspects of radiation from rocket exhausts, although in various places experi- 
mental data (some of which are compared with the theoretical calculations) are 
discussed. In general, the radiation discussed here occurs in the infrared 
{X > 0.8 /i) portion of the spectrum, although occasionally radiation in the visible 
and ultraviolet will be discussed. The details of plume gas dynamics, while 
important in the analysis of radiation, are not discussed except for the effective 
particle temperatures (two-phase flow) associated with radiation from solid 
propellant rocket exhausts. Also, form factors are not discussed in detail, 
although references are included in various sections which describe their 
analysis. 

The equation of transfer of radiation, including the process of absorption, 
emission, local thermodynamic equilibrium, and scattering, is derived first, 
and certain solutions to this equation based upon various simplifying assump- 
tions are discussed. In the analysis of gaseous radiation, radiation from an 
accelerating charge as it relates to the Lorentz line shape (important in the 
analysis of band models) is discussed; the shape and broadening of spectral 
lines is then considered; various band models (such as single-line, Elsasser, 
statistical, random Elsasser, and quasi-random) are described; and, finally, 
various methods of predicting gaseous radiation from rocket exhausts are 
presented. 
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The next major section considers the analysis of radiation from carbon 
particles (such as found in exhausts of H-l and F-l engines) . In this section, 
the formation of carbon particles in luminous flames (including rocket exhausts) , 
the determination of carbon particle sizes and the carbon dispersion equation 
are analyzed; then, the Rayleigh scattering theory as it applies to scattering by 
carbon particles is developed. Various methods of predicting carbon absorption 
coefficients and emissivities are then described and finally calculations of 
radiant heating in luminous flames and rocket exhausts are discussed. 

Radiation from A1 2 0 3 particles (such as found in exhausts of solid propellant 
motors) is analyzed and a detailed discussion of the Mie theory of scatter- 
ing by these particles is given. Alumina (Al 2 0 3 ) particle sizes and distributions, 
includingthe combustion of aluminum powder to Al 2 03 particles, are then consid- 
ered, and various methods of predicting the A1 2 0 3 particle cloud emissivity are 
described. The determination of A1 2 0 3 particle cloud effective temperature is 
discussed briefly. All of these items are then combined to analyze various 
methods of calculating radiation from aluminized rocket exhausts. 

Since this report has been abstracted from a large number of sources 
and references, units of measure used by the authors are retained. Conversions 
to S. I. units are not made in order that the work of these authors be accurately 
presented and remain intact. 
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EQUATION OF RADIATION TRANSFER 


The fundamental equation of radiation transfer as it pertains to radiation 
from both gases and solid particles, such as carbon or A1 2 0 3 , is derived in this 
section. The spectral intensity is discussed, followed by a discussion of absorp- 
tion, emission, local thermodynamic equilibrium, and scattering. The equation 
of transfer is derived, and solutions to this equation with certain simplifying 
assumptions are given. 

This discussion is based upon a number of classic textbooks on radiation 
transfer theory and astrophysics, as well as on several recent documents which 
derive the equation of transfer and/or give some of its solutions. The textbooks 
consulted were Chandrasekhar [15,16], Busbridge [17], Kourganoff [18], 
Ambartsumyan [19], Hopf [20], Wooly and Stibbs [21], Sobolev [22], Rosseland 
[23], Pai [24], and Goody [25]. The documents consulted included Samuelson 
[26], Viskanta [27,28], Love [29], Love and Grosh [30], Viskanta and Grosh 
[31], Seay [32], Brown [33], Bartky and Bauer [34], Lai and Purgalis [35], 
deSoto [36], Goulard [37], Sampson [38], Tellep and Edwards [39], Hansen 
[40], Ueno et al. [41], and Churchill et al. [42]. Most, but not all, of these 
authors included the effects of scattering in their derivations of the equation of 
transfer. In the analysis to follow, the salient features in each of these deri- 
vations are put together to formulate the equation of transfer for use in rocket 
exhaust radiation applications. 


Spectral Intensity of Radiation 


This section discusses the spectral intensity as it relates to the spectral 
energy dE^ or spectral radiant heat transfer rate q^. The spectral energy 
dEx is defined as the radiant energy absorbed, emitted, or scattered in the 
wavelength interval (A., A + dA) , in a time interval (t, t + dt) , passing through 
the solid angle dcu and the unit area normal to the surface dA. This quantity 
is expressed as 

dE = I n • is du> dt dA dA = I cos 0 do; dt clA dA (1) 

XX X 

The spectral heat transfer rate q^ from a rocket exhaust, which is 
analogous to the integral of equation (1) , is also related to the spectral inten- 
sity 1^ and can be written in the form 
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( 2 ) 


q /A = / I cos 0 dcu 
A J A 
00 

The total radiation heat transfer rate would then be the integral of q^/A 
over the entire spectrum. The spectral intensity I A , which is sometimes 
denoted by I„ (frequency dependent) , is also called the specific intensity or 
monochromatic intensity and is usually written in units of Btu/ft 2 -sec -micron - 
steradian or watts/cm 2 -micron-steradian. This intensity at point P, in the 
direction s, is related to the polar angle 0, azimuthal angle cp, element of 
solid angle du>, element of surface area dA, and normal to the surface area 
n, as shown in Figure 1. 


1^ in direction of "s 



FIGURE 1. SKETCH OF SPECTRAL INTENSITY OF RADIATION 

The spectral intensity is in general a function of wavelength A, of time 
t, of solid angle u>, (or direction cosines H, m, n) , and of position coordinates 
x, y, z, and thus can be written as 

\ = \ (x, y> z; m, n; t) . (3) 
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The integrated intensity I is merely the integral of the spectral inten- 
sity I A and can be written as 



Absorption 

As a pencil of radiation travels through matter (such as a rocket exhaust) 
it is usually weakened by absorption. If the decrease in spectral intensity is 
dI A , and the coefficient of absorption is k x , then the following relation (experi- 
mentally proven) holds true: 


dI x = -Vk ds - 

where s represents the distance traversed along a pencil of radiation as shown 
by the vector s in Figure 1. This expression is valid for both line and contin- 
uous absorption. The coefficient of absorption can also be expressed in terms 
of a mass absorption coefficient k as 


where p is the density of a rocket exhaust at a particular point in the plume. 
Also, the notation k for the frequency- dependent absorption coefficient is 
sometimes used. (See the section entitled "Gaseous Radiation, " and subsection 
"Band Models. ") 

The absorption coefficient is a property of the particular rocket exhaust, 
and in general, depends upon temperature, pressure, and exhaust gas compo- 
sition, in addition to wavelength and time (for non-steady-state conditions) . 

The absorption defined in this section is designated as "true absorption," that 
is, absorption which may be converted into heat which may be emitted at another 
wavelength. Some authors prefer to consider another type of absorption as 
scattering, which is merely the redistribution of incident energy into another 
direction but at the same wavelength. Scattering is discussed in detail later in 
this section. 


8 


An example of "true absorption" is the absorption of radiation of fre- 
quency t> mn corresponding to a quantum transition from the lower state m to the 
higher state n of the atoms comprising the particular medium. This transition 
can be represented according to Bohr's Theory as 


hv = E - E , (7) 

mn m n 

where h is Planck's constant, v mn = c/\ mn (where c is the velocity of light) , 

and E and E are energies of the higher and lower states, respectively, 
m n 

A further coefficient, Bj^, which is associated with absorption, is 
known as the Einstein absorption coefficient or the Einstein transition probability 
and is defined in the following manner: 

probability per unit time that an atom in state m 

absorbs a photon of energy hi; and goes to state 

n with intensity I . 

A 

ran 

By assuming a purely absorbing medium and integrating equation ( 5) 
with insertion of equation (6) , a new quantity called the optical depth r ^ may 
be defined. The integration of equation (5) is known as Beer's equation and 
may be written as 


B K 

mn X 


mn 


V s * • 


( 8 ) 


where I^(s) is the intensity after the radiation has traveled a path s through 
the medium. Equation (5) is usually written in the form 


Us) = I (0) e~ 

A A 


(9) 


where the optical depth of the medium is 




pn _ ds 
m, A. 


( 10 ) 


A more general relation for optical depth is derived in the subsection 
entitled "Scattering, " where the effects of scattering are included. 
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Emission 


The process of emission of radiation from a medium such as a rocket 
exhaust can be characterized by an emission coefficient j^. This coefficient 
is defined as the spectral radiant energy rate dq^, emitted by an element of 
mass dm in the solid angle do; during the time interval (t, t + dt) , and in the 
wavelength interval (A, A + dA) , and is given as 

dq A 

^A dm da) dA 

The emission coefficient depends on the wavelength, location of unit 
volume, exhaust gas composition, temperature, and time (for non-steady- 
state conditions) . In general, there are two types of emission processes: 
spontaneous emission and stimulated (or induced) emission. 

The first type of emission occurs when atoms in a higher state, n, 
undergo a quantum transition to a lower state, m, emitting a quantum of energy, 
ht'nm- For this spontaneous emission, another emission coefficient, A nm , can 
be defined. This coefficient is known as the Einstein spontaneous emission co- 
efficient or Einstein spontaneous transition probability, and is defined as the 
probability per unit time that an atom in state n makes a transition to state m. 

Stimulated emission (sometimes referred to as negative absorption) 
occurs when the transmitted intensity becomes greater than the incident inten- 
sity. This process occurs when a condition is artificially created in which the 
number of atoms per unit volume N R in a higher-energy state exceeds the num- 
ber of atoms per unit volume N m in a lower-energy state. This is accomplish- 
ed by adding energy either by absorption or collision. The stimulated emission 
then takes place when the atoms drop to the lower states and emit photons cor- 
responding to an energy difference between the two states. When an atom in 

the state n has been exposed to radiation of intensity K , it is possible to 

nm 

define a stimulated emission transition probability per unit time as B^, such 
that 

^nm \ ~ P robabilit y P er time that the atom makes the 

nm n— m transition and emits energy in the direction 
of I x . 
nm 
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Stimulated emission is very important in the new science of lasers, 
which are discussed by Garbuny [43] , Brown [44] , and Lengyel [45] • Successful 
application of a laser requires the utilization of a system of molecules or atoms 
in excited states in such a manner that the stimulated emission processes 
(obtained through population inversion) are enhanced while the energy losses 
are kept low enough to produce an amplification greater than unity. 


Local Thermodynamic Equilibrium 


Through the concept of local thermodynamic equilibrium (LTE) , estab- 
lished by Kirchoff, it is possible to relate the emission coefficient j A to the 
absorption coefficient k a and to the Planck blackbody intensity of radiation 
B^(T) . This quantity, B A (T) , will be derived in terms of the transition prob- 
abilities A nm , B nm , and B mn , defined above. 

Kirchoff stated that, in general, any closed system which is not experi- 
encing any change with time is in local thermodynamic equilibrium. That is, 
in the system the temperature is everywhere constant, and if the system contains 
more than one substance, it must also be in chemical equilibrium such that there 
can be no sudden change in internal structure. Furthermore, the system must 
be in mechanical equilibrium such that there can be no macroscopic movement 
within the system itself and also between the surroundings and the system. 


Kirchoff also postulated that, if I A is the intensity at any point in the 
medium whose refractive index is n A , then I A /n A is a constant value through- 
out an enclosure containing the medium. He also stated that the quantity 
j x /(n A Km A ) , which is constant throughout any one enclosure, is the same for 
the other enclosures at the same temperature and is a universal function of 
temperature. This relation is expressed by the Kirchoff Law 


= n \ K , B (T) 
X m,A X 


( 12 ) 


The Planck blackbody function B A (T) may be derived in terms of the 
Einstein transition probabilities A nm , B nm , and B mn in the following manner 
as discussed in a previous report [46]. A more detailed treatment of transi- 
tion probabilities is found in the classic paper by Milne [47] . The total emis- 
sion (sum of spontaneous plus stimulated) may be expressed as 
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(13) 


N 

j, = — — (A + B I ) h v 

X p nm nm X nm 

nm K nm 

where p = density of the medium and N n = number of atoms per unit volume in 
the state n. 

The absorption coefficient may be written as 
N 


nm 


m 

P 


B hi, 
mn nm 


(14) 


where N denotes the number of atoms per unit volume in the state m 
m 

Since an equilibrium between energy states and radiation exists such 
that, per unit time, as many atoms are raised from state m to state n with 
the absorption of a photon as undergo transition from state n to state m with 
subsequent emission of a photon, it is possible to write 


N 

n 


(A + B 
nm nm 



= B N I 
mn m X 


(15) 


nm nm 

This implies that the processes n— m and m— n occur at equal rates. 

The resulting equation for I is 

A 


nm 


nm 


nm 


B 


nm 


[ (N B /N B ) -1] 
m mn n nm 1 


(16) 


From Maxwell -Boltzmann statistics, the following relation applies: 


N SL= g ffl<!XP {■ I<V E m ,/<kT, > } = 7 s ex P < ch A mn kT > . (17) 

n n v ; & n 


where g^ and g^ are the degeneracies (or statistical weights) of the m and n 
states, respectively, and k is Boltzmann's constant. 
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By combining equations (16) and (17) , the following equation for the 
blackbody radiation intensity is obtained: 


I 


nm 


= B, (T) 


A nm B nm * B nm 6xp * ch//A nm kT ^ 1 A nm (18) 


as 


Another expression for B% (T) can be obtained from quantum statistics 

nm 


2c 2 h 


2 C 1 


nm 


A [exp (ch/AkT) -1J A 5 [exp(C 2 /AT) -1] A 


= B (T) 


nm 


(19) 


Thus, it may be seen that 


2c 2 h 


A /B 
nm nm A 


B g 

, mn °m 

5- and — = 1 . 

B g 

nm n 


( 20 ) 


By noting that v\ - c and 


B (T) dp = - B (T) dA or B (T) B (T) , 

^ A v c A 


( 21 ) 


equation (19) can be written in terms of the frequency as 


B (T) = 

p' c[exp {hp/kT) -1] 


( 22 ) 


Hence, equation (12) can be written as 


k = K 


| 2c 2 h 


A A m,A I A® [exp (ch/A kT) -1] j 


(23) 


or in terms of the frequency, p, as 
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, 2 ) _2hj/^ 1 

V - n \ K m, v | c 2 [exp (Im/kT) -1] 

The Planck blackbody function behaves differently at the two extremes 
of the spectrum [A— °°(r— 0) or A— 0(r— °°)]» 

For the first case (Rayleigh- Jeans distribution) : 

B^(T) =2kTc/A 4 or B^(T) = 2kTi^/c 2 . (25) 

X— oo v-~ 0 



For the second case (Wein Distribution): 

B (T) = (2hc 2 /A 5 ) e ' hc/AkT or B^(T) = (21m 3 / c 2 ) e" h " /kT . (26) 

A— *■ 0 


By integrating equation (19) over wavelength or equation (22) over 
frequency, the following equation results for the intensity of blackbody radiation 
B(T) : 


a 


B(T) = f B (T) dA T 4 , 


(27) 


where a is the Stefan-Boltzmann constant and is equal to 
B 


<T_ = 


2tt 5 k 4 


B 15 c% 


,2u3 


(28) 


Scattering 


o 

If there are particles in a medium (having a radius of the order of 100 A 
or so) , a pencil of radiation of intensity I A may be weakened not only by absorp- 
tion, but also by scattering. This scattering phenomenon includes the combined 
effects of reflection, refraction, diffraction, and transmission of radiation by 
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the particles. In the exhaust of hydrocarbon-fueled rockets, scattering of radi- 
ation by carbon particles whose radii are much less than the wavelength of radi- 
ation is characterized by the Rayleigh theory of scattering. In the exhaust of 
aluminized solid propellant rockets, scattering of radiationby A1 2 0 3 particles whose 
radii are of the same order of magnitude as the wavelength of radiation is char- 
acterized by the Mie theory of scattering. These theories are described in 
detail in the sections entitled "Radiation from Carbon Particles" and "Radiation 
from A1 2 0 3 Particles," respectively. 

The type of scattering considered in this section is single scattering. 
Problems in multiple scattering, briefly discussed in the section entitled "Radia- 
tion from A1 2 0 3 Particles," involve the scattering of a diffuse field of radiation 
by very large numbers of particles which in turn absorb and scatter radiant 
energy originating or scattered from other particles. Also, only coherent scat- 
tering, i.e. , scattering in which the absorbed quantum is re-emitted in the same 
wavelength, the re -emitted and absorbed quanta differing only in direction, is 
considered in this analysis. 

A scattering coefficient a ^ may be used to characterize a medium if, 
from a pencil of radiation incident on an element of volume of height ds and 
area dA, the amount of energy scattered from it in all directions is 

a I cos 0 dA ds du> dA dt , (29) 

A A 

in which a mass scattering coefficient cr can be defined as 

m, A 


cr = p a 
A m, A 


Since the mass of the element of volume is 


(30) 


dm = p cos 0 dA ds , 


(31) 


expression (29) may be written in the form 


cr I dm do> dA dt . (32) 

m,A A 
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The radiant energy scattered from a pencil of radiation will contribute 
to the radiant intensity in another direction. To assess fully the overall effects 
of scattering, it is necessary to account for the direction of scattering. In 
general, a phase function, p(cos ©) , may be defined to account for the direc- 
tional dependence of scattering. This function may be called a scattering function 
or an amplitude function and may be written differently [such as p(0 , , <p') ] 

by various authors. The coordinate system for the phase function used in this 
analysis is shown in Figure 2. 


Z 



FIGURE 2. SKETCH OF COORDINATE SYSTEM FOR RADIATION 
SCATTERED BY A PARTICLE 

The phase function will be represented in the following manner by observ- 
ing Figure 2. The rate at which energy is scattered by the element of mass dm 
from a pencil of rays of intensity I^( s, p, 4>) in the direction (p,0) and having 
a solid angle dto with a wavelength interval dA into a solid angle dco' in the 
polar direction O' and the azimuthal direction </>' , is 
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(33) 


. I. (s p(cos ©) 

HI , A A 


~ — dA du> dm . 
4 tt 


The angle © is the scattering angle between the incident intensity 1^ 
(s and the scattered intensity 1^ (s , m' , ') whose cosine may be deter- 

mined as 


cos© = cos0 cos0 + sin0 sin0' cos (<p'-<p) - ny.' + (1 -p 2 ) 


i 

, 2\2 


(34) 

1 

(l-jU ,2 ) 2 COS (<+>'-</>) 


where p = cos 0 and ju* = cos 0 ' . 

Furthermore, the quantity p(cos 6) is a probability function which 

expresses the probability that a photon in the wavelength interval dA incident 
on the element of mass dm, from a direction within the solid angle dcu, char- 
acterized by angles 0 and cp, will be deflected into a direction within the solid 
angle dco* characterized by the angles 0* and cp' . 

The rate of loss of radiant energy from the pencil of rays caused by 
scattering in all directions may now be written as 


4tt 


m, A 


dA dm do) I^(s,p, cp) J p(cos ©) 


da) 1 

47T 


This formulation agrees with expression (32) if 


(35) 


. au> 

J P( cos ©> 77 = i; 
0 


(36) 


that is, if the phase function is normalized to unity. However, it must be 
stated that equation (36) applies only for no absorption and for perfect isotropic 
scattering; i. e. , the magnitude of the radiation intensity is independent of the 
scattering angle ©. 

In the general case where both true absorption and scattering are present, 
equation (36) is less than unity. An albedo for single scattering oP Q may there- 
fore be defined as the fraction of radiation lost from an incident pencil of radi- 
ation. Hence, l-co 0 would be the remaining radiation transformed into other 
forms of energy or of radiation of other wavelengths. The albedo for single 
scattering is thus defined as 
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=£ 1 , 


(37) 


-4jr dw' 

/ P(cos ©) — 

0 


= CO 

o 


a 

m,A 

B 

m,X 


where a mass extinction coefficient /3 m ^ has also been added. This mono- 
chromatic mass extinction coefficient is equal to the sum of the monochromatic 
mass absorption and scattering coefficients, as 


/3 = k + a . (38) 

m,\ m,A m, X 

For anisotropic scattering the phase function may be written in terms of 
the Rayleigh phase function (if the radii of the particles are very small com- 
pared to the wavelength) as 

p(cos ®) = 3/4 (1 + cos 2 ®) . (39) 

It may also be written in terms of the Mie coefficients if the radii of the parti- 
cles are of the same order of magnitude as the wavelength (discussed in the 
section entitled "Radiation from A1 2 0 3 Particles") . The phase function may 
also be expanded in a Legendre polynomial series as 

OO 

p(cos ©) = Y co p (cos ®) , (40) 

u n n 
n=l 


such as used by Churchill et al. [48] , Chu and Churchill [49] , and Brown [33] . 


Formulation of the Equation of Radiation Transfer 


The results of the five previous sections will now be put together to form 
the equation of radiation transfer, which describes the radiation field in an 
isotropic medium. This medium absorbs, emits and scatters thermal radiation 
at a wavelength X. In deriving this equation, it is necessary to consider four 
classes of photons, two of which are lost from the radiation field in the direction 
(p, (p) by scattering and absorption, and two of which are gained by the radiation 
field in the direction (p' ,</>') by scattering and emission. Figure 2 should be 
used for the following discussion. 
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The first class of photons considered is that made up of photons incident 
upon the element of mass dm, in time dt, and in the direction (p, 0) contained 
in the solid angle dco which is singly scattered into the solid angle dco' in the 
direction (p 1 , 0 1 ) by interactions with particles in the mass element dm. The 
loss of radiant energy is given as in equation (35) with the change that the phase 
function p(cos ©) will be written p(p,0; ju' , </>’> , denoting scattering through the 
angle between the directions (p,0) and (p',0*). Also, the solid angle dco' is 
set equal to sin o'do'dty' . Hence, the rate of loss of radiant energy caused by 
scattering may be written as 

°m A it 2t 

" 4j J dm dA'dco I^(s,p,<£) J J p(p, 0;p\ 0’)sin 9' dd' d<p' . (41) 

0 0 


This relation may be simplified if equation (36) is used to choose a phase 
function such that 


47T 


27T 

p (p, 0;p\ 0') sin O' d0 d0 = 1 . 

0 0 



(42) 


In this case, equation (41) may be written similarly to equation (32) as 

°m A V s ’ dm ^ dc0 ' (43) 

The second class of photons is that incident upon dm in time dt and 
direction (p, 0) contained in the solid angle dco and absorbed by the individual 
particles in the mass element dm. The rate of loss of energy by absorption 
may be expressed similarly to equation (32) , with the mass scattering co- 
efficient being replaced by the mass absorption coefficient as 


K „ ^ K ( S »M,0) dm dA dco . (44) 

III , A A 

The third class of photons considered is that incident on the mass element 
dm, in time dt and in the direction (p', 0’) contained in the solid angle dco', 
which is singly scattered in the direction (p, 0) into the solid angle dco. This 
change in energy is actually a gain whose rate of change may be expressed as 
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(45) 


°m, A 
47T 


tt „2ir 


dm dA dco J J I ( s, ju', <p') p(p, <p') sin 9 ' d0 ' dcp' . 


0 0 


The fourth class of photons is that thermally emitted from the particles 
contained in the mass element dm, into the solid angle dco in the direction 
in,cp) during the time dt. This change in energy also represents a gain whose 
rate of change is expressed as 


j dm dco dA . (46) 

J A 

Calculation of this contribution may be modified by considering the ele- 
ment of mass to be located in a perfectly insulated isothermal enclosure which 
is maintained at a constant equilibrium temperature, T. Since within the en- 
closure the radiation field is isotropic and is in equilibrium with its surroundings, 
the rate of energy which would be emitted from dm in the wavelength interval 
(A, A + dA) and in the direction (p,0) contained in dco, upon an instantaneous 
removal of the enclosure walls, is 


n 2 k B (T) dm dA dco (47) 

A m,A A 

where equation (12) has been used. 

These four classes of photons or changes in energy may be summarized 
as follows: (1) scattering of a fraction of incident photons into dco', or redirec- 

tion of a fraction of incident energy into dco' (loss) ; (2) absorption of a fraction 
of incident photons coming from dco or decrease in incident energy (loss) ; (3) 
scattering of a fraction of incident photons into dco or redirection of a fraction 
of incident energy into dco (gain) ; and (4) emission of photons from the mass 
element dm into dco or increase in energy (gain) . 

The sum of the above changes in energy should be equated to the net 
change in radiant energy crossing a cylindrical element of cross-section dA and 
height ds. The net change in radiant energy crossing the faces normally in the 
wavelength interval (A, A + dA) , in time dt, and confined to the element of solid 
angle dco is given by 

dl . .. 

_^_( s >Mj <£) dA dco dA dt . (48) 

dt 
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The total derivative 


dl / - 

A (s ,n,(p) 

dt 


may be written as 


dl 0 I 01. 

__^(s,p, 0 )_ __A (s,p, 0) + c A_(s,P,0) . 

dt 9t 9s 


( 49 ) 


The quantity dm occurring in equations (43), (44), (45), and (47) may 
be written as dm = p ds dA, since in this case of radiation normal to the faces 
of a cylinder, the factor cos 0, appearing in equation (31) , may be equated to 
unity. Also, the quantity ds may be written as c dt . 


Thus, equation (49) may be set equal to the sum of equations (43) , (44) , 
(45) , and (47) , with the result as 


01 01 

__A (s,p,0) A(s,m,0) 

at as 


do> dA dt dA = < -(a + k ) I (s,p, 0) 

| m , A m,A A 


cr 7T 2 7T i 

n l K m A B a [T( s)] + ~~ 4 tt f f p(m,0;m',0) sine' de’d^'j 


c dco dA dt p dA . 


(50) 


Cancelling out the quantities cdu> dA dt dA and assuming steady-state 
conditions along the pencil of rays in the direction of the vector s, the equation 
of transfer may finally be written as 


dl 


-< vV y^.o» + 


G ^ 7T 27T 

+ “TT f f M s V,0') p(p,0;p',0') 
4 0 0 X 


n 2 . k B [ T ( s) ] 

AAA 

sin 0 ' do 1 d0 ' , 


(51) 


where the substitutions k, = p k , and a - p o have been made. 

A m, A A m, A 
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Solutions to the Equation of Radiation Transfer 


The final form of the equation of radiation transfer, equation (51) , is 
an integro-differential equation, the solution of which is quite formidable unless 
certain simplifying assumptions are made. The following four simplifying 
assumptions are briefly discussed in this section: (1) use of monochromatic 
effective absorption coefficient, (2) purely absorbing medium, (3) purely 
scattering medium, and (4) plane-parallel case. 

1- Use of Monochromatic Effective Absorption Coefficient. One simpli- 
fying assumption which may be helpful in solving equation (51) is to define a 
monochromatic effective emission coefficient, j e A , as the radiant energy 
leaving the element of mass dm in the medium per unit solid angle, per unit 
wavelength, and per unit of time as 

CT y (S) r 7T 2 7T 

~ + 47r J J I A ( s >/ i »0') P(Ab0;M\c/>') sino'do'd^' 

(52) 

where j A (s) may be defined from equation (12) , if the element of mass dm is 
in local thermodynamic equilibrium. 

The absorption and scattering coefficients are assumed to depend on 
position s only. When the direction of propagation s is specified such that 
there is no dependence of intensity on ju' and 0’, equation (51) becomes the 
following inhomogeneous linear first-order differential equation in the one 
variable , s 

dl ( s) 

= y s) +j e , A (s) . (53) 

This equation may be readily solved if both sides are multiplied by the 
integrating factor exp[ f/3 ( s) ds]. Hence the equation takes the form 


ds I I A (S) ex P[/^ A ( s ) ds J j = J e>A ( s > ex Pl dS] • (54) 
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Integration of this equation and division by exp[ J /^(s) ds] results in 
the following solution: 


I^(s) =Cexp[-//3 A (s) ds] 


+ exp[-J/? A (s) ds] fj & ^(s') exp[ f/3^( s) ds] ds* 


(55) 


where s' is a point along the vector s' between s = 0 and s = s, and the optical 

thickness r (s) is 
A 


t (s) = f B (s) ds = J (k 
0 0 


X + V 


ds 


(56) 


The constant C can be determined from appropriate boundary conditions. 
Equation (55) is still not a true solution of equation (53) unless the effective 
emission coefficient j e ^ is a known function of the intensity I^(s) . Then 
equation (55) would be converted into an integral equation for intensity, for 
which an approximate solution could be obtained by numerical integration. 

2 . Purely Absorbing Medium. For a purely absorbing medium with 
very small particles (of the order of the size of molecules) , the scattering 
coefficient is approximately zero, and the index of refraction n^ of the 
medium is approximately unity; hence, (3 ^ » k ^ and j e ^ ~ j^. Therefore, 
equation (53) reduces to 


dys) 

ds 


- EyT(s)] + j A (s) . 


(57) 


This equation is sometimes written in the form 


1 


K 


X 


dys) 

ds 


- - B a [T(s)] + J A (s) , 


(58) 
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where the function <y s) is Chandrasekhar^ "source function" defined as 

j x < s > 

J X (S)= 7>)' (59) 

A 


The solution to equation (57) or (58) , obtained by using an integrating 
factor similarly to the solution of equation (53) , is 


I A (s) = C exp[-J>ys) ds] 


+ expl-jys) ds] Jy s') exp ( fie (s) ds] ds* . 


(60) 


Using the boundary condition, I (s) = I (o) at s = 0, equation (60) 
becomes A A 


ys) =y0) e 


-t\(s,0) 


/£<•’> 

0 


- t A (s,s') 


ds' 


(61) 


where the optical thickness r (s, s') is defined as 

A 

y(s,s') = f K (s) ds . 
s' 


(62) 


It might be mentioned at this point that another quantity, the emissivity 

at the surface of a medium, e (/x) , may be expressed as 

A 




V°’«> 

Ey t) » 


(63) 


where M - cos0. The emissivity is an important factor in the calculation of 
heating rates from rocket exhausts. In general, the principle of conservation 
of energy holds such that the following relation is true: 
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1 = a (m) + p^(m) + t ^(M) , 


(64) 


where a (p),p (p) , and t (p) are respectively the absorptivity, reflectivity 
\ A A 

and transmissivity of a particular medium. If Kirchoff’s law holds, 


« (/-) = e^(p) . 


(65) 


In the present case [equation (57)] , if no scattering is present, p^(p) = 0, 
and the following expression may be written for the emissivity 


£ A (M) = 1 - t A (p) = 1 



( 66 ) 


Hence, equation (63) may be written as 

I A (0,M) = (1 - e" Tx/M ) B a (T) , (67) 


where this value of H (0) would be inserted into equation (61) to obtain the 
intensity of radiation entering the medium. 

3. Purely Scattering Medium. One other simplified solution of equation 

(53) may be obtained if the absorption coefficient /c^( s) is set equal to zero. 

This may be the case for a purely (isotropically) scattering medium in which 
ft ( s) =ct^(s), k ^( s) =j\(s) =0, and the phase function p( p,0;p',<|>') = i. 
Hence, equation (53) takes the form 


dl (s) %( s ) 

-57- = -V s > V s)+ ^7- 


The solution of this equation is 


n ^ ^ y. 27T tit 

J f /-^( 3 * S ) d s sin0 d0 d<£ . 
0 0 0 


( 68 ) 
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ys) = C exp [ - J <y s) ds] 

( 69 ) 

+ exp[-Ja^(s) ds] Ja^(s') ys') exp[Jcys) ds] ds' . 
Applying the boundary conditions I (s) = I (0) at s = 0, equation (69) becomes 

A A 


ys) = y 0) e 


- t a (s,0) 


+ 


kf B 'f *!%<*'> \ {s ' } e_TA(S,S,) ds' sinO'de'dd)' , 
^000 


(70) 


where the optical thickness t (s,s') is defined as 

A 

T x (s, s') = j <y s) ds . (71) 

s' 

4. Plane -Parallel Case. One final solution of the equation may be 
described briefly. This is the plane-parallel case which has certain applica- 
tions in determining the thermal radiation from a rocket exhaust (by dividing 
the exhaust into various homogeneous isothermal layers) . 

In this case the element of length ds will be represented by dz/p, where 
p = cos 0. The equation of transfer, equation (51) , may now be written in the 
form 


dI 

a A (z) -2ir 


(72) 


+ — — f f I (z ,p',<£') p (p.tfqp'.^') sin 0 ' d0 ' d<£ ' 
47r 0 0 A 


If the scattering is isotropic, the integral on the right becomes 
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( 73 ) 


7T 2ff 1 

f J I (z,p' A) sin0' d0' d<p' = 2it f 1 (z,n'<p') dp' . 

0 0 A _1 A 

If two parallel planes or layers, z = 0 and z = z t , can be drawn close 
enough together, as in Figure 3, in a rocket exhaust such that the composition 
between the two layers is homogeneous and at some averaged temperature, the 
intensities on the boundaries may be specified as 1 


I x i z ,V,<p) = 1^(0, M,</>) at z = 0, n < 0 


and 


(74) 


I a ( z >M, 4>) = I A (z,p,(£) at z = Zj, jn > 0 . 




‘ ' r r r r r 

... ii 

A 



\ ds, 

**-n - — 

r s 

~F 

dz = ^ ds 

i_ 




fi = cos 6 


FIGURE 3. SKETCH OF COORDINATE SYSTEM FOR PLANE-PARALLEL 

CASE 


Hence, equation (72) may be written as 


dl , 

A- ( z.M.fl) 

dz 


Vf! 

M 


I A (z,p) + 



M 


(75) 
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where equation (52) has been used. By using the integrating factor 
T fl (z) dz 

exp J , the solution of equation (75) may be obtained as 


/3^(z) dz 


I (z ,n,<p) = C exp - J 


j8 (z) dz B{ z) dz 

eX P ~f u fW Z ’ >eXP f At 


Applying the boundary conditions of equation (74) , the solution is 


i a (z,m,4>) 


I (z u n,<p) exp -J 


Z\ M z > dz 


Zl Z' P (Z) dz , 

+ f i (z) exp f for 

J J e,X J M M 

z z 


I (O,m,0) exp -J 
0 


z t /3 a (z) dz 


+ jV (z) exp / 


Z ’V Z) dZ | dz' 


for l<p<0 


There are many other ways (including numerical approximations) in 
which the equation of transfer may be solved. However, it is not the purpose 
of this report to go into all the detailed solutions to the equation of transfer by 
various authors. Rather, it is to study the theories of absorption and scattering 
and attempt to predict emissivities of the rocket exhausts so that thermal radi- 
ation from the exhausts may be calculated. 


GASEOUS RADIATION 

This section presents gaseous radiation, including band models and 
absorption coefficients, from liquid propellant rocket exhausts. The principal 
gaseous molecular emitters in exhausts of hydrocarbon-fueled engines such as 
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the H 1 and F 1 are C0 2 , CO, and H 2 0, while only H 2 0 emits appreciably in 
the exhausts of hydrogen-fueled engines such as the J-2 and RL-10. Gaseous 
HC1 is also an emitter in solid propellant motor exhausts in which NH4CIO4 is used 
as the oxidizer in the propellant. Gaseous radiation from molecules is a result 
of electronic energy transitions in the visible and the ultraviolet wavelength 
region, vibration-rotation energy transitions in the near infrared region, and 
pure rotational energy transitions in the far infrared region. This type of 
radiation is emitted in the form of band spectra as opposed to radiation from 
hotter solid particles such as carbon or A1 2 0 3 , which is emitted over a con- 
tinuous spectra. 

The general subject of gaseous radiation has been discussed by hundreds 
of authors in the last 50 or so years, and no attempt will be made to give a 
comprehensive treatment of it here. Instead, a generalized review of its 
important aspects as applied to rocket exhausts (such as determination of 
spectral absorption coefficients) will be given. Excellent references which 
review problems in this area include Yossa [50], Huffaker [4 and 51], Ferriso 
etal. [52], Thompson [53], Simmons [54], and Ritland et al. [55]. ’in this 
section the radiation from an accelerating charge will first be mentioned; the 
shape and broadening of spectral lines will then be discussed; the concept of 
band models will then be introduced; and finally solutions to the equation of 

transfei as they apply to radiation from liquid propellant engine exhausts will be 
treated. 


Radiation from an Accelerated Charge 


hi this section the radiation from a charged particle which is accelerated 
in an external field will be discussed. An expression for the spectral intensity 
will be determined consequently as a function of frequency and electronic charge. 
This expression, known as the Lorentz line shape, will be further discussed in 
the next section on shape and broadening of spectral lines. The analysis in this 
section is based upon classical electrodynamical theory, of which excellent 
treatments are by Jackson [56], Panofsky and Phillips [57], Heitler [58], 

Reitz and Milford [59], Marion [60], and Somme rf eld [61]. 

^The equation of motion for a charged particle which is being accelerated 
in the x direction by a driving or external force may be written as 

—X . 

nix — F F “i* F / ri q \ 

EXTERNAL SELF ELASTIC ' ( 79 > 
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The elastic or restoring force on the charge is 


F ELASTIC 


( 80 ) 


where G is the force constant. The self-force on the charge (the force which 
the field produced by the charge exerts on the charge itself) is 


SELF 


2 

3 



(81) 


where e is the electronic charge on the particle. The external or Lorentzian 
force is 


EXTERNAL 


eE + — (ux B) 
c 


(82) 


where u is the velocity of the charged particle, E is the electric field vector, 
and B is the magnetic field vector. 

The speed of light c has been introduced into the second term on the 
right-hand side of equation (82) corresponding to the Gaussian system of units, 
to show that this term may be neglected when u < < c. 

If the charged particle is moving with simple harmonic motion, the 
quantity x in the self-force term [equation (81) ] may be approximated by 
differentiating the expression 


as 



(83) 

(84) 


where w is the natural frequency of vibration (equal to 2-/r v ) . 
o u 

The external force will thus be written in terms of a harmonic function 
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(85) 


EXT 


- eE 


-lWt 

e 


By making the above substitutions, the equation of motion for the charged 
particle now becomes 


_ _ eE . 

o _k O -lCUt 

x + yx + or x = e 

o m 


( 86 ) 


where 




(87) 


and the second term on the left-hand side of equation (86) is thus a dissipative 
or radiation damping term. 

The solution of equation (86) may be written as 


x = E 


o m 


( - h . V 

V w - lw r / 


-iwt 


( 88 ) 


The acceleration of the charge is thus 


_e / \ -iwt 

X — E (2 2* I 6 

o m y - w - lwy J 


(89) 


The acceleration of the charge can be related to the electric field by the 
following equation 


E = cV irx X) ] ’ 


(90) 


where r is the distance from the radiation field to the accelerating charge. 

Through the use of vector calculus, it can be shown that the term in 
brackets in equation (90) reduces to x r 2 such that the above equation becomes 
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(91) 


E = 



_x 

r 


With the insertion of equation 


(89) for x, equation (91) becomes 


E = 


c 2 rm (o> 2 


(-co 2 ) 

- ci> 2 - iury) 


E 

o 


-icot 

e 


(92) 


The above case of an electronic charge which vibrates harmonically 
because of an external electromagnetic field, i. e. , equation (86) , is important 
in the classical theories of absorption, scattering, and dispersion of radiation, 
and consequently is discussed in more detail in the next section. In the case 
below, the intensity of radiation corresponding to a bound electronic charge 
following a transient disturbing impulse, but not experiencing a forced vibration, 
is derived. This case corresponds to the classical theory of spectral emission. 
In this manner, the equation of motion, equation (86) , may be written as 


x + yx + x = 0, 


(93) 


where the vector notation has been dropped for convenience. 

The solution to this equation yields a harmonic oscillation with amplitude 
x^, decaying in time as e following form: 


x = x 

o 


e -(y/2)t 


-ico 0 t 

6 


(94) 


The acceleration of the charge is thus 


. y .2 -(y/2)t -ioj 0 t 

x = x ( — ~~~ + ici> ) e r e ° 
o 2 o 


(95) 


Hence, the radiation field is 


E = 


~ 2 ~ x = — 5 - 
c r c r 


[«.<- 


V .2 - (y/2) 

-r + io) ) e r 
2 o 




(96) 
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which is a function of the natural frequency co^. 

The above equation may be written in the same form as equation ( 94) , as 


E = E 


o 


e -(v/2) e _iw o t 


(97) 


Since the oscillating charge has been radiation damped, the emitted 
radiation intensity will not be monochromatic, but will be broadened and com- 
posed of an infinite number of harmonic waves. The line width of this intensity 
distribution may be obtained by expanding the electric field, equation (97) , in 
a Fourier series as 


_ r + ^ -icot j 

E = I E e dw , 

J co 


(98) 


where 


E fV^oV^V^dt. 

CO 27T ^ 


(99) 


Integration of the above equation yields the relation 
E 


E = 
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co 2n i (co - co q ) - y/2 


( 100 ) 


The radiation intensity distribution, I , may now be written as 


I » 
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o 2 tt 
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(co - co^) 2 + y 2 /4 


( 101 ) 


The exact relationship between the electric field and the radiant intensity is 
expressed by means of the Poynting vector, as discussed in the next section, 
subtitle "Rayleigh Scattering Theory. " 

Equation (101) may be normalized such that 

4-co 

I = f I d . (102) 
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The parameter y in equation (101) represents the frequency width at 
half intensity (sometimes called "natural" half-width) and is written as in 
equation (87) , as 


y 


= Ato 


h 



(103) 


Equation (101) is known as the Lorentz line shape, further discussions 
of which are given in the next two sections on shape and broadening of spectral 
lines and band models and absorption coefficients. 


Shape and Broadening of Spectral Lines 


Before a description of the several types of broadening of spectral lines 
is given, it is necessary to summarize and discuss briefly the monochromatic 
energy transitions between discrete energy levels of a molecule. The total 
energy, E, of a molecule may be given according to Baker et al. [71] as 

e = e t +e e +e r +e v + e es +e ns +e ie ■ < 104 > 


where E ^ represents the translational energy; E^ is the sum of the orbital 
energies of the electrons; E-^ is the sum of the rotational energies; Ey is the 
sum of the vibrational energies; Egg is the sum of the electron-spin energies; 
Ejsjg is the sum of the nuclear-spin energies; and Ejg is the sum of the inter- 
action energies among these modes. 

Since translational energies in equation (104) occur in a continuum, any 
value is possible. However, all of the other energies are quantized and must 
be specified by quantum numbers. For a polyelectronic atom, these quantum 
numbers are usually designated as the principal quantum number, n, (measures 
extent of electron M cloud M ) , azimuthal quantum number, i , (specifies electron 
angular momentum) , magnetic quantum number, nig, (gives the component of 
angular momentum in a given direction) , and spin quantum number, m g , (gives 
the component of spin, s, in a given direction) . For a molecule, at least two 
other quantum numbers (J, the rotational quantum number, and v, the vibra- 
tional quantum number) must be considered. 
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To understand adequately the fine structure of the infrared molecular 
spectrum, it is necessary to consider three types of molecular models. The 
first of these is the anliarmonic oscillator which is characterized by a parabolic 
potential curve. The second is the nonrigid rotator which represents a rotating 
system composed of two mass points connected by a massless spring. The 
third model is the vibrating rotator (or rotating oscillator) in which it is assumed 
that both rotation and vibration can take place simultaneously. Since the third 
model represents a combination of the other two, only the energy level of this 
vibrating rotator is given. This is 


(„,J) =w .(v + -f) -w x (v + -h 2 + co p y Jv + 4) + ... 


he 


e e 


e ' e 
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(105) 


+ B J(J+ 1) - a (V + -J-) J( J+ !) - D J 2 (J + i) 2 -/3 (v + |)J 2 (J+ 1) 2 +.. 


The terms on the right-hand side of the first line represent the energy 
caused by the anharmonic oscillators; the B e and D e terms on the second line 
represent the energy associated with the nonrigid rotations; and the a e and p e 
terms represent coupling energy associated with the vibrating rotations. The 
quantities x e , y e , B e , D e , a Q , and /3 e are constants defined by Herzberg [62]. 

The analysis of energy quantization as applied to atoms and molecules, 
of course, belongs in the realm of atomic and molecular spectroscopy, the 
details of which will not be given here. Excellent references in this area are 
the three books by Herzberg [62, 63, 64], and the books by Penner [65] , 

Griem [66], Bauman [67], Barrow [68], Bond et al. [69], and Harrison et al. 
[70], 


As the temperature of a molecule rises, the translational energy is first 
increased; then the rotational states, the vibrational states, and finally the elec- 
tronic states become excited. The quantized energy changes correspond to the 
theory of Bohr as stated by equation (7) (repeated for convenience) as 

E -E = h y . 
n m nm 

Since the rotational energy change, E n - E m , is comparatively small 
(of the order of hundredths of an electron volt) , the frequency i^ nm is small; 
hence this type of energy transition occurs only in the far infrared. Vibrational 
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energy changes (and associated frequencies) are slightly larger (of the order 
of tenths of an electron volt) ; therefore, this type of transition occurs at a 
slightly shorter wavelength (near infrared). Electronic energy changes (and 
frequencies) are the highest; hence, this wavelength is the shortest, and these 
changes occur in the visible and ultraviolet region. These three types of energy 
level transitions may be seen in Figure 4, which shows a typical energy level 
diagram for upper and lower excited electronic states in a diatomic molecule. 



FIGURE 4. ENERGY LEVEL TRANSITIONS FOR A DIATOMIC MOLECULE 

For atoms in a monatomic gas (uncombined atoms) , the vibrational and 
rotational energy states are absent, and the spectrum is merely a sharp, dis- 
continuous series of individual lines. For molecules in a more dense gas (such 
as a rocket exhaust) , the spectrum is made up of a nearly infinite number of 
very closely spaced lines (or bands) . This band emission of polyatomic gases 
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is concentrated mostly in the infrared part of the spectrum, representing vibra- 
tional-rotational and purely rotational transitions. The radiative intensity from 
rocket exhausts in the infrared is consequently much higher than the intensity 
observed in the visible and ultraviolet regions. 

Since it is known that gaseous emission and absorption do not take place 
at a single frequency but occur in a broadened portion of the spectrum, it is 
necessary to establish the factors which bring about this broadening phenomenon 
and to discuss briefly some of the types of broadening. In general, there are 
three major types of broadening of spectral lines: (1) natural-line broadening 
(the finite lifetime of initial and final states caused by spontaneous emission) , 

(2) Doppler broadening (random Doppler shift as a result of the thermal motion 
of the molecules), and (3) pressure broadening (perturbations or displacement 
of energy levels during collisions of the molecules with adjacent molecules) . 
Each of the above processes is discussed below based upon the references of 
Marganeau [72,73], Marganeau and Watson [74] , Marganeau and Lewis [75], 
Ch'en and Takeo [76], Breene [77,78], Aller [79], Thompson [80], Mitchell 
and Zemansky [81], Seshadri and Jones [82], Garbuny [43], and Benedict et al. 
[83], 


1. Natural -Line Broadening. A spectral line emitted by an atom is not 
infinitely sharp because of the phenomenon of radiation damping. (See section 
entitled "Gaseous Radiation," subsection "Radiation from an Accelerated 
Charge.") By this damping process, an oscillating, and hence radiating, elec- 
tronic charge is continuously losing radiant energy, and its amplitude of oscilla- 
tion decreases as the charge continues to oscillate about its natural frequency, 
co G . However, the damped oscillation is not monochromatic, as shown by the 
Lorentz profile of intensity, equation (101) (repeated here for convenience) as 
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co 


= I 


o 27 r 


1 

(co - C0 Q ) 2 + y 2 /4 


(101) 


The quantity, y, which equals Aco h , the natural half- width, is also 
equal to the reciprocal of the time constant, r. This equality comes from 
Heisenberg's uncertainty principle in quantum theory 


AEAt = RAcot 


_h_ 

27T 


Acot 


27T 


(106) 
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as applied to the transition from an excited (upper) level E n to the ground 
state Ej. In this manner there is an uncertainty of energy AE (and hence of 
frequency Aw) associated with a finite lifetime t of the upper level. The life- 
time r is very large for the ground state of an atom and very small (~10 -8 sec) 
for an excited state. The energy level corresponding to the ground state Ej 
will be infinitely sharp, but the excited states will have a natural width. The 
actual lifetime in the excited state will not be any longer than the natural life- 
time, and according to equation (106) , the half-width of the observed spectral 
line cannot be any less than its natural width Aw. 

It must also be mentioned that, if a transition occurs from a higher level 
with a lifetime r 3 to a lower level not in a ground state, but one with a lifetime 
T 2 , the resulting uncertainties must add up to 


AE = AE 3 + AE 2 


(107) 


In this case Ej would represent the ground state. Thus, the resulting half- 
width is 


. _ __J_ ^_1 AAA 

Aw 32 - y 32 - + - A 32 + A 31 + A. 

t 3 t 2 


21 > 


(108) 


where the A nm are the Einstein transition probabilities originally discussed in 
the section entitled "Equation of Radiation Transfer. " The natural width v 32 is 
thus affected by the lifetime of the upper state as well as that of the lower state. 

2. Doppler Broadening. Doppler broadening is a major cause of broad- 
ening in rocket exhaust gases (especially at very low pressures) and is associ- 
ated with the thermal motion of the atoms in the gas. When an oscillating 
charge at frequency u Q moves with velocity v x in a direction x to the observer, 
the emitted wavetrains arrive at a higher frequency v. This frequency, by 
virtue of Doppler's principle, is 

v 

V = V o {i ~~T ) ' (109) 
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From equilibrium statistical mechanics, the fractional number of mole- 
cules dn/n moving with linear velocities between v x and v x + dv x may be 
written in a Maxwell-Boltzmann distribution as 


dn 

n 


m 


27rkT ; 
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/ -mv^ 

exp \^F 




( 110 ) 


where n is the total number of molecules per volume, and m equals the mass 
per molecule. The intensity distribution, l p , within the spectral line is pro- 
portional to the number of molecules radiating at a frequency, v, and can be 
written as 


I dr = I exp 
v o 
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( v - v \ 2 

-me 

1 0 1 
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dr = I exp[-i/>(r - r Q ) 2 ] d" » (HI) 


where 


I = 


me 


27rkTr 2 


( 112 ) 


This Doppler broadening of line intensity thus follows a Gaussian distri- 
bution of frequencies of the form exp[-i/)(r - r G ) 2 ] . The Doppler width is thus 
determined as the increment in frequency for which I = I Q /2 and is given by 


Ar h = 2r o [(2 In 2) ^ ] 


(113) 


It might be mentioned that the Doppler shift principle has been used, together 
with a laser beam, to measure velocities of particles in wind tunnels. Discuss- 
ions of this technique are reported by Robertson et al. [84], Foreman et al. 

[85], and James, Seifert and Babcock [86]. 

3 pressure Broadening . Pressure broadening includes several types of 
broadening which involve interactions between randomly emitting atoms and mole- 
cules, the interactions in each type all increasing with increased pressure. The 
two main types of pressure broadening discussed in this section are collision 
broadening (including resonance broadening) and Stark broadening. Two general 
theories, the statistical theory and the impact (or interruption) theory, are 
associated with pressure broadening and are briefly discussed in this section. 


39 


a. Collision broadening. Collision broadening represents the 
shortening of an oscillator lifetime by a collision between some light emitting 
(excited) atom (or molecule) and another atom or "perturber. " In considering 
these interactions between an excited atom and its neighboring atoms, a dis- 
tinction between the type of neighboring atom (or molecule) must be made. The 
surrounding atoms may be either of the same kind as the excited atoms or may 
be of a different kind. 

If the surrounding atoms are of the same kind, the possibility of reso- 
nance (energy of excitation of one being absorbed by the other) may exist. In 
this manner, when an excited atom passes by an unexcited atom of the same 
kind, a finite probability exists that the excitation energy will be transferred 
from the first to the second atom, without the intervention of radiation. This 
consequently results in a reduction of the excited atom's lifetime and a corres- 
ponding broadening of the line emitted. 

.By classical description of this resonance, each atom is usually replaced 
by an oscillator with a single natural frequency, u Q , and a coupling force is 
then introduced between two dipoles. This coupling force has the form 


F.. 

i] 


Const 
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(114) 


where r*j is the distance between the i th and jth oscillator. Since, under the 
same interaction force F y , two unlike atoms result in an inverse 6 th power 
dependence (instead of an inverse 4 th power dependence for the like atoms) on 
the atomic separation, greater broadening will result from interactions with 
like atoms (or molecules) . However, it is nearly always assumed that the 
number of foreign (different) atoms or molecules surrounding an excited atom 
is far greater than the number of like atoms surrounding the excited atom, so 
that resonance broadening can usually be neglected. 


In the rest of this discussion on collision broadening, it will be assumed 
that the broadening is caused by unlike atoms (foreign gas) interacting with the 
excited atoms. The original analysis of this type of broadening was discussed 
by Lorentz [87] . He assumed that, while an atom is being forced to vibrate 
by the incident light wave, its vibrations are suddenly stopped as the atom 
collides with the other atoms and the energy becomes wholly kinetic. This type 
of process will reduce the total intensity of the line in addition to broadening its 
shape. The spectrum of the electrical oscillations is given by the same relation 
for radiation damping, except for the increased factor, r, as 
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The parameter r is defined as 


r T TR + r COLL 


(116) 


where Y TR is the spontaneous transition probability, and r COI T is a collision 
frequency given by kinetic theory as COLL collision 


w = 6 - 8X1010 


P(T) 


(117) 


fnH m dl ^ ^ ai T 016 ° PtiCal collision diameters of the molecules in A- M, 
and M 2 are molecular weights in g/mole; P is the partial pressure in atm- and 
T is the absolute temperature in degrees K. ’ 

For similar intensities and line widths, the collision (Lorentz shape) 

(plotted versus v) in thf wtgs ' 
u , ,® D ppler (Gaussian sha P e ) curve. Also, it is known that at high alti- 
(low pressures) and at high temperatures, the Doppler half-width is 
considerably higher than the collision half-width. Ritland [55] presents a 

“‘ C / f n h0W f C< ^ «» <*>" a, sea lelel’a“ 3 ° K. 

fnn non f f + OP f ler half ' width to collision half-width is only 1/30, while at 
200, 000 feet and at 2000° K, this ratio increases to 360 For H 2 0 in the 2 7-u 

mfo isSvr^vf T K> the D ° PPler half ' Width 1 COllisi0n half -width 

600 The i L ’ "fn* * ’ ° 00 feet at 200 ° S K this ratio increases to 

Figures V* TIT °l °° Ppler half “ width to lorentz half-width are plotted in 
Figures 5 and 6 for H 2 0 and C0 2 , respectively. 

caused bv thA S ! ark br ,° adening - Stark broadening is pressure broadening 

Z l atoms Tht i P f °! ariZ f i0n of radiation ^ electric fields of neigfbor- 

ionized rocket ZhT T br ° ademng ’ whlch is ver y important in plasmas and 
ionized rocket exhausts, was originally investigated by Stark [86]. Stark deter- 

"elds lines ™ SPlit “ the P~ of citric 

, that excited atoms and their neighbors exert a distribution of electric 

TtZTT ° n V 1 ^ er ' He f0lUld Sin0e * electric field w«h 

strength, E, would split a spectral line into its so-called Stark components, a 
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FIGURE 5. RATIO OF DOPPLER TO LORENTZ HALF-WIDTHS AS A 
FUNCTION OF ALTITUDE OR PRESSURE FOR H 2 0 IN 2. 7-jli BAND 
AT VARIOUS TEMPERATURES 


non-uniform field would produce a whole series of Stark components, and thus 
result in broadening of the line. 

Holtzmark [87] , on the basis of the Stark effect, developed a theory of 
line broadening which assumed that the radiating atoms are in the electric fields 
of ions, dipoles, or quadruples. On the basis of this theory, when a spectral 
line is under the influence of a given field, the line is perturbed in some pre- 
dictable manner. That is, it may be split (linear Stark effect) or shifted 
(quadratic Stark effect) . 
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FIGURE 6. RATIO OF DOPPLER TO LORENTZ HALF-WIDTHS AS A 
FUNCTION OF ALTITUDE OR PRESSURE FOR C0 2 IN 4. 3-p BAND 
AT VARIOUS TEMPERATURES 

The linear Stark effect vanishes for all atoms except the hydrogen atom 
For the hydrogen atom a line under the effect of an electric field, E, splits 
into two components with a frequency interval Ar gT of 


A "ST 87r 2 me \ E aE 


(118) 
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where is a quantum number and a is the split width per unit strength. When 
a large number of atoms which have an average field, Ep)> from the surrounding 
dipoles is involved, the broadened line has a half-width of 


Al, h = 2 IT = aE D • 


(119) 


The line shape is similar to that in Lorentz-type broadening: 


y I 21 „ 
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The half-width is thus determined as 


A iv = aE = 4. 54 ann , 
h D 


( 121 ) 


where g is the dipole moment and n is the number concentration of the mole- 
cules. 

For molecules having a quadrupole moment q in a field of strength, E^, 
the Holtsmark theory predicts the half-width to be 


Av = 0. 67 aE =5. 52 qn a. 
h q 


( 122 ) 


For ions of number density n^, in a field of strength, Ej, the 
Holtsmark theory gives for the half-width 


Au = 1.25 aE. = 3.25 en. a. 
h i i 


(123) 


The quadratic Stark effect appears in non-hydrogen-like spectra, and it 
involves a shift of the line (rather than a split) which is determined as 


6v = v - v = bE* 1 
o 


(124) 
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The quadratic Stark effect may be analyzed classically if it is assumed 
that the application of the electric field induces a polarization, P, proportional 
to the electrical field. In this case the polarization is 

P = ’ (125) 


where (3 is the polarizability of the atom. The work performed in producing an 
increment of polarization dP is therefore 

E 1 

AW=//?EdP= — /IE 2 , (126) 

which implies a quadratic dependence on the electric field, E. 

c. Statistical theory. The statistical theory of broadening is used 
to analyze the broadening phenomenon associated with an emitting atom or mole- 
cule which is placed in the midst of an aggregation of broadeners under static 
conditions. This theory assumes that the atoms or molecules are stationary 
and account is taken of only the static distribution of energy levels during en- 
counters. The broadening by coupling (Stark effect) results from the statistical 
distribution of distances which various neighboring atoms may assume with re- 
spect to the emitting atom during times presumed to be long in comparison with 
that of the emission process. 

In analyzing the statistical theory, the Franck-Condon potential curves 
should be considered. These curves merely represent two electronic energy 
states of a radiating atom plotted versus distance r from another atom. Where 
the atoms are at large distances from each other, the two energy levels are 
undisturbed; however, when the atoms are moved closer together, a strong 
repelling force is exerted on each atom. In general, the energy difference AE 
between the two states is not constant so that the emitted frequency depends on 
the instantaneous position, r, of the excited atom. If the atom (or molecule) 
radiates when die separation is r, the net change in energy, AE, is 


AE - hr = E 2 (r) - E t (r) . (127) 


This energy differs from that ordinarily radiated since E 2 - Ej will 
change with r as a result of perturbations caused by other molecules, and the 
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position of the spectral line will be shifted. If a statistical theory is used which 
involves a van der Waals force interaction, the static frequency shift is 


Av 


AE 

h 



(128) 


where /3 is a constant to be determined. 

According to Breene [77], an expression may be written for the prob- 
ability that the broadening atoms (or molecules) will be so distributed in space 
as to distort the natural frequency to a particular value. The intensity in the 
resulting spectral line will be proportional to the desired frequency and may be 
written as 


I(A^) = y(Av) 


-3/2 


Q -Try/Av , 


(129) 


where 




(130) 


Here N is the density of the perturbers and j3 is the van der Waals constant 
in equation ( 128) . 

The half-width may then be written as 


Ar, = 0. 82 7r -y 2 . 
h 


(131) 


The shift of the intensity maximum for van der Waals forces is then determined 
as 



(132) 


46 


d. Impact theory. The impact (or interruption) theory is used 
when the velocities of the perturbing atoms or molecules are so large that the 
times of interaction are small compared to the duration of the emission process. 
This theory leads to the equation for Lorentz or collision broadening. In this 
case it is assumed that an atom or molecule is absorbing or emitting at a sharp 
frequency Vq during the time between collisions. Then at the moment of colli- 
sion the radiation process stops completely, and the energy of vibration is con- 
verted pui’ely into kinetic energy. 

A detuning ot natural frequency co 0 is thus experienced as the perturbing 
particle approaches the emitting atom with a uniform velocity, v. Resulting is 

a net phase shift uoj which depends on the inverse power of the separation dis- 
tance r as 



( 133 ) 


where C is a constant. The total phase shift, r\, is thus found by integrating 
over the duration of the encounter as 
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where a n is a function of gamma functions and p is a "closest approach" 
parameter. 

The intensity distribution of a spectral line may be obtained as a function 
of the time between encounters, t, and mean time between encounters, t 0 , 
according to Aller [79] as ’ 
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where A is determined by the condition that 
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According to Wooly and Stibbs [21], the half -width may be obtained from 
impact theory as 


A, h = 
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= — p 2 v N , 

o > 


(137) 


where 7 C is the collision frequency defined by equation (116) ; N is the number 
of particles per cm 3 responsible for the broadening; and p G represents the dis- 
tance of closest encounter. This distance is given by 


P 


o 
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(138) 


when C is a constant and v is the mean relative speed of the radiating and per- 
turbing particles given by kinetic theory as 



Also, M 0 is the mass of a particle of unit atomic weight, and A t and A 2 are 
the atomic weights of the radiating and perturbing atoms. 

This concludes the section on the shape and broadening of spectral lines. 
Much more information on these phenomena may be obtained by consulting the 
works of other authors [21, 43, and 73-89] . In the next section the various 
band models are discussed. These models are important in the calculation of 
absorption coefficients, which are, in turn, necessary for predicting radiation 
heating from rocket exhausts. 


Band Models 


To predict accurately the gaseous radiation from a rocket exhaust, the 
spectral absorption coefficient, (or k ^) , as a function of wavelength (or 
frequency) for a particular temperature, pressure, and concentration must be 
known. A considerable number of experimental measurements have been made 
of absorptance and transmittance of rocket exhaust gases in the last several 
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years. However, it is not possible to obtain enough experimental data for all 
possible variations in temperature, pressure, concentration and frequency 
which might occur across a particular slant path in a rocket exhaust as viewed 
by a radiometer or spectrometer. 

Because of this lack of experimental data, the concept of band models 
has been developed to predict the absorption in a particular band composed of 
a nearly infinite number of spectral lines. Through the use of band models 
(which are assumed to be mathematically correct) , the absorption, and con- 
sequently the emission of a rocket exhaust may be determined as a function of 
path length f and pressure of the gas, P. 

According to Jamieson et al. [90] , the absorption of a gas represented 
by band models depends on the number of lines in the band, on the relative 
spacing of the lines with respect to each other, on the number of spectral lines 
with a given value of the intensity, on the half-width of the spectral lines, and 
on the line shape (particularly in the wings of the lines) . Theoretically, the 
band model is considered to have an infinite number of absorption lines of 
uniform statistical properties. An interval of this particular theoretical type 
containing many lines is assumed to have properties similar to those of an actu- 
al interval of the band being analyzed. In the band model, each interval will be 
flanked by statistically similar intervals; however, in a real band this is not 
necessarily the case. 

In this section, five particular band models ( including the single-line 
model) are discussed: (1) single-line model (exact method for absorption by 
isolated lines), (2) Elsasser model (band composed of identical, uniformly 
spaced lines) , (3) statistical model (band composed of spectral lines with 
arbitrary intensity distribution and with random spacing between lines) , (4) 
random Elsasser model (band composed of several groups of lines, each of 
which forms an Elsasser band, but the groups are superposed with random 
spacing) , and (5) quasi-random model (band composed of spectral lines whose 
intensities and variation of spacing from line to line are accurately simulated) . 

There have been numerous references on the general subject of band 
models in the last 25 years. The references most referred to in this section 
include the articles by Plass [91-96], Openheim and Ben-Aryeh [97-98], 

Kaplan [99,100], Simmons [101], Penzias and Mac lay [102], Maclay and 
Brabov [103], Green and Wyatt [104], Wyatt, Stull, and Plass [105,106], 
Elsasser [107], Godson [108], Greif [109], Goody [25], and Yossa [50]. 
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1. Single-Line Model. The simplest band model to be discussed is the 
single-line model of Lorentz shape in which the lines are so narrow that there 
is essentially no overlapping of lines. The overlapping or nonoverlapping refers 
to the extent to which absorption in some small spectral interval is caused by 
contributions from a number of lines as opposed to absorption by a particular 
well-isolated line. When the lines do overlap, each line absorbs a smaller 
percentage of the radiation than it does when it is isolated from the other lines. 

A fundamental quantity to be used in band model theory is known as the 
spectral absorptance, A^, equal to the spectral emissivity, when LTE is 
assumed, and defined as 


A„ 




(140) 


where i is the path length through the absorbing gas. An analogous expression 
for Aj, was originally introduced in the section entitled "Equation of Radiation 
Transfer" as the spectral absorptivity. In band model theory this is usually 
written in terms of frequency v and termed absorptance; therefore, the notation 
of equation (140) will be used throughout the rest of this section. 

The spectral absorption coefficient, k v , may be written for a Lorentz 
line shape (the most common shape occurring in actual practice) as 
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where S is known as the line strength and is defined by 


(141) 



(142) 


The half-width y is defined in a manner similar to equation (116) as 


y = y° p + y" P. 
a a b b 


(143) 
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where Y a represents the half-width at unit pressure for self-broadening (colli- 
sions between the absorbing molecules) , and P a is the pressure of the absorb- 
ing gas. The quantity represents the half-width at unit pressure for foreign 
gas broadening (collision between absorbing molecules and other molecules in 
the exhaust gas which are not absorbing radiation at a frequency u ) . The 
quantity represents the pressure of the nonabsorbing line -broadening gas. 

The line strength S is also related to the pressure of the absorbing gas, 
expressed as 


S = S° P , 
a 


(144) 


where S° is the strength at unit pressure. 

The total or integrated absorptance, W, (also known as equivalent 
width) of a single line may be specified at a particular temperature by means 
of the line strength, S, • the line half-width, y, and the path length, i. This 
integrated absorptance (or equivalent width) may be written as 

W= /A dv = J (1 - e~ Kpl ) dv . (145) 

o o 

If equation (141) is inserted into equation (145) and the result integrated 
from zero to infinity, it can be shown, as it is by Ladenberg and Reiche [110], 
that 


W = 27ryf(x) 


(146) 


where 


Si 

X ~ 27ry ’ 

and f (x) may be expressed in terms of Bessel functions as 


(147) 
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(148) 


f (x) = xe X [ J Q (ix) - iJj(ix) ] . 


Two useful asymptotic forms occur for equation (148) . When the param- 
eter x is very small (linear approximation) , f(x) -*x, and the integrated ab- 
sorptance can be written as 


W = 27ryx = Si = S° P l , 

a 


(149) 


which represents a linear variation of W 
Another important parameter, /3, 


as a function of P . 

a 

in band model theory is defined as 


0 = 


2tTY 

d 


(150) 


where d represents the average spacing between the spectral lines. 

The integrated absorptance may thus be written as 

W = /3xd or = fix . 

d 

1 . 

For large values of x (square root approximation) , f(x) — (2x/ir) 2 and 

i i 

W = 2y(2nx) 2 = 2(S4y) 2 . 


(151) 


(152) 


The integrated absorptance may also be written in terms of P as 

1 

W = d(2/3 2 x/7r) 2 , (153) 


and if it is assumed that 7 =7°P a and S = S°P a , equation (153) may be written 
as 
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SQUARE ROOT APPROXIMATION 


(154) 


W = 2(y°S°P 2 i)2 or w° = 2(y 0 S°i) 4 , 
a 


which is known as the square root approximation for the integrated absorptance. 
This square root approximation, together with the linear approximation, is 
shown in Figure 7 in comparison with the Ladenberg and Reiche equation (148) . 

The quantities y° and S° may be determined theoretically as by Howard 
et al. [Ill, 112] or experimentally as by Benedict et al. [83]. Since the above 
theoretical method can predict integrated absorptance readily, provided suffi- 
cient spectroscopic data are available, it is called the exact theoretical method. 
The above equations are valid for absorption by a band of spectral lines, pro- 
vided the lines do not overlap. This assumption holds true when both the pres- 
sure and path length are small. At high altitudes (or low pressures) , the 
Doppler contribution to the line shape can be important (as seen in Figures 5 
and 6) . Plass and Fivel [113] and Struve and Elvey [114] discuss this contri- 
bution and suggest various logarithmic expressions for the integrated absorp- 
tance (or equivalent width) . In the next four sections on band models, various 
approximate methods of determining absorptance which may involve overlapping 
of the spectral lines are discussed. 

2. Elsasser Model . The Elsasser band model, which was the earliest 
band model investigated [115], assumes an infinite array of spectral lines of 
equal intensity which are regularly spaced a distance d apart. This configu- 
ration is represented by the sketch in Figure 8. 



FIGURE 8. ELSASSER MODEL 


Figure 8, compared with Figure 9, shows that there is, greater absorp- 
tion in this model in which the spectral lines are evenly spaced than when the 
line spacing is nearly random. Some portions of the C0 2 spectrum may be 
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represented with a fair amount of accuracy with the Elsasser model; however, in 
the actual CO 2 spectrum there are always weaker lines between the strong, 
regular-spaced lines in the Elsasser model. 

The Elsasser model uses the Lorentz line shape, equations (101) and 
(141) , with a slight modification. The absorption coefficient, « is summed 
over all the spectral lines in the band, and thus represents the total contribution 
to Kj, from all the lines as 


K 

v 


n= + 0 ° 

l 


11= -00 


S 

7T 


Y 

( v - nd) 2 + y 2 


(155) 


It was shown by Elsasser (115) that this summation may be written in a simpli- 
fied form as 


K 

V 


S_ sinh /3 
d cosh /3 - cos £, 


(156) 


where 


- 27riVd . 


With the insertion of equation (156) into the equation for the absorptance, 
equation (140) , and using Elsasser's expression for average absorptance at the * 
band center, 


A 

v 


. +00 


+00 


/ d u/J du 


(158) 


where 


dv = (d/27r) d£ 


(159) 
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the average absorptance may be written as 


A 

v 



, +7r r si 

sinh P 

i exp [- d 

-7T 

cosh p - cos £ _ 


(160) 


Two approximations of equation (160) exist: the so-called "strong line 
approximation and the "weak-line” approximation. The strong-line approxima- 
tion is valid when the spectral lines are far apart compared to their half-width 
(d » y, or /3 — »0) . In this case the absorption is practically complete near the 
centers of the strongest lines in the band. The expression for the absorptance 
thus becomes 


A 

v 


1 Z _ z 2 

erf (z) =2/(7 r) 2 J e dz , 
0 


(161) 


where 


z = (i/3 2 x) 2 



(162) 


The weak-line approximation is valid when the spectral lines are close 
together (d— 0, or 0— °°) . In this case the absorption is small at all frequencies 
in the band even including those near the center of the strong absorption lines. 

The expression for the absorptance then becomes 


A 

v 



r +7r -i S/d 

J 6 


d£ 



(163) 


This expression is analogous to equation (140) and is known as the Beer Lambert 
equation. 

3 statistical Model . The statistical or random band model, originally 
developed by Mayer [116] and Goody [117] , assumes a band in which there is 
no correlation between line intensities and frequencies, and in which there is a 
random spacing of lines. The intensity of the spectral lines may vary in any 
manner as long as it can be represented by a particular distribution function. 
Also, as a result of different quantum transitions, series of spectral lines may 
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overlap in several bands. It is known that the absorption of H 2 0 may be repre- 
sented by the statistical band model over a moderate range of pressures and 
path lengths. Figure 9 shows a sketch of the statistical band model showing 
regions of overlapping lines. 



V 


FIGURE 9. STATISTICAL MODEL 


The basis of the statistical band model is the determination of the most 
probable average transmittance at the center of a frequency interval as a result 

of the absorption by n lines of intensity Sj (i = 1, 2, 3 n) . In this theory 

various probability distributions are introduced in the following manner. If 
N(y) is the probability that a spectral line lies between v and v + dv, and 
diq, ^2 is the probability that one line lies between and + diq 
and that a second line lies between v% and + dy 2 , then the probability of find- 
ing n lines in some frequency interval D is unity, or 


/ N(^ 1 ,j^ 2 ,. 


D 


V ) diq, d^2, .. .dp =1 

n 1 n 


(164) 


Also, on the same basis, the probability that a given spectral line has 
an intensity between Sj and + dSj is P(S^) . This probability function may 
usually be expressed in terms of the Dirac delta function, <5(S - S) , or an 
exponential function, (1/S) exp(-S/S), where S is some mean line intensity, 
described by Malkmus [118] . Normalizing this probability results in the 
relation 


/ P(S.) d(S.) = 1 . (165) 


The weighted average of the transmittance, t u , over the frequency 
interval, D, may now be written in terms of the above probability functions as 
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(166) 


The random distribution hypothesis is now introduced to simplify the 
above equation. This hypothesis assumes that it is equally probable that each 
spectral line has its center at a given frequency v in the frequency interval D 
of n lines, regardless of the position of the other lines. In this manner the 
expression N(iq, • • • • ^ n ) equals a constant; i. e. , it is independent of fre- 
quency, and may be written as 


N(p 1 ,p 2 , • • • » ^ n ) = D 


-n 


(167) 


Since the integral over each line is equal to the integral over any other 
line, equation (166) reduces to 


t = D 

v 


-n 


(D 


f dv f P(S) exp(Si) dS 


-|D 


n 


(168) 


Since the average absorptance, A p , in the case of a nonscattering gas 
(reflectance equals zero) is equal to 1 - T pt the following expression for A„, 
following the mathematical procedure of Plass [91] , may be written as 


A = 1 - 

v 


1 - D 1 f W, 
0 


S.L. 


(S) P(S) dS 


(169) 


If the average value of the equivalent width of a single line is given as 


W S.L. =/ W S.L. (S) P(S) dS ’ 


(170) 
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the expression for the average absorptance finally reduces to 


A - 1 - exp[-(W /d> ] . 

V o. 


(171) 


Since Wg ^ is the same as W inequation (146), a strong-line approxi- 
mation results in which f(x) — (2x/tt) 2, and the average absorptance may be 
written as 


A = 1 - exp 
v 


[-(2/3 2 x/7r) 2 ] . 


(172) 


A weak-line approximation results when f(x) — x as 


A 

v 



(173) 


which has the same form as the weak-line approximation in the Elsasser band 
model, equation (163). 

The regions of validity for the statistical model, as well as the Elsasser 
model, for various band absorption approximations such as strong-line, weak- 
line, and non-overlapping line approximations are given by Plass [94, 95]. 

4. Random Elsasser Model. The random Elsasser model, in general, 
provides a more accurate representation of spectral bands than either the 
Elsasser or statistical model. The random Elsasser band model assumes that 
the absorption may be characterized by the random superposition of Elsasser 
bands, each of which may have a different line intensity and spacing between 
the lines. Through the use of this model, the absorption of HC1 has been pre- 
dicted quite accurately by Stull and Plass [119] and Malkmus et al. [120]. 

Figure 10 shows a sketch of a random Elsasser model. Four different 
Elsasser bands, each with a different intensity and line spacing, are super- 
imposed on the figure. In this manner, as many of the weak spectral lines may 
be used as desired to make a satisfactory prediction of absorption for the par- 
ticular pressures and path length considered in the rocket exhaust. 
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V 


FIGURE 10. RANDOM ELSASSER MODEL 


For this model it is assumed that d i represents the spectral line spacing 
of the i^ 1 Elsasser band and n is the number of superimposed Elsasser bands. 
The average spacing 6 between the spectral lines is thus 



(173) 


An equivalent width for the i 


th 


Elsasser band may be written as 


W„ .(x.,/3,.) =d. 
E,i l di i 









(174) 


where /3 is the parameter defined by equation (150) , x is defined by equation 
(147), and tg ^ and Ag ^ are the average transmittance and absorptance, re- 
spectively, of ’the ito Elsasser band. By a procedure similar to that in the 
preceding section, the average absorptance from a random superposition of n 
Elsasser bands can be written as 


A 

v 


= 1 - 


n 

n 

i=l 







P £ (S.) dS. , 


(175) 


where Pg(S^) dSj represents the probability of the i^ Elsasser band having 
an intensity in the range dS^. 
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An average equivalent width W £ { may be defined in the same manner 
as equation (170) . The average absorptance from a random superposition of 
n Elsasser bands, each of which may have a different intensity, half-width, and 
line spacing, can be written as 



If the line spacings and half-widths are all assumed equal (d = n <5) , the 
average absorptance may be written as 


A = 1 

v 


'W 


1 - 


EJ_ 


n 6 


n 


(177) 


In terms of the quantities /3 and x (or of S and P) , the strong-line 
approximation for the Elsasser model is 


— n 

A = 1 - n i 1 - erf 
v i=i 


[d<to 4 ] 


(178) 


and the weak-line approximation is 


A = l - 

V 


(179) 


5 - Quasi-Random Model . The quasi-random band model is possibly the 
most accurate model of the five mentioned in this report. In an actual band the 
spectral lines are not as regularly arranged in frequency or intensity as in the 
Elsasser or random Elsasser models, nor are they arranged in as random a 
manner as in the statistical model. The quasi-random model is thus better 
able to simulate the actual intensity distribution (including many weak lines) 
and the relative spacing of all the spectral lines. Furthermore, this model is 
able to take into account the effect of the absorptances of the adjacent frequency 
intervals on the particular frequency interval being analyzed. However, because 
of certain complexities in its usage, this model has not yet been programmed 
for rocket exhaust radiation calculations. 
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Using the quasi-random model, the absorptance is first calculated over 
frequency intervals 6 which are smaller than the particular interval size of 
interest. These intervals are assumed so small that the spectral lines lying 
in any one of them may be considered to lie at random positions, without any 
appreciable errors being introduced. This confines the stronger lines to a 
narrow interval around their actual position and prevents overlapping. Theoret- 
ically, if the interval spacing were decreased to zero, this band model would 
locate each line exactly, with no approximations needed at all. However, a 
large enough spacing is, of course, taken to simplify calculations. Also, the 
lines in each frequency interval 6 are divided into intensity subgroups according 
to the respective magnitudes of their line strengths, 

The average absorptance at the frequency v resulting from n lines in 
the intervals 6 k with line centers at the frequencies i^(i = 1,2, . . . ,n k ) may be 
calculated similarly to the statistical model as 




= 1 - 


1 

6 n k 



(180) 


According to Wyatt et al. [105] , it has been determined that only the 
first five intensity decades for each frequency contributes appreciably to the 
absorptance (or transmittance). Hence, equation (180) may be written as 



i 

6 




(181) 


where is the average value of the line intensity, and the total number of lines 
n k in the frequency interval can be written as 



(182) 


To calculate the absorptance at a frequency v which is influenced by the 
wings whose centers are outside of the interval which contains the frequency v , 
the following equation is used: 


A = 1 - n t. , ( 183 ) 

r j=i W 
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where tj„ is the transmittance at v as a result of nj lines in the frequency 
interval 6j. It may be restated that equation (181) represents the total absorp- 
tance at v for both direct and wing contributions. 

According to Jamieson et al. [90], the strong-line approximation for the 
absorptance using the quasi-random model can be written as 

A -1- n !e”^i - [1 - erf (4 ) ] [ , ( 184 ) 

v i=l / 1 1 ) 

where 

t 2 _ 8y^ x./6? . (185) 

s i Aik 

Equation (184) represents only the absorptance from the spectral lines 
whose centers lie within the small interval 6^. The concentration from the 
wings of the lines in adjacent intervals must be added by using equation (183) . 

Several other band models and approximations have appeared in the 
literature. One of the approximations, the Curtis-Godson approximation, is 
discussed in the next section, as well as certain apparent emissivity relations. 
Also, various programs for predicting gaseous radiant heating to the base of 
Saturn space vehicles caused by the exhaust of H-l, F-l, and J-2 engines are 
described. Some of these programs are based on the above applications of band 

models. 


Methods of Predicting Gaseous Radiation from Rocket Exhausts 


In this section various methods of predicting the radiation heating to the 
base of large rocket vehicles (of the Saturn class) are described. Three gen- 
eral methods are discussed, each of which was either generated under NASA/ 
Marshall Space Flight Center (MSFC) contract or was an MSFC in-house effort. 

The first method, an MSFC in-house effort described by Heatherly et al. 
[121] and [122], used the total emissivity method of Hottel [123]. The second 
method was a University of California method which involved apparent emissivity 
and mean path length relations and was developed under contract to MSFC’s 
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P&VE Laboratory. This method is summarized by Tien and Abu-Romia [124- 
126] and by Giedt et al. [127] . The third general method was originally pro- 
grammed for rocket exhaust applications at Rocketdyne under P&VE contract 
and is described by deSoto [128-130] . This method, which uses spectral 
absorption coefficient data, was modified to include the Curtis-Godson approxi- 
mation by Warner-Swasey Company and General Dynamics/C onvair and further 
documented by Brown Engineering Company — all under MSFC Aero-Astrodyna- 
mics Laboratory contract. A detailed review of the third method is given bv 
Huffaker [4, 51] . 

i - Total E missivity Method . The total emissivity method of Hottel was 
used to predict the radiation heating from the five J-2 engines of the Saturn V 
S-II stage to the S-II stage thrust structure and heat shield as reported by 
Heatherly and Dash [121] . The complete gas dynamic and radiation computer 
program using this method for the intersection regions of hydrogen-oxygen 
plumes for a multi-engine space vehicle is given by Heatherly et al. [ 122] 

By the total emissivity method, the radiation heating from a non-isother- 
mal, isobar ic gas can be written as 


( l/ A = n I I °‘ B T (de/dX) cos0 dco dX, 
co X 


where the parameter X is the optical beam length which is defined as 


X = P i . 
w 


(187) 


P w 1S the partlal Pressure of the water vapor in the exhaust (H 2 0 assumed to 
be the only emitter in the LH 2 -lox system) , and i is the path length through 
the gas. The quantity o^T 4 is the Planck blackbody function, B^(T). 

A solid angle shape factor, F, may be defined as 


F - ~ J cos 0 dco , 
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where 6 is the angle which the line of sight to the rocket exhaust makes with 
the normal to the area where the radiation heating is to be determined. 

Equation (186) may now be written, with the insertion of equation (188) , 
as 


“ /a = f ‘’b/ t ; 1(F1) d( V> • (189) 

w 

In this program, total emissivity, e, versus optical beam length, P i, 
for various values of temperature, T, was obtained by curve fitting with 
polynomials, the curves as presented by Hottel [123]. This total emissivity 
program, while simple in scope, is limited to systems of one molecular emitter 
(i. e. , H 2 0 in LH 2 -lox systems) . The program does have its advantages in 
that integration over all wavelengths (use of total emissivity rather than spec- 
tral emissivity) is avoided. 

2. Apparent Emissivity and Mean Path Length Method. The apparent 
emissivity method of predicting gaseous radiation is based upon the following 
equation from Eckert and Drake [131] for apparent spectral emissivity, e^, of 
a semi-infinite cylindrical gaseous body. 


e = 
A 


_1_ 

IT 


//d- 

P <P 


~K\S 

e A ) 


sin/3 cos/3 d/3 dcp 


(190) 


hi this equation, /3 represents the polar angle and <p the azimuthal angle 
as seen in Figure 11. The quantity s is the path length (analogous to i in 
equation (187) . Equation (190) may be nondimensionalized by using the follow- 
ing parameters: 


H = 



R = , K = k r , and S = — 

r A A e r 

e e 


(191) 


where r_ is the radius of the rocket nozzle exit, h is the height of the nozzle 
exit as measured from the base plane, r is the radial distance in the base plane, 
and is the spectral absorption coefficient. 
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FIGURE 11. SKETCH OF SEMI-INFINITE CYLINDRICAL GAS BODY 
USED IN APPARENT EMISSIVITY AND MEAN PATH 
LENGTH METHOD 

Equation (190) may now be written in nondimensional form as 
1 K S 

- e ^ ) sin/3 cos/3 d/3 d (p . (192) 

/3 <p 

This may be further arranged in the form 



(193) 


where 
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F - 


~ J J sin/3 cos/3 d/3 dtp and = 


i -K S 

~ f f e A sin/3 cos/3 d/3 d<£>. 

P <P 

(194) 


The quantity F is the configuration factor or the blackbody apparent 
emissivity (analogous to equation (188)) , and e^ c is the contribution to the 
emissivity caused by the finite absorption coefficient of the gaseous body. Equa- 
tions (192) and (193) may be further modified by making the substitution: 
tan/3 = sec cp tan/?' (see Figure li). In this manner the quantities F and e 
become Ac 


and 

(195) 

d/3' dp . 

Using these relations, curves for the apparent spectral emissivity as a 
function of R for various values of H and have been given by Tien and 
Abu-Romia [132] and [133] for both semi-infinite cylindrical and conical bodies. 
Figure 12 shows the curves for H = 1 as taken from these references. It may 
be seen that the apparent spectral emissivity for a conical body of 10-deg cone 
angle is much higher than that for a cylindrical body for the same value of R. 

The above expressions apply only for a gray-gas assumption. In this 
manner the total nondimensionalized absorption coefficient, K, equals the 
spectral nondimensionalized absorption coefficient, K^, and the total emissiv- 
ity, e, equals the spectral emissivity, e^. However, by making use of the 
mean path length, L, (defined as the path length averaged over all path lengths 
viewed at different angles from a certain location in the base plane) , an expres- 
sion for total emissivity may be obtained without making the gray-gas assump- 
tion. This dimensionless mean path length is a function of H and R and may 
be defined in the expression 


i_ r r cos 2 <j> tan/3 1 sec 2 /3' 

if 3, J (cos 2 <p + tan 2 /3 1 ) 2 


d/3' d</> 


Ac 


P 4> 


-K^S cos 2 p tan/3 1 sec 2 /3 1 
(cos 2 p + tan 2 /3') 2 


e =F{1 - e KaL ) . 
A 


(196) 
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Apparent Emissivity in Base Plane 



FIGURE 12. APPARENT EMISSIVITIES IN THE BASE PLANE OF 
A CYLINDRICAL AND A CONICAL GRAY-GAS BODY FOR H = 1 


The total emissivity may now be written, with the insertion of equation 
(196) , as 


e(H,R) = 

f e B dX 

l x x 

TC t F 

«, [1 - exp(-K L) ] d\ 
r A 

(197) 

00 

J B x dX 
i x 

<V r L 

^ A. 5 [exp(C 2 /AT) -1] 
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ch 
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(198) 
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In these equations, h is Planck's constant and k is Boltzmann's constant. 

If the infrared spectrum is known (experimentally or theoretically) in a 
narrow band width, equation (197) may be approximated by the relation 


e(H,R) [1 - exp(-K^L) ] (D 2i - D^) , (199) 

where the quantities D ni (A ni ,T) are known as the relative cumulative spectral 
radiance, tabulated by Pivonsky [134], and are defined as 

1 ^ n i ^CjdA 

D ni (A ni ,T/ “ a fi T 4 / A>[exp(C 2 /AT) -1] ‘ (200 > 


For a cylindrical semi-infinite body with shielding (H # 0) , Giedt et al. 
[127] and Tien and Abu-Romia [133] suggest the following semi-empirical 
equation for L(H,R): 


h(H,R) 


4 sin/3' 
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tan -1 ( sin /3 1 tan <6 ) 
o o 


E 2 R 1 

( TT 
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/R 2 - 1 
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{ 2 
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’ R / 


TT cos 2 / 3 1 - 
o 

16 R 3 


( 201 ) 


where Ej and E 2 are the elliptic integrals of the first and second kind, respec- 
tively. 


Other expressions for the mean path length based on the experimental 
data of Howard, Burch, and Williams [135] have been given by Tien and Wang 
[136]. For conditions in which a power law type of correlation of integrated 
absorptance can be used, such as 

W = / A p dy = c w a pt> , (0. 5 < a < 1, 0 < b < 0. 5) , (202) 

v 
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where w is the absorber concentration and P e is an equivalent pressure defined 
by Burch, Singleton, and Williams [137], the mean path length may be repre- 
sented as 


L = -jf J r 3 * * & cos cp dci> , (203) 

CO 

where r is the distance between the differential area on the confining surface 
and the differential gas volume. 

For large values of w in which the following logarithmic form of inte- 
grated absorptance may be used: 

W = f A dr = E + F In (wP m ) , m < 1 , (204) 

J V 6 

V 

where E and F are constants defined by Howard et al. [135] , the mean path 
length may be expressed as 

L = exp |"~y J ( In r) cose p dcoj . (205) 

co 

Tien and Wang [136] and Rolfe [138] present ratios of mean path length 
of a non-transparent gas to that of a transparent gas for spherical enclosures, 
for an enclosure composed of two infinite parallel planes, for an infinite circu- 
lar cylindrical enclosure, and foi’ radiation to the center of the base of a finite 
circular cylinder. 

In summary, the general method of apparent emissivity and mean path 
length, while yielding more accurate values of radiative heating than the total 
emissivity method of Hottel, is felt to be inferior to the method described 
below in which band models and/or spectral absorption coefficients are used. 

3. Spectral Absorption Coefficient Methods . The following three meth- 

ods (two of which use band models) of predicting radiation heating from liquid 

propellant rocket exhausts are based upon the use of spectral absorption co- 

efficients. These methods will be termed (a) "weak-line or exact absorption 
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coefficient method," (b) "statistical band model method with Curtis-Godson 
approximation," and (c) "generalized method with modified Curtis-Godson 
approximation. " 

a. Weak-line or exact absorption coefficient method. The weak- 
line or exact absorption coefficient method does not make use of band models, 
since for this method the absorption is assumed not to depend on the spectral 
arrangement of lines in the band. This method sometimes provides more 
conservative values of radiation heating than methods (b) or (c) described 
below. However, one objection to the method is that to obtain a sufficiently 
accurate numerical calculation in certain regions, the wave number (1/A.) 
interval must be made extremely small (~0. 01 cm -1 ) . 

The basic equation to be solved with the weak-line method is equation 
(57) , given in the section entitled "Equation of Radiation Transfer. " The 
integrated form of this equation is equation (61) with 1^(0) , the incident inten- 
sity upon the plume assumed equal to zero. Equation (61) can thus be written 
in the form 


s s 

I (®) = / j, ,(s) exp - f K (s') ds' 
p J Q v t) " 


ds 


(206) 


where the s and s' notation has been reversed from that used in equation ( 61) 
and the subscript v has been used instead of A. 

Assuming that Kirchoff's Law holds, the emission coefficient j p can be 
written as in equation (12) with the refractive index of the medium n^ assumed 
equal to unity. The mass absorption coefficient v will be used instead of 
k v . Equation (206) now becomes 


, °m , s 

I (s ) = f p/c B (T) exp - fpK ds 
p J 0 m,i/ v i 


ds , 


(207) 


where s m is the upper limit of the variable s in the plume along all lines of 
sight beyond which there is zero contribution to the radiation. 

The radiative heat transfer per unit area may now be obtained by inte- 
grating equation (207) over all frequencies v and solid angles dco = sin# d0 d <p 
as 
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q/A = f f f f P K B (T) exp “ / P* ^ s | sin0 cos 6 d0 dcp ds dv. v 
v s <p 9 ’ L 0 ’ J 

(208) 

The absorption coefficients in equation (208) may be obtained experimen- 
tally or from band model calculations such as tabulated by Ferriso et al. [52] . 

Dash [139] describes calculations of radiative heating from the Saturn S-IVB 
stage J-2 engine to the S-IVB thrust structure using equation (208) . The 
numerical solution of equation (208) was discussed by deSoto [128-130] , but 
this program was limited to a single rocket plume with axisymmetric geometry. 
However, Tarbell and Alligood [140] and Alligood [141], using local partial 
pressures and spectral coefficients evaluated at standard conditions in equation 
(208), described a numerical solution to equation (208) which could be used 
for any number of rocket plumes. 

The numerical integration of equation (208) , as discussed by Tarbell 
and Alligood [140] , is performed in the following manner. The integral of the 
absorption coefficient is first treated as 

n 

A(s .,v,0,(f>) = Y p K As. ; (209) 

' 1 r P m,i^ l 

J i=l 

the integral over the line of sight is next evaluated as 
n 

B{v,6,(j>) =Y J exp[-A(s., v,6,4>) ] p« m v B^(T) As. ; (210) 


the integral over frequency v is next obtained as 
n 

c = Yj b (m,$) ; ^ 211 > 

i=l 

the integral over (p is next evaluated as 
n 

D(0) = C(6,<p) sine COS0A0. ; (212) 

i=l 
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and the final integral over 6 is obtained as 


n 

q/A = E = £ D(</>) A (p. . 
i=l 


( 213 ) 


b. Statistical band model method with Curtis-Godson approxima- 
tion. The statistical band model method with the Curtis-Godson approximation 
has been described in detail by Krakow et al. [142] and by Huffaker [4, 51]. 
This method is felt to be superior to the weak-line or exact absorption coeffici- 
ent method described above in that it does not depend upon such a small wave 
number (or frequency) interval for integration. Also, it has its advantages in 
that the band model data such as line strength, S, line spacing, d, and line 
half-width, y, can be obtained from homogeneous samples of gas and then 
applied to an inhomogeneous gas such as a rocket exhaust. 

The equation for spectral intensity used in this method is actually the 
same as deSoto's equation (206) , but has been modified as described by 
Krakow et al. [142] , Krakow [143] and Tourin and Krakow [144] in the follow- 
ing finite-sum approximation: 


n 



where t = exp (~pn m p s) is the transmittance of the section of sample in zone 
i. The quantity v now represents the wave number (1/cm) instead of the fre- 
quency as originally defined. 

The heat flux equation, (208) , now becomes 

m 

q/A = f V I(y.,AiO Av. cosd dco . (215) 

] J 1 


In this method of predicting radiation heating, the statistical band model 
(See this section, subsection entitled "Band Models," part 3. ) with equal line 
strengths and equal line widths was used. Equation (171) was then used for 
the average integrated absorptance. This equation now takes the form 


\ 

\ 
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(216) 


A „ " 1 " 6XP [-< W S. L. /d) ] = 1 - \ = 1 - [- 2*3- f(x) J 


where tj, denotes the average transmittance in the wave number interval Av, 
and f(x) is the Ladenberg-Reiche function in equation (148) . 


If tj, in equation (216) represents the transmittance of an inhomogeneous 
gas, such as a rocket exhaust, then it is possible that some hypothetical homo- 
geneous sample, such as obtained in a laboratory burner, would have the same 
T u if it had certain values of y/d and x. The Curtis-Godson approximation 
[142] as given in equations (217) and (218) is thus used to obtain the band 
model parameters of the hypothetical homogeneous gas sample: 
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(218) 


By dividing equation (218) by equation (217) the following value of 
(y/d) j, which is the value of y/d for the hypothetical homogeneous gas sample, 
is obtained: 
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l (S°/d) 
h=l 
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h h 
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By dividing equation (217) by 27ry/d, the following value of x i# which is 
the value of x for the hypothetical homogeneous gas sample, is obtained: 

[Z,< s ° /d > h p ^ h 

L h=l 


x. = 

l 


2 '| 1 <s ° /d, h p h f h [ <r >h p h + <^ /d >h p h b ] 


( 220 ) 
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In the above two equations the subscript a stands for absorbing gas and 
the subscript b stands for broadening gas (such as N 2 ) . The quantities Pu 
and i h are the pressure of the absorbing gas and the path length of zone h, 
respectively. Ferriso et al. [52] give details in the determination of the band 
model parameters S°/d and y°/d for H 2 0, C0 2 , and CO. 

The transmittance in equation (214) may now be obtained as follows: 


-In t.(i/ f A»/ ) = 27r(y/d).f(x.) , 


where (y/d) . and x. are evaluated from equations (219) and (220). 

Two other expressions for W/d may be used instead of the one in 
equation (216) , which involves Bessel functions. One of these expressions is 
the exponential probability distribution described by Ferriso et al. [52]* 

W S°l ( s°l \ ~2‘ 

d d y + 7 iy°P J ’ (222) 


which has a more gradual transition region between the two asymptotic regions 
(weak- and strong-line regions) as a result of intensity distribution among the 
spectral lines. 


The other expression for W/ d is the effective absorption coefficient 
expression of Ferriso, Ludwig, and A bey ta [145]: 


W 

T = V 


1 + 


*eff^ V 

4 (y/d) 


(223) 


where 


— — _ 2tt y 

K eff ~ d~Pi d * (224) 

a 

Plass [91] shows that equation (223) and the W/d expression in equation 
(216) agree to within 10 percent. 
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c. Generalized method with modified Curtis -Godson approximation. 
The generalized method uses an integrated form of the Curtis-Godson equations, 
(217) and (218) , to obtain radiation from inhomogeneous gases in which the 
Beer-Lambert equations, equations (140) and (163) , are not applicable. This 
method, which has been proven to be accurate for strong concentration and 
temperature gradients, also takes into account the Doppler broadening effects 
which are important at low pressures (~0. 1 atm) . Ferriso et al. [52] and 
Huffaker [4,51] describe this method in detail, and Alligood [146] and Conway, 
Yossa, and Alligood [147] describe a computer program involving this method. 

The equation for heat transfer rate per unit area in the generalized pro- 
gram, which is similar to equation (208) except for the transmittance factor, 
t„, is 


dt 

q/A - f f f f B (T) ds sin0 cos0 d0 d<p dr . 
p <p 6 s 


(225) 


The transmittance t„ is equal to n t (i) , a product over the various 

i v 

radiating species in the rocket exhaust. The quantity t^(i) may be calculated 
by two methods. The first method, described by Alligood [146] , based on the 
investigation of Ferriso et al. [52], considers only the larger of collision and 
Doppler broadening (but not both as does the second method) . For the first 
case, the transmittance may be written in terms of the statistical band models 

as 


t ( i) = exp 
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(226) 


where the notation in braces implies the larger of D c (i) and D D (i) and the 
subscripts c and D stand for collision and Doppler broadening, respectively. 


The quantity D„(i) may be written similarly to equation (223) as 


D (i) = F(i) 
c 


, 1 [ F(i )l 2 

4 A c (1) 


(227) 
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and the quantity D D (i) may be written as 


1.7A D (i) 
D D (l) = F(i) 


in 


1 + 


0. 589[F(i)l 

A D (i) 



(228) 


The quantity F(i) is defined as 

s r s S(i) , 

F (i ) =/ P (i) v,„ ds = / -duT ds> 


(229) 


where S(i) is the line intensity originally defined in equation (142). 

The parameters A c (i) and A D (i) represent the modified Curtis-Godson 
approximation and are equal to 
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and 


A D (i » 


V 1 * 

J d(i) 


dF(i) 

ds 


ds , 


where 


(231) 


dF(i) 

ds 


p(i) K^(i) 


(232) 


The fine structure parameters (Doppler half-width, Y D , collision 
half- width Y c , and line density, 1/d) are given by Alligood [146] for H 2 0 (all 
bands) , CC^ (4. 3 ju and 2. 7 p bands) , and CO (5 p bands) . The tables of k„ 
as a function of wave number and temperature (300° K, 600° K, 1200 K, 

1500° K, 1800° K, 2400° K, and 3000° K) are found in Ferriso et al. [52]. A 
more accurate listing of absorption coefficients and fine structure parameters 
for H 2 0 has been tabulated by Ludwig et al. [148] for temperatures between 
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1200° K and 3000° K and for path lengths of 2, 5 and 10 feet. Herget, Muirhead 
and Golden [149] present revised tables of band model parameters for H 2 0 for 
temperatures between 300° K and 1400° K and for pressures from 0. 1 to 2. 0 
atm. 

The other method of evaluating the transmission t is given by Ferriso 
et al. [52] which combines both the Doppler and collision broadening in the 
expression 

i 

t„(D =F(i) (1 -y"V , (233) 

where 

2 " 2 2 -2 

y = { i - [D c (i)/F(i)J J. + { 1 - [D D (i)/F(i)J J - 1 . (234) 

It is felt that the generalized method (including the modified Curtis- 
Godson approximation) of predicting radiation heat transfer from rocket ex- 
hausts is the best of all of the methods described in this section. Huffaker 
[51] states that the difference between the measured inhomogeneous gas trans- 
mittance and that calculated by the modified Curtis-Godson approximation is 
within experimental error and not greater than 2 percent. 

Using the generalized method with the modified Curtis-Godson approxi- 
mation, calculations of spectral intensity have been made for the J-2 engine 
exhaust, as described by Huffaker [51]. Figure 13, taken from Huffaker, shows 
a comparison between calculations and measured values of spectral intensity 
from a J-2 engine exhaust during a sea level static firing. Figure 14 shows 
theoretical radiation calculations compared with experimental data obtained 
from the exhaust of a 1/45-scale F-l engine fired at Cornell Aeronautical Lab- 
oratory which used 0 2 and C 2 H 4 as propellants. The radiation calculations for 
this engine show better agreement with the data than for the J-2 engine (Fig. 13) ; 
however, it is known that the spectrometer used for the J-2 test was not cali- 
brated properly after the firing. Additional calculations will be made by Ther- 
mal Environment Branch to compare with spectral intensity measurements 
taken during altitude firings (in which there is no atmospheric absorption) of 
the J-2 engine. These measurements will be performed in the AEDC J-4 test 
cell at Tullahoma, Tennessee, in the spring of 1967. 
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J-2 Nozzle 


P 0 * 715 
O/F * 5.5 



FIGURE 13. COMPARISON OF MEASURED AND CALCULATED 
INTENSITY FOR J-2 ENGINE FIRED AT SEA LEVEL 

This concludes the section on gaseous radiation. It is felt that the basic 
points of the theory such as radiation from accelerated charges, shape and 
broadening of spectral lines, band models, and methods of predicting gaseous 
radiation from rocket exhausts have been covered sufficiently. For additional 
detail on gaseous radiation the more than 90 references mentioned in this 
section should be consulted. For details on the rocket exhaust plume programs 
which must be generated before the radiation heating can be predicted, it is 
suggested that references such as Farmer, Prozan, Ratliff, and McGimsey [1], 
Prozan [150], and Ratliff [151] be consulted. 

RADIATION FROM CARBON PARTIClfS 


Thermal radiation from carbon particles is the dominant means of radi- 
ant heating from exhausts of hydrocarbon-fueled rocket engines such as the H-l 
and F-i, which power the first stages of the Saturn I and Saturn V vehicles. 
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FIGURE 14. COMPARISON OF MEASURED AND CALCULATED 
INTENSITY NEAR THE EXIT OF A 1/45-SCALE F-l ENGINE 
FIRED AT CAL AT 120, 000 FEET 

respectively. These engines burn RP-1 (kerosene) as the fuel as opposed to 
other hydrocarbon-fueled engines which burn ethyl alcohol (ethanol) or aerozine- 
50, and consequently produce little or no carbon in their exhausts. 

The carbon emission spectrum from H-l and F-l engines is practically 
continuous (follows a graybody or blackbody distribution, depending on optical 
depth) near the exit plane in the 1-ju to 4-p wavelength region. This is in con- 
trast to the emission spectrums of the H 2 0, CC >2 and CO gaseous molecules 
which are frequency dependent (emit in certain bands only) in the infrared 
wavelength region (as shown in the previous section) . At distances far down- 
stream of the exit (at altitudes less than about 100,000 feet) when afterburning 
is still present, the carbon spectrum is no longer continuous. Also, if a spec- 
trometer which views the plume is mounted a considerable distance from the 
exit, atmospheric absorption bands can absorb the carbon continuum as seen by 
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the spectrometer in the 1. 4— , 1 . 9— pt , 2.7-pand 6.3-p H 2 0 bands and 2.7 -p and 
4.3-/U C0 2 bands. This carbon continuum spectrum can be seen in such reports 
as Simmons [152] , de Bell, Simmons, and Levin [ 153] , Levin, Wagner, and 
Thomson [154], an Aerojet-General report [155], Wagner [156], and in Figure 
15 of this report. Figure 16, presented for comparison with Figure 15, shows 
the spectrum of an exhaust from an N 2 0 4 /UDMH propellant engine. The lack of 
carbon continuum for the exhaust of this non-kerosene liquid-fueled engine is 
easily seen. In Levin, Wagner, and Thompson [154] evidences of chemilumi- 
nescence can be seen by noting the peaks of C 2 , CN, and CH emission above the 
carbon continuum. 



4 .8 1.2 1.6 2.0 2.4 2.8 3.2 

Wavelength, X, microns 


FIGURE 15. SPECTRUM OF LOX/RP-1 PLUME SHOWING 
CARBON CONTINUUM 

Accurate predictions of radiant heating from carbon particles is difficult 
for a number of reasons. In the gas generator, carbon is a nonequilibrium 
product of combustion because of the rather low temperatures (~900° K) involved; 
hence, its concentration in this region must be calculated with finite reaction- 
rate chemical kinetics. It is possible, however, that carbon is in equilibrium 
in the combustion chamber, nozzle, and near-field of the exhaust because of the 
higher temperatures in these regions. The concentration of carbon in these 
regions should be calculated by present equilibrium thermochemical programs; 
however, this has not yet been accomplished satisfactorily because of the un- 
certainties in the O/F distribution in these regions. A plot of carbon concen- 
tration versus O/F ratio for various temperatures shows that, for a small 
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FIGURE 16. SPECTRUM OF N 2 0 4 /UDMH PLUME SHOWING LACK 

OF CARBON CONTINUUM 

increment in O/F ratio (for a given temperature) , a very large increment in 
carbon concentration will result. To treat radiation from carbon particles 
adequately, therefore, the mechanism of carbon formation must be fully under- 
stood, the size distribution of carbon particles must be known, the optical 
properties of carbon ( dispersion theory, scattering phenomena, etc.) must be 
analyzed, and the emissivity of carbon particles must be calculated (or meas- 
ured) . 

Scattering of radiation by carbon particles is frequently assumed negligible, 
although this phenomenon may not be legitimately neglected if the size of the 
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particles is large and the wavelength under consideration is small. Scattering 
by carbon particles may thus be calculated by Rayleigh scattering theory in which 
“ « 1* where a = 2 n r^/X (r p is the radius of the particle) , and by Mie 
scattering theory in which a ~ 1 . Since carbon particles usually fall in the 
realm of Rayleigh scattering, this theory is discussed in this section; the Mie 
theory is discussed in the following section since it applies mainly to larger 
(AI 2 O 3 ) particles. This section discusses in succession the formation of car- 
bon particles in luminous flames, the determination of carbon particle sizes, 
the dispersion equation for carbon, the Rayleigh scattering theory, and methods 
of predicting absorption coefficients and emissivities for carbon particles. 


Formation of Carbon Particles in Luminous Flames 


In recent years there has been a tremendous amount of literature written 
regarding the mechanism of carbon formation in luminous flames (including 
rocket and turbojet combustion chambers and exhausts) . Some of the more 
detailed reviews pertaining to carbon formation include the books by Gaydon 
[157], Gaydon and Wolfhard [158] and Palmer and Cullis [159] , and the articles 
by Parker and Wolfhard [160], Porter [161] and [162], Thomas [163], Street 
and Thomas [164], Daniels [165], Singer and Grumer [166], Tesner [167], 
and van der Held [168] . Other less detailed articles regarding carbon formation 
which also contain much experimental data were written by Gaydon and Fairbairn 
[169], Arthur and Napier [170], Scully and Davies [171], Fenimore, Jones, 
and Moore [172], Ferguson [173], Lee, Thring, and Beer [174], and Stehling, 
Frazee, and Anderson [175]. Some of the many books and articles which dis- 
cuss combustion in luminous flames and combustion chambers ( some of which 
also contain descriptions of carbon formation) were written by Godsave [176] , 
Khitrin [177], Wicke [178], Yagi and Kunii [179], Behrens [180], Paushkin 
[181], Minkoff and Tipper [182], Fristrom [183], Penner [184], McCafferty 
and Hibbard [185], Penner and Datner [186], Wise and Agoston [187], Spalding 
[188], Williams [189], Fristrom and Westenberg [190], and Bahn [191]. 

In a rocket combustion chamber such as that of the H-i or F-l engine, the 
fuel (RP-1) and the oxidizer (liquid oxygen) are pumped into the chamber by 
a turbopump, which is driven by a gas turbine. The gas turbine is, in turn, 
driven by the exhaust of a gas generator. Inside the combustion chamber, 
combustion takes place at a relatively high O/F ratio (~2. 5) and at a high 
temperature (~3000° K) . However, in the gas generator, combustion takes 
place at a low O/F ratio (~0. 4) and at a low temperature (~900° K) . Before 
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the RP-1 fuel is injected into the combustion chamber, however, it is routed 
around the chamber and nozzle to cool these components regeneratively. Carbon 
may thus be formed in the relatively high O/F ratio, high temperature region of 
the combustion chamber or in the relatively low O/F ratio and temperature 
regions of the gas generator and the walls of the combustion chamber and nozzle. 
It is expected that more carbon will be formed in the latter regions because of 
the excess fuel and lower temperature experienced in these regions. 

Paushkin [181] states that the formation of carbon in combustion cham- 
bers depends upon a number of factors such as degree of atomization and tem- 
perature of the fuel and oxidizer, design of the combustion chamber (including 
its chamber pressure) , and surface tension, viscosity, and composition of fuel 
(aromatic hydrocarbon content) . Since RP-1 fuel (kerosene) contains a large 
amount of aromatic (cyclic) components (~20 percent) compared to aliphatic 
(straight-chained) components, it is usually assumed that the majority of the 
carbon formed comes from this aromatic portion of the fuel. 

In a given combustion chamber, Penner and Datner [186] state that the 
amount of carbon deposited, for a wide variety of hydrocarbon fuels, can be 
correlated by an empirical equation of the form 

T 

C = ~ In (K* R - K 2 ) + ~~ - K 5 , (235) 

iS.3 lv 4 


where C is the rate of carbon deposition, R is the carbon-to-hydrogen ratio, 
T b is the boiling temperature of the fuel, and Kj to K 5 are empirically deter- 
mined constants such as found in the article by Starkman, Catteneo, and 
McAllister [192] . Hie rate of combustion of a cloud of carbon particles in 
luminous flames, according to Lee, Thring, and Beer [174] may be determined 
from the following equation based upon experimental data: 
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where P c is the density of carbon particles (soot) , d Q is the initial diameter 
of the particles, m Q is the initial mass flow, m is the mass flow after time, 

t, P_ is the partial pressure of oxygen, T is the particle (or flame) temper- 
°2 

ature, and R is the gas constant. According to Howard and Essenhigh [193], 
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the complete sequence of events for carbon particle combustion takes place in 
four steps: heating, ignition, pyrolysis, and heterogeneous combustion. Gas- 
phase combustion will occur simultaneously with heterogeneous combustion on 
the surface of the particle if the surface flux of volatile pyrolysis is small; 
however, if the surface flux is large, the reaction zone will be forced away 
from the solid surface, thereby screening the carbon particle from oxygen 
attack. 


Flames occurring in combustion chambers of liquid-fueled rocket en- 
gines are neither diffusive nor completely premixed; hence, the general theories 
of carbon formation in both diffusion and premixed flames must be considered. 
According to Gaydon and Wolfhard [158] , in a pure diffusion flame, the com- 
bustion processes are assumed to take a fairly long time for their completion; 
however, there is little or no contact between the oxidizer and fuel. For purely 
premixed flames, the oxidizer and fuel are in direct contact, but the combustion 
processes are limited by the very short time of passage through the reaction 
zone. Carbon formation theories differ, depending on what type of flame 
(diffusion or pre mixed) occurs, but also these theories differ from author to 
author, so that there is no one general theory available which will adequately 
describe carbon formation. In the remainder of this section some of these 
various carbon particle (and ion) formation theories are briefly analyzed. 

1 • The C», Atomic Carbon, and C 3 Condensation Theories. The C 2 
condensation theory was originally developed by Smith [194] , and has been 
discussed by Gaydon and Wolfhard [158] and Parker and Wolfhard [160], 

Smith, after performing a spectroscopic study of ethylene burning in air, sur- 
mised that carbon was produced by a polymerization of C 2 molecules. He also 
discussed the experiments of Klemenc, Wechsberg, and Wagner [195] , who 
investigated the decomposition of carbon suboxide, C 3 0 2 , and found C 2 bands in 
absorption. After passing through a reddish stage in color, the C 2 rapidly 
polymerized to normal solid carbon. Gaydon and Wolfhard [158] and Parker 
and Wolfhard [160], however, believed that this C 2 condensation theory is not 
valid in practice because of the low times of reaction and concentration usually 
involved. Gaydon [156] states that semi-quantitative estimates of the amount 
of C 2 in flames, studies of flash photolysis, and the relative positions of C 2 and 
carbon in flames also tend to make this theory unbelievable. 

Gaydon and Wolfhard [158] state that it is possible that free carbon 
atoms may be condensed directly from flames since, in equilibrium with graphite, 
the concentration of C atoms is several orders of magnitude greater than that of 
C 2 . Cabannes [196] believes that solid carbon in flames could be formed by the 
condensation of carbon vapor whose main constituent is assumed to be C 3 . 
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2. The Hydrocarbon Polymerization-Condensation Theory . This theory 
involves the general pyrolysis (or decomposition) of pure hydrocarbons as dis- 
cussed by Gaydon and Wolfhard [158] , Porter [161, 162] , Palmer and Cullis 
[159], Thomas [163], Street and Thomas [164], Gordon [197], and many 
others. According to Gaydon and Wolfhard [158] , during the pyrolysis of 
hydrocarbons, a number of processes are found to occur, such as hydrogenation, 
dehydrogenation, cracking to simple hydrocarbons, polymerization, and some- 
times condensation to aromatics. The particular type of pyrolysis depends 
upon the temperature and the presence of catalysis. The basis of all the hydro- 
carbon polymerization theories is that polymerization of the fuel occurs first, 
aided sometimes by oxidation, followed by dehydrogeneration (in which oxidation 
may be a factor) , finally resulting in carbon being the end product. Street and 
Thomas [164] present an elaborate diagram which shows probable interrelation- 
ships among the possible routes to the formation of carbon involving polymeri- 
zation. 


Thomas [163] gives a description of "soot" formation from hydrocarbon 
fuels using the polymerization-condensation theory, and states that "soot" is 
not carbon, but has a polybenzenoid hydrocarbon stiucture. Rummel and 
Veh [198] were the first to suggest that aromatics and polycyclic hydrocarbons 
preceded the formation of carbon (or soot) . A current investigation at Rocket- 
dyne [199], performed under Thermal Environment Branch contract, is under 
way to attempt to determine the chemical formulation of the soot molecule. 

3. The Acetylene Theory. Porter [161,162] and Anderson [200] state 
that, while low temperature (< 1000° C) pyrolysis of hydrocarbons tends to 
produce higher polymers, at high temperatures (> 1000° C) decomposition of 
hydrocarbons predominates over polymerization. Because of these high tem- 
peratures, Porter [161] concludes that, when the time for half reactions is of 
the order of one second or less, the pyrolysis of hydrocarbons does not result 
in polymerization, but in decomposition to smaller molecules. Portei [161] 
further concludes that the thermal decomposition of hydrocarbons results in 
dehydrogenation and cracking to smaller molecules, and the last stable hydro- 
carbon to be observed before carbon formation is acetylene. It is known that 
acetylene, which is a very endothermic substance, is then very easily decom- 
posed to carbon and hydrogen. Discussions of experimental investigations of 
carbon formation from acetylene have been reported by Stehling, Frazee, and 
Anderson [175] and by Westbrook, Hellwig, and Anderson [201]. 

4. The Surface Decomposition Theory . This theory, which was dis- 
cussed by Tesner [167] , involves the formation of carbon through decomposition 
of hydrocarbons on the particle surface. Tesner studied the formation of 


86 


carbon black in an externally heated reaction tube and concluded that the process 
took place in two stages: (1) nuclei formation and (2) particle growth. He 
found that when a hydrocarbon is thermally decomposed in a flow system, the 
spontaneous formation of nuclei of carbon black particles in the main part of the 
gas is accompanied by the simultaneous surface formation of a solid carbon 
layer on the reaction chamber walls. 

5. The Boudouard Reaction Theory. It has been proposed by Behrens 
[180] and van der Held [168] that the Boudouard reaction, which is the reaction 
predicted in an equilibrium chemistry analysis, 


2 CO ^ C0 2 


+ '"SOLID 


(237) 


may contribute to carbon formation in premixed flames. This reaction, which, 
in the presence of hydrogen, proceeds through the stages 


CO + OH ^ C0 2 + H and CO + H ^ c SOL1D + OH > (238) 


is usually so slow, however, that in typical reaction times in combustion cham- 
bers (~ 5 x 10~ 3 sec) , the Boudouard reaction would be far from complete. 
Foster [202] states that carbon (soot) is not likely to appear in equilibrium 
in the products of fuel-rich premixed combustion because the decomposition of 
hydrocarbons to carbon proceeds more rapidly than oxidation reactions of hydro- 
carbon or soot with C0 2 or H 2 0. These latter reactions would cause the system 
to be in equilibrium if it were allowed to proceed for a sufficient length of time. 


Carpenter, Foreman, and Goldstein [14] state that to determine the 
amount of free carbon with considerable accuracy, it would be necessary to 
consider the dissociation of each constituent in the exhaust gas and to match 
the dissociation of each constituent subject to the constraint that the total num- 
ber of atoms remain constant in the resulting mixture. Such calculations are 
given by Gaydon and Wolfhard [158], Huff, Gordon, and Morrell [203], Penner 
[204] , and Dodge [205] who states that carbon deposition for the Boudouard 
reaction is possible only if ec t ua l 1° or greater than Kp/P, where 

Kp is the equilibrium constant for the process, P is the total pressure of the 

mixture, and n and n are the number of moles of CO and CC^, respec- 
, CO C Oo 

tively. 
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6. "Vapor Pocket" Theory. Boynton [206] relates a possible theory for 
carbon formation in rocket combustion chambers using the "vapor pocket" 
theory. He believes that the individual drops of fuel and oxidizer in the chamber 
evaporate soon after injection to form small vapor "pockets” whose size depends 
upon the characteristics of the particular injector. These vapor pockets would 
thus interact with each other, burn at the edges, and then break up, coalesce, 
and circulate very rapidly. The reaction zone would transfer heat by conduction 
into the fuel vapor pocket interior, breaking the vapor into fragments which 
react to form carbon particles. The carbon would be more apt to form from 
the larger vapor pockets in which the fuel vapor volume -burning surface ratio 

is larger, and would also be more apt to form (Boynton and Neu [201]) when 
less efficient injectors are used. This theory is compatible with the combustion 
inefficiency calculations of Farmer et al. [ 1] . 

7. New Carbon Formation Theories. Several new theories on carbon 
particle formation were discussed at the Eleventh Symposium (International) on 
Combustion held at Berkley, California, August 14-20, 1966. Homann and 
Wagner [208] discussed experiments of acetylene and benzene premixed flames 
and stated that for acetylene flames, the polycyclic aromatic hydrocarbons 
form some distance behind the oxidation zone and are similar to those evaporat- 
ing from soot samples when heated in a mass spectrometer under reduced 
pressure. Their concentration profiles in acetylene flames showed that they 
were by products instead of intermediates for carbon. For benzene flames, 
however, it was found that these polycyclic aromates appeared during fuel 
oxidation such that their concentration was much larger in the region where 
solid carbon formed. Echigo, Nishiwaki, and Hirata [209] disucssed spectro- 
scopic analyses of ten hydrocarbon flames and proposed a new theory that the 
combustion process results in dehydrogenation and polymerization reactions in 
the fluid phase. They found the resulting, unsaturated hydrocarbons ("pre-soot 
substances") emitted band spectra rather than a continuum and that eventually 
soot particles were agglomerated after dehydrogenation and polymeration were 
complete. 

Also in the Eleventh Symposium, Tesner, Snegyreva, and Soorovikin 
[210] discussed formation of carbon by thermal decomposition and incomplete 
hydrocarbon combustion, describing a two-step process: the generation of 
carbon nuclei and the growth of the carbon particles. They measured rates 
of carbon particle formation during thermal decomposition of acetylene in 
laminar diffusion flames and decomposition of aromatic hydrocarbons when 
mixed with a turbulent flow of hot combustion products. They found that, be- 
cause of chain branching, soot particle formation occurred with a tremendous 
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speed, but because of quadratic terminations, was rapidly suppressed; hence 
the curve of rate of particle formation acquired a sharp peak. Cullis, Read, 
and Trimm [211] discussed another theory at this symposium and stated that it 
is unlikely that carbon formation always involves the intermediate production 
of acetylene. They conducted experiments on the pyrolysis of acetylene and 
vinylacetylene (a polyunsaturated compound) at 500° C and 700° C and found, 
in contrast to a wide variety of acetylenic products formed during acetylene 
pyrolysis, vinylacetylene was converted directly to solid carbon products. 

8 carbon Ionization Theories. Certain investigators such as Singer 
and Grumer [166], Carpenter, Foreman, and Goldstein [14], Einbinder [212], 
Shuler and Weber [213] , Smith [214] , and Sugden and Thrush [215] have shown 
that carbon particles can be responsible for ions observed in flames and rocket 
exhausts. Singer and Grumer [166] believe that the ionization in hydrocarbon- 
fueled flames must be caused by accumulating species which are either carbon 
nuclei or carbon precursors of low ionization potential. This ionization poten- 
tial (or work function) for carbon particles has a low value (between 4. 35 eV 
for solid carbon and 13. 3 eV for gaseous carbon) compared to the reaction 
intermediates and products of flames (H 2 0, CO, C0 2 , OH, O, and H) which 
have ionization potentials of between 12 and 16 eV. 


Shuler and Weber [213] studied the ionization of carbon particles in 
acetylene -oxygen flames burning in air and proposed that the free electron 
concentration be given by the Saha [216] or equilibrium thermionic emission 
equation 


N = K = 
e 


2 (27r m k T) 
e 


3/2 


.-** T 


(239) 


where <p is the carbon particle work function, m g is the mass of the electron, 
T is the temperature in degrees K, and k and h are the Boltzmann and Planck 
constants, respectively. Sugden and Thrush [215] also studied carbon particle 
ionization in acetylene flames and proposed the following modified equation: 

2 (27T m k T ) ^ ^ 

N = 7T exp[- (c/> + N e 2 /r N )/kT] , (240) 

eh e P p 


where r and N„ are the radius and concentration, respectively, of the carbon 
particles and e is the electronic charge. Einbinder [212] presents a further 
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modification to equation (239) by writing the ratio of the 
to the carbon particle concentration as 


electron concentration 


N e /N D = < r n k T /e 2 ) In (K/N ) + ~ , for 
c P p e 2 


(N /N > 3) 
e p 


(241) 


where K is determined from equation (239). 

Pergament [217] and Pergament and Calcote [218] , Kurzias [219] and 
other workers at AeroChem have investigated chemi-ionization effects in rocket 
combustion chambers and exhausts when solid carbon is present Thev state 
that the electrons and chemi-ions such as CHO + and H 3 0 may be formed if 
carbon gas is m equilibrium with hydrogen atoms in the process 


C GAS + H 


CH 


(242) 


the CH combines with free oxygen to produce CHO + 
process 


ions and electrons in the 


CH + O — CHO + + e 

The CHO + ions are then combined with water 
in the process 


to produce hydronium 


ions, 


(243) 

h 3 o + , 


CHO + H 2 0 — H 3 0 + + CO 


(244) 


The Thermal Environment Branch at MSFC is sponsoring an experimen- 
tal program at altitude at Cornell Aeronautical Laboratory involving the use of 
resonant cavity probes to determine electron concentrations in an ethylene 

Zs e a n rner e ^ "*** ^ engine) * Later > 11 is P^ed to make 

type measurements in the plumes of model solid propellant ullage and 

retro motors of the Saturn vehicle. The theoretical programs at AeroOhfm 

and at Boeing [220,221] will be used to check the experimental data. 
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In summary, the problem of determining a suitable theory for formation 
of carbon particles and ions in flames and rocket exhausts is very difficult. The 
Thermal Environment Branch, with assistance from one of its contractors, 
plans to initiate a program soon to determine carbon particle formation and 
concentration associated with H-i and F-l engine exhausts. This program 
would involve basically the study of the mixing processes of RP-1 and liquid 
oxygen, the study of the chemical composition of RP-1 and its pyrolysis products, 
and the chemical reactions leading to soot formation. Another program spon- 
sored by this branch which is being performed at General Applied Science 
Laboratories by Edelman [222] involves the analytical investigation of plume 
afterburning associated with H-l and F-l engines. In this analysis the diffusion 
equations are combined with the finite-rate chemical kinetic equations to handle 
the mixing and combustion along the plume boundary. The radiation to the 
vehicle base from carbon particles present in the mixing layer is then analyzed 
in the program. A preliminary report has been completed on this subject [223] . 


Carbon Particle Size Determination 


Accurate knowledge of carbon particle sizes in flames and rocket exhausts 
is an important factor in determining the scattering and emission processes and 
consequently the thermal radiation from these flames and exhausts. Measure- 
ments of the size of carbon particles in luminous flames and rocket exhausts 
vary, depending upon the various methods used to obtain the particular samples. 
In general, however, most investigators have found that the radius, r , of 
carbon particles lies between 50 A to 1000 A- Wolfhard and Parker [224] have 
obtained electron micrographs for carbon particles in hydrocarbon flames 
showing this 50 A to 1000 A range for particle radius. Parker and Wolfhard 
[160] found that in non-smoky acetylene flames all particles were of about the 
same size (r p ~ 50 A) ; however, in a smoky flame large particles of approxi- 
mately 250 A were also found. Lee, Thring, and Beer [174] , using a gaseous 
hydrocarbon mixture of 70 percent C3H 8 , 10 percent C 2 H 4 , and 20 percent 03!!^ 
found an initial average particle radius of 200 A, but found that the size of the 
particle decreased as the particle underwent combustion. Scully and Davis 
[171] experimenting with various aromatic hydrocarbons, found that the average 

o o 

particle radii ranged from 100 A to 400 A- 

Erickson, Williams, and Hottel [225] performed light- scattering meas- 
urements on soot in a benzene-air flame and also obtained electron micrographs 

o 

of samples on a probe showing an average soot particle radius of 125 A. Meas- 
urements of particle radius obtained by light-scattering measurements were 
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considerably higher (r p ~ 700 A) than the electron micrographs; however, it 
was determined that the soot in the flame was partially agglomerated such that 
the experimental data were strongly affected. Kunugi and Jinno [226] also 
used light- scattering techniques to obtain sizes and concentrations of soot 
particles in diffusion flames. These investigators calculated particle size by 
comparing the measured ratio of the two light- scattering intensities for perpen- 
dicular and parallel polarizations, with the intensities calculated by Mie scatter- 
ing theory, and obtained a value of r p of between 600 A and 1000 A- These 
values were 3 to 5 times higher than those obtained for soot particles collected 
on a solid surface inserted into a flame; however, these authors state that the 
light- scattering method used tended to measure only the larger particles which 
scatter the light more effectively than the smaller ones. Singer and Grumer 
[166] obtained particle samples on stainless steel screens in rich, flat flames 
of propane-air and ethylene-air. These investigators also obtained evidence of 
agglomeration in the shape of long filaments with opaque granules at the tips or 
bends of the filaments. Particle radii as large as 2500 A were observed on the 
filaments. Millikan [227,228] attempted to circumvent the agglomeration 
problem by sucking the particle-containing gas through a small orifice into a 
quartz tube, expanding it to a low pressure and catching the particles upon a 
cold target. The carbon particles obtained in this manner were about 50 A in 

o 

radii 6 mm from the burner exit and grew to about 200 A in radii at 12 mm 
above the burner exit. 

Tesner [167] also obtained carbon particle size distributions in various 
hydrocarbon flames and found virtually no particles at all below a radius of 
50 A, and an average particle radius of about 200 A. He also found carbon 
concentrations ranging from 2 x 10 10 to 2 x 10 12 particles per cm 3 . Johnson 
and Anderson [229] recently investigated carbon particles formed by pyrolysis 
at 500° C to 1000° C for a large number of hydrocarbon fuels. They found that 
the various types of solid carbon particles formed varied from small, isolated 
specks of material, about 50 A in radius up to large dense particles about 
2500 A in radius, which were approximately spherical in shape and which were 
hooked together into networks of chain-like structures. It was found that, by 
varying the pyrolysis temperature, contact time, and hydrocarbon concentration, 
there appeared to be a gradual transition between one type of particle to another 
with no sharp dividing line between particular types at any stage. 

Boynton [206] obtained carbon particle samples on flat plates from the 
exhaust of small rocket engines (P c = 500 psia, nozzle exit diameter = 1.5 
inches, A e /A t = 8) using RP-1 as the fuel and liquid oxygen as the oxidizer. 

The flat plates were mounted on a water-cooled pendulum which was moved 
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rapidly through the plume. The pendulum was swung through the plume at various 
distances downstream of the exit, with the amount of deposited carbon determined 
by weighing the disks on an analytical balance. It was found that the amount of 
carbon on the flat plates decreased with increasing distance downstream and that 
the average radius of the particle was approximately 175 A, the size distribution 
being Gaussian. 

In a recent investigation performed at Rocketdyne [199] under a Thermal 
Environment Branch contract, carbon particle sizes were determined from full- 
scale F-l engine exhausts. The particles were pumped through a water-cooled 
probe shown in Figure 17 which was immersed 1. 12 inches and 4. 12 inches into 
the plume for two engine firings. It was shown that the average particle radius 
ranged from 100 A to 300 A, a range which agreed well with previous samplings 
taken from model F-l engine exhausts. 



FIGURE 17. ROCKETDYNE F-l HOT GAS SAMPLER SHOWING 
4. 12 -INCH PROBE IMMERSION 

Ferriso et al. [230] and Boynton et al. [231] also report on measure- 
ments of carbon particle size in rocket exhausts. They, as did Boynton [206], 
collected the particles on small water-cooled plates which were rapidly passed 
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through the plume at four inches downstream of the nozzle exit. The samples 
of carbon particles were taken at various oxidizer-fuel ratios for three engines 
of area ratio, A e /At, of 5.25, 3.0 and 1.5 (See this section, subtitle "Methods 
of Predicting Carbon Absorption Coefficients and Emissivities, " part 4. ). How 
ever, it was found that at oxidizer-fuel ratios smaller than 1. 6 the carbon 
deposits were too heavy and the particles agglomerated considerably during the 
collecting process such that it was impossible to distinguish individual particles. 
Using stereographic photomicrographs with shadow casting, it was found that at 
O/F ratios greater than 1.6 the mean value for r p was approximately 200 A, 
a value which is representative of all the mean particle radii mentioned in this 
section. 


Dispersion Equation for Carbon 


In this section the dispersion equation for carbon is discussed and the 
relations for conductivity, a, polarization, P, dielectric constant, e, and index 
of refraction n = ^ - in 2 are determined. The Rayleigh theory of scattering by 
carbon particles is dependent upon this dispersion theory, as will be seen in 
the next section. 

The classical theories of dispersion in gases, liquids and solids, and 
metals (conducting materials) are described by Stratton [232], Slater and Frank 
[233,234], Born and Wolfe [235], Panofsky and Phillips [57], Sommerfeld [236], 
Bottcher [237] , Fleagle and Businger [238], and Wood [239] . A dispersive 
medium is defined as one in which the index of refraction of the medium varies 
with frequency of propagation. Another definition sometimes used is that a dis- 
persive medium is one in which the phase velocity v of a superposition of 
harmonic wave trains of infinite length and duration varies with the frequency 
of propagation. In a nonconducting (dielectric), nonmagnetic medium, a wave 
is propagated with a phase velocity v = c/n, where c is the velocity of light in 
empty space and n is the index of refraction equal to (e/n) 2 « e^, since the 
permeability p is usually assumed close to unity (except for ferro -magnetic 
materials) . 

For gases, the dispersion theory was actually discussed in the previous 
section, subtitle "Radiating from an Accelerated Charge," to predict the radi- 
ation (absorption) by a classic dipole. In this manner, an electron of charge 
e was acted on by a force eE of an external electric field E = E o e- ic0t and 
received a linear restoring (elastic) force, -Gx = -iiicOq x, proportional to the 
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displacement x of the electron, and also a damping or dissipative force, -myx, 
proportional to the velocity of the electon. This equation, (86) , will now be 
rewritten in the form 


mx + me o x + myx = eE 
o 


(245) 


which has the solution 



(246) 


It may be noted that in the remainder of this section (and the following section) 
the quantities, E, J, P, D, x, x, and x are all vector quantities. 


The polarization, P, will now be defined as 


P = e 


E 


N k\ 


N. e 2 /m 

E y -j— 

Y or - or - icoy 
k o r k 


(247) 


where is the number of electrons per unit volume characterized by the con- 
stants c and y^. 

The dielectric constant e will now be written as 


e = 1 + N -fj— , (248) 

o 

where (3 is the polarizability and is equal to the dipole moment ex divided by 
the electric field E. 


Hence, equation (248) may be written as 


N e 2 /me 

1 + T _ ,.,2 TUZ = < n l - in 2> 2 • 


, u ) 2 - u ) 2 - icoy. 

k o, 'k 

k 


(249) 


95 



For an absorbing gas the dielectric constant is a complex quantity and 
may be written in terms of its real and imaginary components, n t and n 2 . 
Furthermore, the second term in equation (249) is usually small in gases with 
respect to unity so that the square root of equation (249) can be written as 


1 

€ 2 = nj - in 2 = 1 + 



N. e 2 /me 
k o 

O 9 • 9 

or - CO - icoy 
\ k 


(250) 


in which the real part rq and imaginary part n 2 are, respectively, 


n i 


1 (N k eVme o ) ( "°k - 


(c o 2 . - CO 2 ) 2 + co z y 


,2 .,2 


(251) 


and 


n 2 




(N k e 2 /me o ) uy k 


(co 2 - w 2 ) 2 + co 2 y, 
o, k 

k 


(252) 


As the frequency increases, the index of refraction n experiences the 
phenomenon of "anomalous dispersion" in the vicinity of each of the resonant 
frequencies, co 0k . This phenomenon can be seen in Figure 18 as at o> 0 the 

quantity n t - 1 changes sign while the quantity n 2 (and hence the absorption) 
is a maximum. 
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For solids and liquids it cannot be assumed that the force acting on an 
electron is merely equal to the electric vector of the electromagnetic wave in 
free space. Instead, the added force on the electron caused by the polarization 
of the body must be considered. According to Slater and Frank [234] , a small 
sphere with its center at the position of the electron in question may be consid- 
ered to have been cut out of the medium. Charges are then induced on the sur- 
face of the sphere producing a force at the center of the sphere which now has 
a total electric field at the center of the sphere of 


E' 



(253) 


This contribution to the field because of polarization of the body is also 
known as the Lorentz field. The polarization, P, as given in equation (247) 
now becomes 



(254) 


and an expression for the dielectric constant may be written as 


£ - 1 (llj + in 2) 2 - 1 = J_ y 

e + 2 (n t + m 2 ) 2 + 2 3 y- 


N, e 2 /me 
k o 


u) 2 - co 2 - icuy. 
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k 


(255) 


For conducting materials such as metallic or carbon particles, a slightly 
different approach to the dispersion theory is developed. It is assumed that not 
all of the electrons are bound to the atoms (i. e. , bound electrons, N^) , but 
some are free to move about within the material (i. e. , conduction electrons) . 
Both the bound and conduction electrons of the material are forced to oscillate 
by the incident, periodic electromagnetic wave. The total number of electrons 
per unit volume, N, may be thus written as 


N 


k 


-- Zd 
. + N - — 

k cm 


(256) 
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where Z is the atomic number, m a , the mass in grams of one atom of the con 
ducting material, and d is the density of the conducting material. 

It is usually assumed that the fixed, positive ions of the conductor com- 
prise a region of constant electrostatic potential. Within this region a cloud 
of free or conduction electrons wanders about and thus carries the current in 
the direction of an applied electric field, E. The motion of free electrons is 
opposed by collisions at the lattice points which are occupied by the heavy ions. 
The momentum transfer from the drifting electrons to the lattice points results 
in thermal vibrations of the ions and a damping of the electron motion. The 
resulting equation of motion is thus similar to equation (245) except that there 
is no restoring force; hence, u> o is zero. This equation takes the form 

mx + my x = eE , (257) 

c 



However, since the static conductivity a' may be defined as J - cr'E, a complex 
conductivity a may thus be defined as 

N e 2 /m i(e 2 /m) N 

(260) 

y - ico co + iy 

' c c 

The dielectric constant may now be written as 



e 

o 



(261) 
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where E is now E - E f - P/3e . 

o 


Hence, the dielectric constant may be written as 


e fa 
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(262) 
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But the Lorentz field may be assumed negligible in the visible and infrared 
compared with E; then E 1 = E, and equation (262) becomes 


e = € + ( e 2 / m) Yj 
k 




- 


(263) 


The complex propagation constants of a particular conducting sphere 
(carbon particle) and the surrounding medium (kj and k 2 , respectively) can be 
written in the form 


kj k 2 = k 2 = pt co 2 + iupco , 


(264) 


where p is the magnetic permeability of the medium and n equals ki/k 2 , 
according to Stull and Plass [240]. If equations (260) , (263), and (264) are 
solved simultaneously, the following equations for the real and imaginary parts 
of the index of refraction result: 
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Stull and Plass [240] obtained values for the constants in equations (265) 
and (266) from Halpern and Hall [241], who investigated fast-charged particles 
in carbon, and from Senftleben and Benedict [242], who measured nj and n 2 
for amorphous carbon at 2250° K in the visible part of the spectrum and the 
reflection coefficient (or reflectivity) in the infrared. This reflection co- 
efficient is defined by Fresnel's equation as 


( n i - l) 2 + 4 

P X (n 4 + l ) 2 + n 2 


(267) 


Stull and Plass [240] then solved equations (256) , (265) , (266) , and 
(267) for the four unknowns y c , 74 , N c and N 4 and obtained for a temperature 
of 2250° K the dispersion equations for amorphous carbon to be 

2 2 6.448 x 10 32 3.224 x 10 32 3. 224 x 10 32 

111 n 2 “ 1 + 4 . 062 x 10 35 - co 2 + 9. 549 x 10 33 - co 2 + 5. 217 x 10 33 - co 2 

6.348 x 10 32 (1. 966 x 10 32 - co 2 ) 3.05 x 10 31 

+ (1.956 x 10 32 - co 2 ) 2 + 1.369 x 10 33 co 2 2. 323 x 10 31 + co 2 (268) 


and 


_ 6.347 x 10 32 x 3.70 x 10 16 cu 

ni ” 2 “ (1. 956 x 10 32 - co 2 ) + 1. 3 6 9 x 10 33 co 2 

(269) 

3 . 05 x 10 31 x 4 . 82 x 10 15 
+ (2.323 x 10 31 + co 2 ) co 


These two equations have been used by several investigators to obtain 
emissivities of carbon, as will be discussed later in this section. 


Rayleigh Scattering Method 


This section describes the general theory of scattering by small particles 
such as carbon particles (in which a « 1) known as Rayleigh scattering. This 
type of scattering involves the radiation from an oscillating dipole as described 
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in the previous section. Among the many references which discuss Rayleigh 
scattering are Oster [243] , Cadle [244] , Green and Lane [245] , van de Hulst 
[246], Sinclair [247], Sinclair and La Mer [248], Dallavalle 1249], Orr and 
Dallavalle [250], Sommerfeld [236], Humphries [251], Tverskoi [252], and 
Bullrich [253]. 

According to the basic theory of Rayleigh [254-256], the oscillating 
electric field of the radiation incident upon a transparent optically isotropic 
particle in which a equals (2ir r^/X) « 1 induces an electric moment in the 
particle. The statement that the particle is optically isotropic implies that 
the polarizability and index of refraction are independent of direction in the 
particle. Since the particle acts like a linear electric oscillator, it does not 
radiate light along other directions, however, and the scattered light is plane- 
polarized if seen at 90° to the incident ray, whether the incident ray is polar- 
ized or unpolarized. This polarization of radiation scattered by a particle cloud 
(such as the exhaust of a hydrocarbon-fueled rocket engine) is seen in Figure 19. 



FIGURE 19. POLARIZATION OF RADIATION SCATTERED BY 

A PARTICLE CLOUD 

If the particles were anisotropic, the direction of the incident radiation's 
electric field E might not coincide with the direction of the induced dipole 
moment p. hi this case the radiation scattered at right angles would not be 
completely plane-polarized with vertical vibrations (as for an isotrpic particle) , 
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but would have a weak horizontal component. This isotropic versus anisotropic 
scattering phenomenon is also shown in Figure 19. 

When the incident radiation is plane-polarized on an isotropic particle, 
the intensity scattered at an angle ip can be written according to Sinclair [247] 
as 


I 
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(270) 


where R is the distance from the sphere (carbon particle) to the point of obser- 
vation, ip is the angle between the propagation direction of the scattered radiation 
and the electric vector in the direction of incident radiation, and V is the volume 
of the sphere equal to 4/3 tt rp. The scattered radiation is thus plane-polarized 
and has its electric vector lying in the plane formed by the electric vector in the 
diiection of incident radiation and the direction of propagation of the scattered 
radiation. 


When the incident ray is unpolarized, the intensity of radiation scattered 
in the direction 0, according to Sinclair [247], is 
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(1 + cos 2 9) = ij + i 2 


(271) 


where 0 is the angle between the direction of propagation of the incident and 
scattered radiation. This equation shows that the radiation scattered by the 
paiticle consists of two incoherent, plane -polarized components whose planes 
of polarization are mutually perpendicular. The component q (the quantity in 
parentheses multiplied by unity) has its vibrations (electric vector) at right 
angles to the plane ol observation. The component i 2 (the quantity in paren- 
theses multiplied by cos 2 0 ) has its vibrations parallel to the plane of observation. 

Figure 20 shows the angular distribution of intensity and polarization of 
the radiation scattered from a small isotropic scattering particle if the incident 
radiation is unpolarized. The solid line represents the total intensity according 
to equation (271) , and the dotted lines the intensities of the polarized components. 
Curve 2, which represents the cos 2 0 term in equation (271) , is the component 
which is horizontally polarized, and curve 1, which represents the unity term in 
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FIGURE 20. RAYLEIGH SCATTERING DIAGRAM 
SHOWING SCATTERED INTENSITY, I , AS A FUNCTION OF 0 FOR 

UNPOLARIZED INCIDENT RADIATION OF INTENSITY, I 

o 

equation (271) is the term which is vertically polarized. From this figure it 
is thus seen that light scattered at 90° is vertically plane-polarized. 

To obtain the scattering cross-section a s , which is necessary in pre- 
dicting emissivities of carbon particles, the methods of Stratton [232], 
Sommerfeld [236], Born [257], Slater and Frank [234], or van de Hulst [246] 
should be consulted. The electric field of a scattered wave may be written 
according to van de Hulst in terms of the distance R from a dipole (R assumed 
much greater than X) as 
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k 2 p sin 7 
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-ikR 
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(272) 


where y is the angle the scattered ray makes with the induced dipole moment 
p = (3 E q (where /3 is the polarizability which is equal to ex/E) and k is the 
wave number which is equal to 2n/X. The intensities of the incident and scatter- 
ed radiation, I G and I, respectively, may be taken from the time average of 
the Poynting vector, S' 1 ' , of the plane wave 

= V Re (E H ) , (273) 

2 6 <p 
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where is the average magnetic field corresponding to the current element 
in the dipole. These intensities, using Gaussian units, reduce to 


c 


2 c 



E 

and I = 0 
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=3 

0 

8 7 r 



(274) 


If the scattered intensity I is integrated over a sphere, the total energy 
scattered in all directions, according to van de Hulst, is 


W = k 4 c p 

U 


(275) 


The scattering cross-section <t s may now be defined according to Slater 
and Frank [234] as the area on which enough energy falls from the plane wave 
to equal the scattered intensity. Hence, by dividing equation (275) by I Q , the 
scattering cross-section is obtained as 
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For spherical particles small with respect to the wavelength (a « 1) , 
the polarizability j3 may be written according to the Lorentz- Lorenz equation 
as 
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where x is given by equation (246) . 

Hence, the scattering cross-section may be written as 
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showing an inverse wavelength to the fourth power dependence. 
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The fractional decrease of the incident radiation scattered in all direc- 
tions may now be written in the following form: 


I e = I 
o o 


-NpCTgi r 128 7 r 5 

e P = I„exp[- — ^ 


n - l 


n 2 + 2 


v]. 


(279) 


where N p is the concentration of scattering particles per unit volume and f is 
the path length through the medium. 


Another method for obtaining the scattering cross-section for Rayleigh 
scattering is sometimes used [79, 105]. In this manner a solution to the macro- 
scopic wave equation, 


8 2 E ± 9 2 E 

9x 2 " v 2 ’ (280) 


where the phase velocity v equals c/n, may be found as 


E E q exp £ i27rr(t - y )J - E Q exp £ -27rrn 2 x/c J exp j~i2 7 r^- ( t - n x x/c)J . 

(281) 

Since the intensity varies as the square of E, the intensity becomes 


I - I exp -47Tim 2 x/c = I exp(-N a x) , 
° L Jo p s 


(282) 


where x is the same quantity as i in equation (279) . 

With the introduction of equation (252) for n 2 (put in terms of frequency 
r) and the insertion of equation (87) for y, it is shown by these authors that 
for ( v - v Q ) » y the scattering cross-section becomes 


a 

s 



(283) 
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(284) 


When v » v the cross-section becomes 
o 


a 

s 


_ 8 TO 4 

3m 2 c 4 



87r e 4 A 4 , 
o_ 1 

3m 2 c 4 A 4 


This is the Rayleigh scattering equation which implies that the scattering 
varies as the inverse of the wavelength to the fourth power (see equation (278)) . 
It might be mentioned that for v Q « v the scattering cross-section reduces to 


( 7 

S 


_ 87re 4 
3 m 2 c 4 


(285) 


which is the Thompson scattering cross-section for free electrons, and shows 
no dependence upon wavelength. 

Equations (283) and (284) are actually the same form as equation (276) 
if it is remembered that the polarizability in equation (276) can be written as 
ex/E and equation (246) is used for x. In the next section these expressions 
for scattering cross-section are used to predict emissivities from carbon 
particles. 


Methods of Predicting Carbon Absorption Coefficients 

and Emissivities 

This section discusses several methods which have been used to predict 
absorption coefficients and emissivities of carbon particles in luminous flames, 
in certain types of carbon suspensions, and in exhausts of hydrocarbon-fueled 
rocket engines. The first general method to be discussed involves the theory 
of Rayleigh scattering discussed in the last section; the next method involves 
calculations made using Mie theory (which is discussed more thoroughly in 
the next section for A1 2 0 3 particles) ; the third method is concerned with the 
inverse wavelength method; and the final method involves various experimental 
measurements. 

1. Rayleigh Theory Method. Using the Rayleigh theory of predicting 
scattering and absorption cross-sections, Main and Bauer [258] have calculated 
absorption coefficients for particulate carbon in carbon-air mixtures. These 
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calculations were performed for particle radii of 50 A to 1600 A and for wave- 
lengths of 2000 A, 4000 A, and 8000 A. 

Using the Rayleigh (or Mie) theories, the absorption coefficient for a 
particle cloud may be written in the form of the total or extinction cross-section 
a as 


K 


A 


a. 



(286) 


where equals a g + <r a and cr a is the absorption cross-section. The cross- 
sections a s , a a , and are related to the efficiency factors Q g , Q a , and Qj., 
respectively, in the following manner: 
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cr 

s 

r 2 

P 




and 




(287) 


By using equation (276) for the scattering cross-section, the following 
relation for the scattering efficiency factor is obtained: 


Q 


s 



n 2 - 1 
n 2 + 2 


(288) 


where a = kr p = (2 tt/A) r p . 


The absorption efficiency factor, Q a , derived in 
the next section, may be written as 


Q a = - 4a Im ( 5fi ) 


(289) 


The concentration, N p , of carbon particles of mean radius r p , which 
makes up a dispersion of mass density p can be written according to Main and 
Bauer [258] as 


N 

P 



7rr 


P 



(290) 
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where 6 is the density of solid graphite. The absorption coefficient may thus 
be written in the form 


m 3 ,£-\ 1 

k. = a N = — (-f~) — 

A t p 4 6 r 


i t f « - 1 

-4a Im I 5 — r 
\ n 2 + 2 


n - 1 
n 2 + 2 


(291) 


However, for particle radii of 200 A or smaller (close to the median 
size of carbon particles in flames or rocket exhausts as discussed in the second 
part of this section) and for wavelengths greater than about 2 /it, Main and Bauer 
found that the scattering cross-sections were several orders of magnitude 
smaller than the absorption cross-sections. Hence, the absorption coefficient 


k reduces to 
X 


n 2 - 1 \ i 
n 2 + 2 ) J 1 


(292) 


This relation for absorption coefficient thus depends only upon the index 
of refraction n (and hence the wavelength of radiation, A) and is thus independ- 
ent of radii of the particles. Using values of lij and n 2 of graphite at 2250° K 
based upon the calculations of Stull and Plass [240] and the measurements of 
Senftleben and Benedict [242], Main and Bauer [258] calculated values of k ^ 
for a carbon-air suspension in which the amount of carbon was varied from 0 
to 100 percent by weight. These calculations were made for temperatures of 
3000° K to 10,000° K and for pressures of 0. 0316 atm and 1 atm. 


2. Mie Theory Method. The detailed derivations of the equations involved 
in the Mie scattering theory are discussed in the section entitled "Radiation from 
A1 2 0 3 Particles," where it is shown that the Rayleigh theory is merely a limiting 
case of the Mie theory. Although the Mie theory has usually been used for 
larger particles of 1 to 5 p in radius and in which the parameter a is approxi- 
mately 1, various investigators such as Stull and Plass [240], Krascella [259], 
Beheshti [260], Foster [261], Erickson, Williams, and Hottel [225], Howarth, 
Foster, and Thring [262], and Yagi [263] and Ruedy [264] (in two very early 
analyses) have used Mie theory for carbon particles in which a is generally 
assumed to be less than unity. 
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In Stull and Plass's [240] investigation, carbon particle emissivities 
were obtained for a range of particle radii from 0. 005 p to 0. 1 p and a wave- 
length range of 0. 4 p to 20 p. Krascella [259] considered particle radii ranging 
from 0. 05 p to 1 . 0 p and a wavelength range of 0. 25 p to 4. 0 p. Beheshti [260] 
investigated particle radii from 0. 005 p to 1, 0 p and wavelengths varying from 

0. 2 p to 40 p. Foster [261] treated particle radii from 0. 05 p to 0. 4 p and a 
wavelength range of 0. 5 p to 6. 0 p. 

Erickson, Williams, and Hottel [225] made both theoretical calculations 
and experimental measurements to determine intensity of radiation scattered 
from a premixed laminar benzene-air flame. In the experimental investigation, 
relative intensities for both the parallel and perpendicular components were 
measured at a number of angular positions in the range of scattering angle 
e = 30 to 150 deg. For the theoretical investigation using Mie theory, values of 
intensity were calculated for refractive indices n = 1. 71 - 0. 76i and n = 1.40 - 

1. OOi for values of a of 0. 8 to 1.30 and for wavelengths of 4358 A and 5461 A. 

Howarth, Foster, and Thring [262] evaluated the extinction coefficient 
Qt for carbon pai tides less than 600 A in diameter in the 1-p to 4-p wavelength 
range and over a temperature range of 1000° K to 2000° K. These authors 
calculated the damping constants y u y 2 , y 3 , and y c appearing in equations (265) 
and (266) as a function of temperature and found that in the 1-p to 4-p wave- 
length range and 1000° K - 2000° K temperature range the value of varied 
from 0 to 30 percent. It was found that the mean value of extinction coefficients 
measured over a non- isothermal optical path varied less than 10 percent. 

hi all of the above investigations mentioned, the carbon particles were 
assumed to be spherical and embedded in an infinite, homogeneous, non- 
conducting medium. The general equations for the Mie total and scattering 
cross-sections are given by Mie [265], Stratton [232], van de Hulst [246], 

Born and Wolfe [235], and are derived in the section entitled "Radiation from 
AI2O3 Particles, ” of this report as 

2?rr 2 <=o 

a t = rv 2 P He Yj (2m +1) (a +b ) (293) 

1 a m m ' ' 

m=l 

and 

271 r 2 <=o 

% = X < 2m + « <|vf + | b J 2 > • (294) 
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where the ^ and b m are the Mie coefficients which may be expressed in terms 
of spherical Bessel and Hankel functions, J m and hJ", respectively. 

The equations for a m and b m may be simplified if the carbon sphere 
and the external medium are assumed to have the same magnetic permeability 
(Mi = P 2 ) , if the recursion relation is used for J m and its derivatives, and if 
the Neumann function Y m - i (J m - H^) is used. The discussion of the above 
method to obtain Mie coefficients, together with another method, the logarithmic 
derivative method, is given in the following section. 

The eight groups of investigators mentioned [225, 240, and 259-264] 
found, as Main and Bauer [258] , using the Rayleigh theory, that the absorption 
cross-sections o a (and hence Q a ) were much larger than the scattering cross- 
sections a s (and hence Q s ) for low values of a. Figure 21 (taken from Foster 
[261]) shows this phenomenon at a wavelength of 1. 8 p. It will be shown in 
the following section that with further increase in a the total or extinction 
efficiency factor will oscillate about the value 2.0. 



FIGURE 21. EXTINCTION AND SCATTERING EFFICIENCY FACTORS 
FOR CARBON PARTICLES AS A FUNCTION OF a FOR A = 1. 8p 

After the scattering and absorption cross-sections have been obtained, 
the emissivity of a cloud of particles may be obtained as 


e = 1 - e 


-N d (<t, 


+ °a) 1 = 1 _ 


-N 


p 


(295) 
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Both Stull and Plass [240] and Beheshti [260] present curves of emissiv- 
ity for carbon particle clouds as a function of wavelength for various values of 
L = Npi and particle radii. Stull and Plass also present curves of intensity 
(discussed in the following subsection) and emissivity for two distributions of 
particle sizes. One of these was a symmetrical distribution of particle radii 
about 200 A and was calculated as 


N(r ) dr - (N /159. 52) exp < 
P P P ) 


j . 


(r - 200) /90 


1 2 


dr 


(296) 


where rp was in angstroms. Another distribution, which was asymmetrical, 
was calculated as 


N(r ) dr = 4. 75 x 10 5 N r " 3 e 640//r P dr (297) 

P P P p p ' ’ 

for 50 A < r < 1000 A and was assumed to be zero outside of this range. 

r 

Figure 22, taken from Stull and Plass [240] , shows emissivity as a 
function of wavelength using Mie theory and equation (297) as the particle 
distribution. It is apparent from this figure that scattering is important only 
for the small wavelengths (X < 2 ju) . For L = 10 13 particles/ cm 2 , it can be 
seen that emissivity is that of a blackbody over practically all of the infrared 
wavelength region. Thompson [266], who made radiation measurements in the 
exhaust of an F-l turbine, showed that, for a value of L of 2. 72 x 10 13 (cor- 
responding to a path length of 4 3/4 in.) , using equation (297) , an 800° K 
blackbody curve resulted for the emissivity. Figure 23 (taken from Beheshti 
[260]) shows emissivity as a function of wavelength for various particle radii 
and values of L. These values of rp and L are felt to be representative of 
exhausts of RP-1 fueled rocket engines. 

3. Inverse Wavelength Method. In addition to the Rayleigh and Mie 
theories of predicting emissivities of carbon particles, numerous investigators 
have used approximate relations based on an inverse wavelength dependence 
which is derived from experimental data. Some of these investigators include 
Becker [267], Rossler and Behrens [268], Pepperhoff [269], Naeser and 
Pepperhoff [270], Pepperhoff and Bahr [271], Plyler and Humphreys [272], 
Rossler [273], Heidman and Priem [274], Yagi and Lino [275], Hottel and 
Broughton [276], Millikan [227,228], Schack [277], and Siddall and McGrath 
[278]. 
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Emissivity 



FIGURE 22. CARBON PARTICLE EMISSIVITY AS A FUNCTION OF 
WAVELENGTH FOR THE PARTICLE SIZE DISTRIBUTION 
OF EQUATION (297) 



FIGURE 23. CARBON PARTICLE EMISSIVITY AS A FUNCTION OF 
WAVELENGTH FOR PARTICLES OF VARIOUS RADII 
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These empirical expressions for mass absorption coefficient x and 
spectral emissivity have usually been expressed in the following form: 

K m,A = CX_A and € A = 1 " eXp(K m,A Pi) = 1 " ex P(C X _A pi). 

(298) 

In these equations C is a constant independent of wavelength, but A is not 
necessarily independent of wavelength. The total emissivity is usually found 
from equation (298) by integrating over all wavelengths as 

°° A 

f[ 1 - exp(C A pi) ] B dA 
0 A 

, (299) 


where is the Planck blackbody intensity function defined in equation (19) . 

Hottel and Broughton [276], who analyzed various experimental data 
of soot deposited on glass, concluded that the exponent A could be represented 
by a value of 0. 95 in the infrared down to 0.8 In the visible part of the 
spectrum they recommended an average exponent of 1.39 for various amyl 
acetate, acetylene, and city -gas soots. Pepperhoff [269] and Rosslcr [273] 
state that values of A higher than unity may be attributed to the effect of par- 
ticle size. Millikan [228], however, states that for the particle diameter 
range 600 to 800 A the exponent A is independent of particle size, state of 
aggregation, or temperature. He comments that the exponent A is determined 
mainly by the hydrogen content of the sample, since he found that A increased 
from 0. 9 to 1. 9 as the hydrogen content varied from 0. 8 to 4. 1 percent by 
weight. 


Siddall and McGrath [278] summarize results of various experimental 
investigations involving determination of the exponent A. They state that, in 
the visible part of the spectrum, A can be assumed independent of wavelength, 
but in the infrared, it should be written in the form A = a - b In A. They also 
comment that the parameter A is independent of soot layer thickness but is 
dependent upon the fuel from which the soot is formed. 
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Sato and Matsumoto [279] calculate (and measure) radiation from a 
carbon particle cloud using equation (298) for spectral emissivity as 

e - 1 - exp(-C pt/\ A ) = 1 - exp[-(C 2 /TA A ) (CipT/C 2 )] = 1 - e 
A, 

p 

(300) 


where C 2 is given in equation (198) as ch/k. 

The total emissivity from a carbon particle cloud is then found as 


e 

P 





1 ) 


(301) 


where the quantity ip M, (X + 1) is defined in terms of the gamma function 

T(x + 1) as 

e -xl d|= -di- ln r(x+» =<T(x+i> • < 3 ° 2 > 

J o e« i d)< 


The exponent A was considered to be unity, and the mass absorption coefficient 

k was taken from Schack [277] to be 
m, A. 


m, A 


-CA 


5. 7 x 10 5 /A 


(303) 


where X is measured in microns and k ^ is in m The expression for total 
emissivity, e [equation (301)] was then used to predict radiation from a cloud 
of carbon particles, as will be discussed in the following section. 

Siddall and McGrath [278] and Thring et al. [280] discuss another 
expression for mass absorption coefficient K m ^ = CA which is based upon 
the Mie theory calculations of Hawksley [281] io obtain the total or extinction 
efficiency factor Q t . In this manner the absorption coefficient may be 
expressed as 
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(304) 


k — pK = Q N A 
A m,A t p 


where A is the projected area of a particle (equal to 7rdp/4) , and the total 
or extinction efficiency factor may be written as 


Q 


t 


24 7rd o 

£ 

A [ (n^ + nf n 2 ) 2 + 4(n? - n 2 ! n§ + 1)] 


24 7r d 

= F(X) 


(305) 


where the index of refraction n is now written as n = n x (l - in 2 ) . 
The absorption coefficient may now be written as 

36 7r V* F (A) 

K \~ A 


(306) 


where V* is the average volume of particles per unit volume of the cloud and 
may be defined as 


V' = N 7rd®/6 (307) 

r l -1 

where is the average number of particles per unit voliane. 

Ferriso et al. [230] used equation (306) in the following form 


- 367rpF(A)/p o A , (308) 

where p is the carbon particle mass density (g/cm 3 of cloud) and p Q is the 
density of carbon in bulk form, to plot the parameter x^A/p for values of 
temperature from 150 to 3000° K and values of wavelength from 1 to 10 p. 

This curve is reproduced as Figure 24 of this report. 
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FIGURE 24. LINEAR ABSORPTION COEFFICIENT FOR CLOUDS 
OF VERY SMALL CARBON PARTICLES 

To evaluate the parameters and n 2 (originally given in equation 
(268) and (269)) , Ferriso et al. [230] made the following approximations for 
wavelengths greater than about 1 y. 

n 2 - n 2 n| « 4. 36 - 1. 31 ( 309 ) 

and 

, 0.58 1. 675 A 

n l n 2 w “ + p (310) 

where 

P = a(2250° K)/a(T) 

6 = N (2250° K)/N (T) ( 311 > 

c c 

a(T) = e 2 N /me yy 
' c o c 
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where A. is the wavelength in microns, <j is the dc conductivity [ equation (260)] 
evaluated according to Stull and Plass [240] based upon measurements of 
Senftleben and Benedict [242], and N c is the number of conduction electrons 
per cm 3 . 

4. Experimental Methods . Various experimental methods of obtaining 
absorption coefficients, emissivities, reflectivities or transmissivities of 
carbon samples or carbon particle clouds are discussed in this section. This 
discussion is in addition to the brief one in the last section in which experimental 
data were mentioned in connection with the inverse wavelength calculation of 
absorption coefficients and absorptivities [equations (298) and (299)]. In the 
last 50 years or so, there have been numerous experimental programs per- 
formed to determine the emissivity of carbon or graphite samples placed in 
furnaces or other devices, carbon particles placed in various suspensions such 
as H 2 0, He, or N 2 , and of carbon particles in luminous flames. Only in the 
last few years, however, has there been an attempt to determine experimentally 
carbon emissivities and absorption coefficients from rocket exhaust gases. 

Plunkett and Kingery [282] , Rutgers [283] , and MacPherson [284] 
review the literature for emissive and reflective properties of carbon or 
graphite samples (polished or unpolished) placed in such devices as carbon 
tube furnaces and heated by electrical current or electrical resistance heaters. 
These investigators state that spectral emissivity generally decreases with 
increasing wavelength above 1 p; however, in the visible region the spectral 
emissivity closely approaches that of a graybody (i. e. , is essentially con- 
stant for all visible wavelengths) . These authors also comment that there 
appears to be no temperature dependence on the spectral emissivity, while the 
total or integrated emissivity appears to have a positive temperature dependence. 

Lanzo [285] and Lanzo and Ragsdale [286] discuss transmissivity 
measurements made on four carbon powders of particle diameters 0.08 p, 

0. 15 p, 0.45 p, and 1.40 p suspended in water for the wavelength range 0.2 p 
to 1.0 p. Spectral transmissivity measurements were obtained with a spectro- 
photometer, and total transmissivity measurements were obtained with a 
xenon-arc source and a radiation detector. It was found that transmissivity 
(and hence emissivity) was essentially independent of wavelength in the visible 
range, thus agreeing with Plunkett and Kengery [282], Rutgers [283], and 
MacPherson [284]. 

Marteney [287] measured extinction (absorption plus scattering) of 
carbon particles of size r p = 0. 0045 p suspended in He and N 2 carrier gases 
in the wavelength range 0.255 p - 0. 545 p. Extinction coefficients were again 
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found to be essentially constant with wavelength in the visible range. It was 
found that the application of de-agglomeration aerodynamic shear forces in a 
nozzle through which the carrier gas-particle mixture was passed caused an 
increase by a factor of nearly six in the extinction coefficient for the carbon 
particles. 

Numerous investigators such as Millikan [227, 228] , Plyler and 
Humphries [272], Silverman [288], Tourin [289], and Sato and Matsumoto 
[279] have made experimental measurements of emissivity in luminous flames. 
In the last several years there has been a considerable amount of literature 
involving the experimental determination of emissivities, spectral radiance, 
and temperatures in luminous rocket exhausts. Some of the unclassified refer- 
ences in this area include Rossler [273], Simmons and de Bell [290-292], 
Hergetetal. [293], Boynton et al. [231], Ferriso et al. [230], and Hendershot 
and McCaa [294, 295] . Some of the classified references in this area (per- 
taining mainly to H 1 and F-l engines) include Thompson [296], Simmons [152, 
297], de Bell and Speiser [298], de Bell and Simmons [299], Herget, 
Schumacher, and Enloe [300], Wagner [156], Levin, Wagner, and Thompson 
[154], and Tourin [301]. 

The measurement of emissivity for both flames and rocket exhaust gases 
depends upon the principles of monochromatic gas pyrometry as discussed by 
Silverman [288], Hecht [302], Hill [303], Hornbeck [304], Brenden [305], 
Branstetter [306], Tourin [307, 308], Brewer and McGregor [309], Harrison 
[310] , and de Bell et al. [153] . In this manner the intensity I F , as seen by 

the radiometer or spectrometer viewing the flame or rocket exhaust may be 
written in the Planckian form as 


2C t 1 

V " e F A X [exp(C 2 /X T ) - 1] 

A A r 



(312) 


and the intensity Iq of the continuous source (graybody) also seen by the 
radiometer or spectrometer is 


2 Ci 1 

" € G X 5 [exp(C 2 /XT G ) - 1] 



(313) 
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However, when a monochromatic beam of radiation of known intensity 
I G passes through a gas, part of it is absorbed by the gas such that only a 

fraction, I-p^, is transmitted through the gas and received by the detector. 

The fraction lost is thus the spectral absorptivity of the flame, a jy which, 

by Kirchoff's law, must equal the spectral emissivity, Cjy These quantities 
are expressed as 


a 




(314) 


This spectral emissivity may thus be measured as a function of wave 
length for a known value of for the graybody source. In determining the 

flame temperature, T F , as seen by the detector, equation (312) is divided by 
equation (313) , and the following equation results: 


T = 
F 


A In 


{ 1 + (W \) ( 1 ' W) IexptC!AT(3) ' 1 1 } 


(315) 


Ferriso et al. [230] and Boynton et al. [231] used this equation to 
determine the gas temperatures of a small RP- 1/gaseous oxygen engine ex- 
haust at 4. 4 m, the CO 2 and CO vibration-rotation bands, and at 2. 2 p (region 
of no molecular emission) to obtain the carbon temperature. 


Simmons and de Bell [290,291] discuss the use of two path pyrometry 
to determine simultaneously the spectral emissivity and temperature of a 
luminous rocket exhaust. The main advantage to this method is that a com- 
parison source behind the rocket exhaust is not needed. Instead, a corner- 
cube reflector (glass prism) is placed behind the exhaust in effect to double 
the apparent depth of the plume. The radiant intensity thus obtained when a 
reflector is placed behind the plume can be written as 




(316) 


119 


where is the reflectivity corresponding to the wavelength of the 
Assuming thermal equilibrium, e F = I F , and the absorptivity 

the emissivity, which is then written as 


reflector. 

a equals 
* 


e 


= 1 - 




(317) 


Assuming that the Wien law for blackbody radiation is sufficiently 
accurate, the temperature at the nozzle exit may be calculated as 
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where Tbr is the single-path brightness temperature at the nozzle exit and 
is thus the ratio of the double -path spectral intensity to the single - 

path spectral intensity. 

Most of the spectral and total emissivity data from the exhausts of the 
1,500,000-lb thrust F-l engine and the 200,000-lb thrust H-l engine are 
classified. However, most engineers who have calculated radiative heating 
to the base of the Saturn vehicles are inclined to use a value close to unity 
for the emissivity in the infrared wavelength region. Simmons [297] devel- 
oped an empirical relation (which is classified) for spectral emissivity in the 
visible spectrum for an RP-l/lox rocket exhaust. This empirical relation 
was a function of wavelength, mixture ratio, chamber pressure, and thrust 
chamber size. 

Perhaps the most extensive experimental analysis of emissivity from 
carbon particles has been determined by Ferriso et al. [230] and Boynton 
et al. [231] at General Dynamics/Convair under NASA/MSFC contract. In 
this investigation, the dependence of the carbon spectral absorption coefficient 
on wavelength in the 0. 7-p to 5-p range and on temperature from 1045° K to 
2600 K was determined. Three small contoured-nozzle rocket motors of 
nozzle area ratio 5. 25, 3. 0, and 1. 5 were used in the investigation. The 
motors used the Foelsch [311] nozzle to produce a constant temperature near 
the exit. By varying the chamber pressure and, hence, O/F ratio, various 
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values of constant temperature were obtained for each motor. These tempera- 
tures were 1045° K, 1290° K, 1415° K, and 1680° K for the 5.25 area ratio 
motor; 1230° K, 1470° K, and 1760° K for the 3. 0 area ratio motor; and 
1700° K, 1900° K, 2370° K, and 2600° K for the 1. 5 area ratio motor. 

All of the above temperatures were calculated from a relation such as 
equation (315) . The emissivity of the entire plume, was obtained from 
experimental measurements as shown inequation (314). This emissivity can 
be expressed in terms of carbon, H 2 0, and C0 2 absorption coefficients as 


e 
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- exp 


ru ,p) i 

m,A c 


( k . u ) + 

m,A H 2 0 


(k . u) 
m,A CO: 


>]} 


(319) 


where u represents the quantity pH taken at standard temperature and pressure. 
The contribution to the emissivity of CO and other molecular emitters is 
neglected. 


The equation for the quantity (^ m , A j-j q usin § ^ ie statistical band 
model, was given by Ferriso et al. [230] (see equation (227)) as 


(k u) _ = uk (1 + uk /4a) 
m, A H 2 0 m, A m, A 


(320) 


where a represents a fine structure term represented by 


a = 





(1 


h 2 o' 



(321) 


In this equation, a is the ratio of half-widths due to self -broadening, a x is the 
ratio of foreign gas half- width to H 2 0 half-width, and C^ is the concentration 

of water vapor. For C0 2 , the average absorption coefficient can be determined 
by the Beer-Lambert equation. In this case, equals S/d, where S/ d was 
given as a function of temperature and wavelength by Malkmus [312] . 
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Figure 25, taken from Ferriso et al. [230], shows the quantity (K~ m ^p) c 
for the 5.25 area ratio motor. The calculated curves are based upon bulk 
carbon properties as determined by equations (308-311) . Good agreement with 
experimental data was obtained up to 1700° K, after which the experimental 
data were practically independent of wavelength. At temperatures below 1200° K 
and at wavelengths less than 0. 8 p, the observed emitted intensity was greater 
than the blackbody intensity. However, it was felt that the effects of scattering 
or a non-equilibrium distribution of non-emitting states in the rocket exhaust 
could have contributed to this phenomenon. 


Radiant Heating Calculations for Luminous Flames 
and Rocket Exhausts 


This section discusses various methods of predicting radiant heating 
based upon the emissivity relations of the preceding section. The general 
methods are applicable to both luminous flames or rocket exhausts in which 
carbon is the principal emitter. The problem of thermal equilibrium between 
carbon particles and the gas is also described. The base radiant thermal 
environment used as design criteria for the Saturn I, S-I stage, and Saturn IB 
(now Uprated Saturn I) , S-IB stage caused by H-l rocket exhausts and for the 
Saturn V, S-IC stage caused by F-l rocket exhausts is also discussed. 

Stull and Plass [240] present an equation for the radiant intensity from 
a distribution of carbon particles which absorb, scatter, and emit radiation 
isotropically and which are at a constant temperature equal to that of the gas. 
Integration of this equation over wavelength dA and solid angle do) would thus 
yield the radiant heating rate per area. This equation, which includes higher- 
order scattering terms, is 


z 

I (0 ,<f> ) = B (T) N cr f maX exp(-N a z sec 0 ) 
A o ,r o A p ad p t o 


S 1 + ("^r) ex Pt'N p a t R') (322) 

j 1 + (“ir) /,5772 e X p(-N /t R") { l + ...}dV"j dV* j ds , 
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FIGURE 25. VALUES OF (k p) FOR CARBON PARTICLES 

m,A c 

DETERMINED FROM MEASUREMENTS ON THE 
GD/C 5. 25:1 AREA RATIO MOTOR 
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where 0 O and <p Q are the polar and azimuthal angles, respectively, from the 
element of area being heated to the emitting carbon particles in the volume 
element dV, and R is the path length from the element of area to the volume 
dV. The element of volume dV is R 2 du>dR, where R = zsec0 o and du> is the 
element of solid angle. The dV 1 terms involve the amount of radiation emitted 
in the volume dV' scattered within a volume dV which is separated from the 
volume dV' by the path length R'. In this case the element of volume dV is 
dR 1 dA, where A = R ,2 du>. Also in this equation, the quantity cr a N p is equiva- 
lent to the quantity k x = pK m> \ in equation (207) , and cr s N p is the scattering- 
coefficient a x = pa m x defined in equation (29) . 

For the case in which scattering by the carbon particles may be neg- 
lected (i. e. , for r p < 300 A and A > 2 p) , the emitted intensity may be written 
by taking only the first term in equation (322) (a s = 0 in the other terms) as 

a , x 

I (fl ) = B (T) — cose 1 - exp [<7 L(0 <p )] , (323) 

A o A a. o ( too) 


where 
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(324) 


The quantity z max represents the total path length projected in the 
direction z, where z is measured in a direction at right angles to the front 
surface of the flame or rocket exhaust. 

Stull and Plass [240] calculated spectral intensities for carbon particles 
with radii equal to 200 A, and for the two distributions, equations (296) and 
(297) . They found that the maximum intensity for these carbon particles 
occurred at smaller wavelengths than the corresponding maximum of the black- 
body intensity curve. The intensity was found to increase as the number of 
particles increased until the limiting blackbody intensity was reached. For 
particle radii of 200 A, about the average size reported by most investigators 
(See this section, subsection entitled "Carbon Size Particle Determination. ") , 
the intensity was practically negligible when there were less than 10 11 particles/ 
cm 2 , and closely approached that of a blackbody when there were more than 
10 14 particles/cm 2 . This characteristic can also be seen in the emissivity 
curve (Figure 19) for the particle distribution of equation (297) . 
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Sato and Matsumoto [279] present another radiant heating equation 
which is based upon the inverse wavelength dependence of the carbon coefficient. 
This equation holds for luminous particles imbedded in a gas; however, the 
effects of scattering are neglected. It is assumed that the nonluminous gases, 
such as H 2 0 and CO 2 , are uniformly distributed in the flame. Furthermore, it 
is assumed that the gas temperature equals the temperature of the carbon 
particles. The expression for the spectral emissivity may thus be written as 


= 1 - exp { -l<’‘ m> A P V U m,X P, g l< } 


(325) 


where the subscript c stands for carbon and the subscript g stands for non- 
luminous gas. 


The radiant intensity may be written as in equation (63) for no scatter- 


ing as 


\ = £ A B A (T) = 


— 7 -^ 7 [ 1 - exp i - [ (k p) + (k p) 

exp(C 2 /A T) -1 \ m,X" c m,A g 

(326) 


By integrating over wavelength dA and solid angle dc o> the expression 
for radiant heating rate per unit area becomes 


-5 


^ A= // exp(C 2 AT)TT - eX P{- [ ( K m ,X P) C +( Vx P) g ] 1 ) dWdX 

(327) 


A u) 


By using the variables 4 = C 2 /A T and x - pT/C 2 as defined in 
equation (300) and using the symbol rj for («xP\) 4, equation (327) becomes 


Ci 

q/A = 

C 2 


h) |_ 0 


s 3 -(xl + n) 

e«-l 


di- 


do) . 


(328) 


By making the substitution of equation (302) and using the approximation 
e w 1 - €g, equation (328) now becomes 
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(329) 
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Tuttle [313] presents another method of calculating radiant heating from 

carbon particles in a luminous gas. This particular calculation was for the 

exit plane of an H-l engine on the Atlas vehicle. A form of equation (207) was 

used to obtain the spectral radiance, I , as 

A 

S r 

r m r r s 

I = f (pK ) B (T) exp - f (PK ) ds 
A m , A c A L 0 m> A c 

(330) 



The dimensionless optical depth may thus be written as 
s m s m 

T x‘j <pK m,X ) c dS = / 'Va’c ds • (331 > 


For the temperature range considered in this analysis, Tuttle states 
*m,A is independent of temperature, and hence is not a function of path 
length in the gas. Therefore, the optical depth may be written as 


T = (p K ) f 
A m m, A e J 
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P C (S) 
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ds 


(332) 


where p m ^ is the mean particle density or the average number of carbon 
particles/cm 3 . 

If the carbon is assumed to be amorphous with a specific gravity of 1.8, 
a radius of 500 A, and a ratio of the mass of carbon to the total mass flow of 
0. 005, the mean particle density may be calculated to be 5. 0 x 10 8 particles/cm 3 . 
Tuttle presents a temperature distribution curve based upon Simmons' [297] 
data as 
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where a A , b A , and c A are constants tabulated by Tuttle [313] and r c repre- 

A 

sents ( the optical depth at the center of the plume. Inserting these values of 
“ lto equation (330) results in the following equation: 


\ ■ V 1 
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+ (c A + 1) c' 2Tc a] . 


(335) 


This equation was compared with the spectral data of deBell and Speiser 
[298] in the 0. 5-p to 3-p wavelength range. Very good agreement with the 
experimental data was obtained except at low wavelengths where scattering 
effects should be included and in the 1.4-p, 1. 9-p, and 2.7 -p bands in which 
atmospheric absorption was present. Tuttle [313] also included expressions 
for upper and lower limits of radiation from the cone of the H-l engine. The 
upper limit took into account the contribution of the H 2 0 emission, and the 
lower limit neglected this emission. 

Several new investigations of radiation from luminous flames and rocket 
exhausts have appeared in the literature recently. Thring, Beer, and Foster 
1 314] discussed calculations of carbon emissivities based on Mie theory and 
also measurements of soot concentration from hydrocarbon-air turbulent 
diffusion flames. These authors found that in certain cases the measured 
emissivity was 2 to 3 times larger than the calculated value; however, they 
believed this to be because the refractive index of soot containing hydrogen was 
possibly very different from that of pure carbon. Wolfhard and Hinck [315] 
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discussed the excitation (non-equilibrium radiation) of OH molecules and alkali 
metals in luminous flames and rocket exhausts. Zirkind [316] presented a 
general review of radiation processes in rocket exhausts, discussing band 
models, scattering processes, and spectral profiles. 

In all of the radiation analyses mentioned in this section, it has been 
assumed that the carbon particles were in thermal equilibrium with the gas 
(i. e. , T p = T g ) . Simmons and Spadaro [317, 318] have investigated this 
thermal equilibrium assumption by determining the amount of thermal lag of 
carbon particles in rocket nozzle flow. These authors set up a heat balance 
equation for a single carbon particle consisting of (1) the rate of heat transfer 
by collision with the gas molecules, (2) the rate of radiant heat transfer to the 
external surroundings, and (3) the rate of change of the heat content of the 
particle. 

Based upon the kinetic theory of gases and the thermal accommodation 
coefficient, a, as discussed by Wiedmann and Trumpler [319] , the heat trans- 
fer in time dt between the carbon particle and the gas, dQ t may be written 
according to Dillon and Line [320] as 



where (47r r 2 ) is the surface area of the spherical carbon particle; C y , MW, 
R, P, and 'f'g are the specific heat at constant volume, molecular weight, gas 
constant, pressure, and temperature, respectively, of the gas; and Tp is the 
temperature of the carbon particle. 

The radiant heating emitted per time from a single carbon particle may 
be expressed in the form 


dQ 2 = (4? r r^) <T fi £ p (1 - a g ) dt , (337) 

where eg is the Stefan— Boltzmann constant, Cp the total hemispherical 
emissivity of the particle, and ctg is the average total absorptivity of the gas 
stream which is defined in terms of the Beer-Lambert equation as 


128 


( 338 ) 


a ^l-e-^i-e-^p/ 4 * V 

g 

In this equation the quantity i represents the average path length (radius of 
the stream) , and Np is the concentration of the particles in the gas. 

The heat content of the particle may be written as 

dT 

dQ 3 - (-J *r 3 ) p p C dt , (339) 


where pp and C are the density and specific heat, respectively, of the particle. 

Combining these heat inputs (dQ 1# dQ 2 , and dQ 3 ) the following equation 
is obtained: 
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Simmons and Spadero [317] solved this equation numerically for T p 
as a function of axial distance x down the nozzle for three different rocket 
nozzles. This was done for accommodation coefficients, a, of 0. 5, 0.75, 
and 1.0, particle diameters 10 3 , 10 4 , 10 -5 cm, and for total gas emissivity 
of zero and unity. Figure 26, taken from Simmons and Spadero ( shows that 
the higher the accommodation coefficient, the closer the particle temperature 
was to the gas temperature, and that for carbon particles of 500 A or less in 
radius the particle and gas temperatures were practically the same. 

A recent investigation at Rocketdyne [199] using the Abel inversion/ 
zone radiometry technique as discussed by Herget et al. [293] has shown that 
it is possible for the temperature of the carbon particles in the region between 
the plume core and the afterburning mantle to be actually several hundred 
degrees Kelvin cooler than the gas temperature. When the carbon particles 
pass into the afterburning region, however, the temperature was found to be 
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FIGURE 26. VARIATION IN CARBON PARTICLE TEMPERATURES 
ALONG THE AXIS OF THE ROCKETDYNE 2:1:8 AREA RATIO NOZZLE 


equal to or slightly greater than the gas temperature. From the same investi- 
gation it was speculated that carbon particles may be formed just downstream 
of the first normal shock (in a fuel-rich region) in the plume, where the 
particles are again hotter than the gas. 

In the last several years there have been numerous NASA/MSFC, 
Chrysler, and Boeing reports published that describe the thermal environment 
to the base of the S-I, S-IB, and S-IC stages, respectively. Because of the 
uncertainties in obtaining accurate spectral absorption coefficients for carbon 
particles, the emissivities of the H-l and F-l plumes have usually been 
assumed to be those of a blackbody (i. e. , equal to unity) , and scattering effects 
have usually been neglected. Plume shapes and properties have been calculated 
by the methods described by Farmer et al. [1], Prozan [150], and Ratliff [151], 
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The results of scale model data (from both short- and long-duration tests) and 
of all the in-flight radiation measurements on the unmanned Saturn flights have 
been extrapolated to the case of the manned and unmanned Saturn IB and V 
flights. It has been shown that radiation to the base of these vehicles is a 
maximum at sea level and drops off considerably at altitude because of the 
expansion of the plumes. 

For the Saturn I vehicle, a summary of the base thermal environment 
caused by heating from the H-l rocket engine exhausts has been presented by 
Payne and Jones [9], Jones [321], and by Hartley and Fricken [322]. The 
Saturn IB base thermal environment caused by heating from the H-l engine 
exhausts is discussed by Francis [323], Fricken [324], Hartley [325-327], 
Fricken [328], and by Hartley and Fricken [322] . Most of these analyses 
assumed that the carbon density was five percent of the gas density (a . value 
representing an upper limit as obtained from Rocketdyne data) . Geraghty 
[329] presented a computer program for obtaining the form factors for the 
S-I and S-IB plumes. Figure 27, taken from Hartley and Fricken [322] , shows 
the design criteria for incident radiation, total heating, recovery temperature, 
and wall temperature for the S-IB flame shield, the hottest part of the base of 
this stage for all vehicles up to vehicle AS-202. All vehicles after AS-202 had 
the inboard engine turbine exhaust dumped in this region, thus cooling the flame 
shield considerably. The general configuration of the base region of the Saturn 
IB vehicle S-IB stage is shown in Figure 28. 

Mullen [330] and McEntire et al. [331] presented an analysis of the 
Saturn V S-IC stage base heating environment based upon plume temperature 
models obtained from radiant heating measurements taken during F-l engine 
static firings. Hughes and Reid [332] presented a discussion of the radiation 
data obtained from these F-l full scale static (sea level) firings. Jacobs [333] 
and Wasko [334] discussed the heating data (convection as well as radiation) 
obtained from S-IC scale model altitude tests with external flow using the 
Cornell short-duration shock tube technique. Patrick [335] discussed the 
radiation and total heating data obtained from exhausts of l/20th scale F-l 
engine firings at MSFC's Test Lab. Figure 29, taken from Mullen [330] 
shows the theoretical axial distribution of the sea level radiation heating on 
the base of the S-IC for various radial positions. Figure 30 shows the overall 
configuration of the base region of the S-IC stage. 

The Thermal Environment Branch is continuing to improve its base 
heating analyses of radiant heating from H-l and F-l rocket engine exhausts. 
Radiation heating rates are currently being calculated for the launch umbilical 
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FIGURE 27. FLAME SHIELD DESIGN LEVELS FOR INCIDENT RADIATION 
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FIGURE 28. SCHEMATIC OF THE S-IB BASE REGION 

tower and various components near the launch deflector at Pad 39 at Cape 
Kennedy which will be exposed to heating during lift-off of the first Saturn V 
vehicle. Radiometer and radiation calorimeter measurements were made in 
November 1966 near the flame bucket for one of the static firings of the 
S-IC stage at the MSFC Test Laboratory. Figure 31 shows the 200- to 300-foot 
luminous plumes of the S-IC stage during an MSFC static firing. Finally , more 
accurate predictions of radiant heating from carbon particles are currently 
being obtained with the incorporation of the General Dynamics/Convair 
[230, 231] carbon absorption coefficient program. This program is 
currently being used at Chrysler in an attempt to match the onboard spec 
trometer data obtained in the base region of the S-IB stage for the AS-203 
flight. The description of the radiation computer program which describes 
the use of the combined carbon, H 2 0, C0 2 , and CO absorption coefficients for 
a single, axisymmetric plume has been written recently by Dash . 


*M. J. Dash, NASA/MSFC Memorandum to be published, 1967. 
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FIGURE 29. S-IC RADIATIVE BASE HEATING FROM F-l ENGINES 

AT SEA LEVEL 

This concludes the section on radiant heating from carbon particles. 
The general theory of carbon particle radiation, including formation, size, 
scattering, and emissive properties of carbon particles and discussions of 
radiation calculations have been presented. Although there is considerable 
room for improvement in the general theory, it is felt that the present state 
of the art has enabled MSFC and its associated contractors to predict 
satisfactorily radiation thermal environments caused by carbon emission 
from the Saturn I, S-IB, and S-IC stage engine exhausts. 
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FIGURE 30. SCHEMATIC OF TIIE S-IC BASE REGION 


FIGURE 31. S-IC STATIC FIRING 




RADIATION FROM Al 2 0 3 PARTICLES 


Thermal radiation from the exhausts of aluminized solid propellant rocket 
motors is characterized by graybody (continuum) emission from the AI 2 O 3 
particles and by A1 2 0 3 particle temperatures which are larger than the gaseous 
temperatures at the same point in the plume. The metallic particles in the 
propellant, which may be magnesium, beryllium, or iron, as well as aluminum, 
are added to the propellant to increase its specific impulse and consequently 
its performance. However, the corresponding thermal radiation from the 
metallic oxides in the exhaust may be several times higher than the correspond- 
ing gaseous radiation emitted in the various molecular bands. Excessive 
radiative heating from exhausts of solid propellant motors may thus damage por- 
tions of the base region if the motor (s) are used as the main source of pro- 
pulsion, or may damage certain adjacent components and structures if the 
motor is used as an ullage or retro motor for stage separation purposes (such 
as for the Saturn vehicles) . The problem of microwave attenuation through 
the exhaust of a solid propellant motor is also more critical than through the 
exhaust of a liquid propellant engine of comparable size and performance. How- 
ever, this problem is not treated in this report. 

In this section the radiative properties of A1 2 0 3 particles are surveyed 
and analyzed in detail. It may be mentioned that MgO (magnesia) , BeO 
(beryllia) , and Fe 2 0 3 (ferric oxide) also act as emitters in exhausts of certain 
types of solid propellant engines. However, A1 2 0 3 (alumina) is found in the 
exhausts of the majority of metallized propellants in existence at the present 
time, and since the radiative properties of A1 2 0 3 are more fully understood 
than those of the other oxides, this section is devoted almost exclusively to 
radiation from A1 2 0 3 particles. 

It has been shown by many investigators such as Talbert et al. [336] , 
Dimmock and Courtney [337], Wagner, Cramer, and Borough [338], Briscoe, 
Bullara, and Bressler [339], Sutton [340], Schindler and Penzias [341], 

Miller [3421, Baker and Allport [343], Carlson [344], Carlson and Du Puis 
[345] , Launstein et al. [346] , two Aerojet-General reports [ 155, 347] , Miller 
and Sternchak [ 348] , Lai and Purgalis [ 35] , a Bamirac report [ 349] , and 
Condron [350] that the spectrum of the exhausts of solid-fueled engines follows 
a graybody continuum from about 1 p to around 3 p. It is known that the 
exhausts of double-base propellants of the nitroglycerine-nitrocellulose type, 
such as the Minuteman third stage, are similar to the carbon spectrums as 
discussed in the previous section. The double-base propellants generally 
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produce more carbon in the exhaust than do the composite solid propellants 
which are composed of a fuel, binder, and a finely ground oxidizer such as 

nh 4 cio 4 . 

The composite propellants such as those in the Minuteman first and 
second stages, Polaris, UTC-120 in. solids, and Saturn ullage and retro motors 
produce a continuous spectrum (on which gaseous band emission is superimposed) 
with a lower emissivity but with a higher particle temperature than do the double- 
base propellants. The composite propellants also produce an emission band of 
low intensity for the HC1 molecule around 3. 5 p. It is believed that the con- 
tinuum in composite propellant exhausts is caused mainly by A1 2 0 3 emission, 
although it is possible that carbon could also be emitting in the 1 - 3-p region. 
Upon analysis of some of the references mentioned in the above paragraph, it 
was found by Carlson et al. [351] that when the A1 2 0 3 mass fraction was approxi- 
mately doubled, the ratio of continuum to gas emission was approximately 
doubled, even though the carbon mass fraction was the same. This implies 
that the A1 2 0 3 was the major contributor to the continuum rather than the carbon. 
Also, spectrometer measurements at Aeronutronic [344] of the plume of an 
H 2 - 02 -Al 2 0 3 water slurry motor confirm the presence of an A1 2 0 3 continuum 
for 1-3 p. 

This A1 2 0 3 continuum emission in the 1. 9 -m to 2. 6 -p region (where 
there are no molecular emitters) may be seen in Figure 32, which shows the 
spectrum of the exhausts of scale models of the S-II ullage and S-IC retro motors 
fired at Cornell at 100,000 feet. This continuum would have existed down to 
1 ju or shorter, but the filter on the spectrometer cut off around 1.7 m. The 
higher temperature of the A1 2 0 3 particles (higher than that of the carbon for 
the same gas temperature) produces a very luminous core of particles, the 
larger particles of the order of 3 to 5 m tending to hug the plume axis and the 
smaller particles of the order of 1 m less tending to flare out beyond the gas 
plume. Some of these luminous solid propellant plumes are shown in Figures 
33, 34, and 35, which are respectively, the plume of the Minuteman first stage 
motor fired at sea level, the plume of the S-II ullage rocket fired at 121,000 
feet in the AEDC J-4 test cell, and the plume (Fe^ instead of A1 2 0 3 ) of a 
scale model of the S-IB retrorocket fired at 200,000 feet in an altitude cell at 
Cornell. 

It is possible that another type of radiation from solid propellant motors, 
the so-called "searchlight" radiation, may occur. This type of radiation involves 
liquid particles radiating at very high temperatures in the combustion chamber 
through the nozzle into the particle plume where the radiation is then scattered 
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FIGURE 32. RADIATIVE INTENSITY FROM EXHAUSTS OF S-II ULLAGE 
AND S-IC RETRO MOTORS FIRED AT CAL 


out of the plume. The radiation takes place in the visible or ultraviolet region 
(A < 0.8 m) because of the high temperatures, resulting in a blackbody curve 
which has its maximum in the visible or ultraviolet. The maximum radiation 
emitted would thus be equal to agT^Aj., where T c is the combustion chamber 
temperature and A t is the throat area. This searchlight effect has been shown 
by Carlson et al. [352] to be an important factor with small nozzles or light 
aluminum loadings. Recent investigations at Rocketdyne [199] have indicated 
that this searchlight effect can also occur in hydrocarbon-fueled engines in the 
ultraviolet and visible regions. 

Calculations of thermal radiation from particle plumes (such as A1 2 0 3 
plumes) are difficult because the particles can scatter radiation considerably 
as well as emit radiation. The Mie theory of scattering for these large AL0 3 
particles, in which the parameter a = 2 7r Tp/A is of the order of unity, must 
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FIGURE 33. FIRST STAGE M1NUTEMAN FLUME AT SEA LEVEL 


be nsec! to obtain scattering and absorption cross-sections, a s and a a , 
respectively (or efficiency factors Q s and Q a , respectively) from which 
the optical thickness, t, of the particle plume may be obtained. The cross- 
sections depend on the scattering angle, 6, on the complex index of refraction 
of A1 2 0 3 , n = nj - in 2 (which depends on temperature of the particles and wave- 
length of radiation) , and on the A1 2 0 3 particle size distribution. The cmissivity 
of the particle cloud is then obtained from the ratio of absorption to extinction 
efficiency factors, Q a /Q^, and the optical thickness, which is a function of 
particle concentration, N , and path length, i, as well as of the extinction 




FIGURE 34. S-II ULLAGE ROCKET PLUME AT ~ 121,000 FEET 
IN AEDC J-4 TEST CELL 

cross-section or efficiency factor. A two-phase plume program is used to 
compute the particle concentration, Np, for use in the optical thickness ex- 
pression and also the particle temperature, Tp. In this manner, the radiation 
from an A1 2 0 3 particle cloud may be determined with a reasonable degree of 
accuracy, depending upon various assumptions made, such as particle size and 
distribution. 

This section discusses the Mie theory of scattering, the determination 
(both experimentally and theoretically) of A1 2 0 3 particle sizes and distributions 
and the determination of A1 2 0 3 particle cloud emisssivities (including effects of 


140 



FIGURE 35. S-IB RETROROCKET PLUME AT ~ 200, 000 FEET 
IN CAL ALTITUDE CHAMBER 

anisotropic and multiple scattering) . The determination of A1 2 0 3 particle 
temperatures (two-phase flow phenomena) is then briefly treated, and finally 
various procedures for calculating radiation from aluminized rocket exhausts 
are analyzed. 


Mie Theory of Scattering 


The Mie theory of scattering of plane, monochromatic electromagnetic 
waves by spherical homogeneous conducting particles (in this case A1 2 0 3 
particles) in which O' ~ 1 is described in this section. The end result of this 
theory is that expressions for extinction (or total) and scattering cross-sections 
<7^ and ct s , respectively (or efficiency factors and Q s ) , are obtained as a 
function of or. These quantities are then used in the emissivity expressions as 
discussed in the third part of this section. This section shows that one of the 
limiting cases of the Mie theory in which a « 1 results in the Rayleigh theory 
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of scattering as discussed in the following part. The other limit in which a » 1 
belongs to the realm of geometrical optics and will not be analyzed here. 

The Mie theory has been used by hundreds of authors since Mie [265] 
first developed it in 1908. Perhaps the most extensive lists of references on 
the general subject of Mie theory and light-scattering calculations and experi- 
ments are given in Ihe bibliography compiled by Boudreau and Stone [3531 (494 
references) , and in the reports and articles by Oster [243] (227 references) , 
Hawksley [281] (227 references), Love and Wheasler [354] (167 references), 
and Pendorff [355] (117 references). The general derivations of the equations, 
starting from the Maxwell equations, are perhaps developed best by Stratton 
[232], Van de Hulst [246], Born and Wolfe [235], Havard [356], Svaton and 
Winer [357], Wyatt [358], Panofsky and Phillips [57], Goody [25], Morse 
and Feshbach [359], La Mer [360], Bromwich [361], Nawrocki and Papa [362], 
and Newton [363] . In this section the principal parts of the derivations as 
discussed by the above authors are given to determine the extinction and scatter- 
ing cross-sections for use in emissivity predictions. This section covers the 
following areas: (1) determination of Mie coefficients, a m and b m ; (2) deter- 
mination of extinction and scattering cross-sections, cr t and o s , and efficiency 
factors, Qt and Q s ; (3) limiting case of Rayleigh scattering; and (4) discussion 
of the Mie theory calculations. 

1. Determination of Mie Coefficients, a m and b m . In analyzing the 

problem of scattering of plane waves by homogeneous spherical particles, it is 
usually assumed that the incident radiation is linearly polarized and that the 
particles are randomly distributed and separated from each other by distances 
that are large compared to the wavelength of incident radiation. The Mie co- 
efficients, a m and b m , (and consequently the cross-sections and efficiency 
factors) may then be obtained by solving Maxwell's equations and applying 
appropriate boundary conditions. Maxwell's equations in Gaussian units are 
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and 


V x E 
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(342) 
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A third equation, the equation of conservation of charge, will be given 
as 


V • 



= 0 , 


(343) 


where p e is the charge density and J is the current density. If the time 
dependence of the steady-state fields may be described by the factor e ia;t , 
equations (341) and (342) reduce to the form 


V x H = i kn 2 E 


(344) 


and 


V x E = - ikH 


(345) 


where the quantity kn is the propagation constant in a medium of complex 
refractive index n (the quantity k represents the propagation constant in a 
vacuum) . Tne quantities k and n may be defined as 


k 
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This notation differs slightly from that used by Stull and Plass [240] in 
equation (264) . 


For a homogeneous medium, it can be shown that the electric and 
magnetic field vectors, E and H, respectively, satisfy the vector wave equation 


V 2 A + k 2 nA = 0 


(347) 


where the vector A may be either E or H. This equation may be solved 
(Havard [356]) by solving the following scalar wave equation: 
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V 2 ip + k 2 nip = 0 . 


(348) 


If $ is the solution of equation (348) , then the vectors and N^, defined by 
van de Hulst [246] as 


M^Vx (rip) and nkN^ = VxM^ , 


(349) 


will satisfy the vector wave equation (347) . These vector quantities are related 
by the equation 


nkM = V x N . (350) 

4> <P 

According to van de Hulst, if u and v are two solutions of the scalar 
wave equation (348) , and' if M u , N u , M y , and N y are the derived vector fields, 
the equations (344) and (345) may be satisfied by the equations 


H = n(-M + iN ) 
u v 


(351) 


and 


E = M + iN 
v u 


(352) 


By constructing three vector solutions L ^ M^, and N^ for 

equation (347) , it may be shown by Stratton [232] that the full components of 
and in polar coordinates may be written as 


M = 0, 
r 


nkN 

r 


9 2 ( rip) 

dr 2 


+ n 2 k 2 rip 


(353) 


M = 

0 r sin 9 dcp 


nkN 0 


1_ a 2 ( rip) 

r 9r 90 


(354) 
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and 



i_ djrijj) 

r 80 


nkN = — — 

<p rsinS drdd 


(355) 


Stratton and van de Hulst show that the magnetic and electric fields, 

H and E respectively, may be completely described by choosing solutions 
u and v to the scalar wave equation (348) , for the incident wave outside the 
spherical particle, the scattered wave outside the spherical particle, and the 
field inside the spherical particle. For the outside, incident wave, the solutions 
are 


u = aY K J (kr) , 
— , m m 
m=i 


= B V K J (kr) 
Lj m m 


m=l 


(356) 


For the outside scattered wave, die solutions are 


u = aV -a K (kr) , v = B Y -b K H (2 ^(kr) 

u ■ ill ill ill Lj m m m 


111=1 


111=1 


(357) 


For the inside wave, the solutions are 


u = A V nC K J (nkr) , v = B V nd K J (nkr) . (358) 

u m m m Lj m m m 


m=i 


m=l 


In these equations the quantities A, B, and K are 

m 


a iwt ^ ^ icut . 

A = e cos <p, B = e sin<£ 


(359) 


and 


K 

m 


2 m + 1 
m(m + 1) 


P' (cosO) 
m 


(360) 
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The quantities P' m are the associated Legendre polynomials with 1=1. 
The a m and b m are the Mie coefficients to be determined subsequently, and 
the c m and d m are related to the a m and b m . The quantities H n ^ are the 
Hankel functions related to the spherical Bessel functions J m (z) and N m (z) 
as given by Morse and Feshbach [359] as 


H 


(2) 


m 


(z) 


= J (z) - iN (z) 
m m 


(361) 


By using the Riccati-Bessel functions, 

( 2 ) 

ip (z)=zJ (z), X (z) = -zN (z) , and £, (z) = zH (z), (362) 

m m m m m m 

and equation (361) , the following relation may oe written: 


t, ( z) — ip (z) + ix (z) . ( 363 ) 

m m m 

The components of the electric field vector in equation (351) , E 0 and E^, 

1 0 

both contain the quantities v and ~ ( ru ) > which are different inside and 

outside the spherical particle. The components of the magnetic field vector in 
equation (351) , H 0 and H^, both contain the quantities nu and ^ (rv) which 


are also different inside and outside the spherical particle. All of these quan 
tities must be equal on the boundary of the spherical particle. Van de Hulst 
[246] presents boundary conditions which equate the quantity in brackets on 
either side of the boundary of the spherical particle. 

[nu]: 

ip (a) -a l (a) = nc ip (p) 
v mm mm 

(364) 

J* 1 9(ru)~ 

L n dr _ 

: ip' (a) -a t,' («) = c ip' (/3) 

m mm m r m 

(365) 

[v]: 

ip (o) -b C («) = d ip (fi) 
r m m m m r m 

(366) 
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(367) 



where the primes denote derivatives and ft - na - n 


2 7r r 



If c m is eliminated from the first pair of equations and d m from the 
second pair, the Mie coefficients a^ and bj^, the electric and magnetic 
wave coefficients, respectively, may be obtained as 


a 

m 


»m (g> »m (a> ~ a *m m <n (a) 
t m , “ ) ’ ^m ( “> 


(368) 


and 


b 

m 


* m {a) ~ ^m (/J) C (a) 
" Sn (a) ‘ ^m (/j) 


(369) 


where equation (362) would then be used. The Mie coefficients are sometimes 
expressed in other ways, such as with Neumann functions, 


Y (z) 
m 


J (z) 
m 



m 



(370) 


as written by Stull and Plass [240] as 

-l 


a 

m 


1 + i 


nY..(«) (P) - JJP> Y m+l {a) 


111 


m+1 


m 


"•>-<“> J m+1 (,i) - J m ,W 


111 


111 


(371) 


and 


b = - 


m 


i + i 


(m + 1 ) (n 2 - 1 ) Y m (q) J m (g) + P Y m <°> J m+1 Q) ~ »*« J m W Y m+ l (a) 
(m + 1) (n 2 - 1) J m (o) J m O) + «« J m+i ^ 


1 


(372) 
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or in terms of the logarithmic derivative functions as used by Infield [364] , 
Aden [365], Krascella [259], Beheshti [260], Tibbodeaux and Beheshti [12], 
and Smith [366]. These logarithmic derivative functions simplify the calcu- 
lations when the index of refraction is complex and when the Bessel and Hankel 
functions are small and complex. In this manner the two functions 77 (z) and 
4 m (z) may be defined as ' 


v <z) 




H <2) <z) 
m-1 

H ( 2 ) (z) 

m 


m 

z 


(373) 


and 




In 


z J (z) 
- ni 




m 


(374) 


Using these quantities, the Mie coefficients may be written as 


a 

m 


and 


J m ( “> 
H m <“» 


v (t> * - 

i m («) -«v» 


b 

ni 


j m < a) 

H (2) (a) 
m 


1 m (W - n, m ( tt ) 

V (W 


(375) 


(376) 


2. Determination of Scattering and Extinction Cross-Sections, <t s and a t , 

and Efficiency Factors, Q and Q . The scattering and extinction cross-sections 
s 

may be obtained by either of two methods. The first method involves the use of 
the following scattering matrix as discussed by van de Hulst [246] : 


( 




S 3 

Si 


ikr 



(377) 
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where the parallel and perpendicular components of the electric field of the 


incident wave, 


E 


and E 


°X 


, respectively, are 


E = cos <p 

r\ i i ! 


and 



- sirup 


(378) 


and the parallel and perpendicular components of the scattered wave are 


E|! = E e and E^-E^ . (37 

In the matrix, equation (377) , the quantity Sg = S 4 = 0 for spherical 
particles, and the quantities Sj and Sg depend only on Q. Substituting the 
asymptotic form (i n+1 /kr) e* lkz for H^(kr) for large distances r from 
the particles in equation (357) , for the scattered wave, the resulting field 
components may be written as 


E ,, ^ E„ - H 
II 0 <p 


— exp(-ikr+ icut) cos <p Sg(0) , 


(380) 


and 


E , = - E = H = 


(p = H 0 = " Itr exp * ~ ikr + iuJt ^ Sin ^ S i < d ) 


(381) 


where the functions S t (0) and Sg ( 6» ) are, thus, 


Si «) - T . , .. fa (cosO) +b r (cos0) (382) 

u . m(m + 1) m m m m 

m=l L J 


and 


S2<e> = £ 

m-1 


2 m + 1 

m(m + 1) l m m 


[ b 7 t (cosfl) + a t ( cos 6 ) 
m m ill m J 


(383) 
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The quantities 7r m (cos0) and r m (cos0) are defined as 


7T (COS0) 
m 


1 

sin 0 


P 1 ( cos 0) = 
m 


dP (cos0) 
m 

dcos0 


(384) 


and 


r (COS0) 
m 


j f~ 

— P 1 ( cos 9 ) = COS0 7T (COS0) 
d0 m v m 


-sin0 


d 7T (cos 0) 

m 


d cos 0 


(385) 


The total or extinction cross-section is related to the total or 
extinction efficiency factor, Q t , as a t - Q t ^r p . The efficiency factor Q t may 
be obtained from equations (382) and (383) by setting 9 = 0 so that S^(0) — 

^ (0) - S(0) for 0=0. The total efficiency factor is then found, as described 

by van de Hulst [246] , as 

OO 

Q =A- Re[S(0)] = -K Yj (2m+ 1) Re (a m + bj , (386) 

1 “ m=l 

where the expression Re (a m + b m ) denotes the real part of a m + b m . 

The scattering efficiency factor, Q s = Ug/Trip, may be written in 
terms of the intensity functions ij and i2 which are defined as 


ii = 



and i 2 



(387) 


The quantity q is proportional to the intensity whose electric oscillations 
are perpendicular to the plane of scattering (polarized light in the horizontal 
plane) . The quantity i 2 is proportional to the intensity whose electric oscilla- 
tions are in the plane of scattering (polarized light in the vertical plane) . The 
polarization, P, is then related to these intensity functions by the following 
equation: 

_ h (0) ~ *2(0) (388) 

ii ( 0 ) + i 2 ( 0 ) 
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According to Churchill et al. [42], the quantities q(0) and i 2 (0) can 
be related to i(0,<p ) , which represents the power scattered by a spherical 
particle from a linearly polarized beam of unit intensity into a unit solid angle 
in the direction ( 9,<p ) . This quantity i(e,<£) is written as 






cos 2 <p + i 2 (0) sin 2 



(389) 


For a randomly polarized beam of unit intensity, the power scattered 
into a unit solid angle may be obtained by averaging i(e,0) over cp, as 



[h(6) + i 2 (0)] 


(390) 


The scattering cross-section may now be defined according to Churchill 
et al. [42] as the total power scattered by a particle from a beam of unit 
intensity as 


7T 

a - 2n f i(0) sind d 6 
S 0 


(391) 


This definition is in agreement with that of van de Hulst [246] , who 
states that if the total energy scattered in all directions is equal to the energy 
of the incident wave falling on the area <r s , the scattering efficiency factor may 
be written in the integral form: 

a 1 ^ 

Q = — -j- = — y J [ij(0) + i 2 (0)] sine de . (392) 

s n P “ 0 


The conversion of equation (392) into the following one in terms of the 
Mie coefficients a m and b m is quite difficult and has been given by Debeye 
[367] as 



(393) 
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The absorption efficiency factor Q a or cross-section a a may thus be 
found as the difference, respectively, between Q t and Q s or and ct s . 

The second method of obtaining the total and scattering cross-sections 
or efficiency factors involves the use of the Poynting vector as discussed by 
Stratton [232], Svaton and Winer [357], Panovsky and Phillips [57], and 
Goody [25] . The energy absorbed from an incident electromagnetic wave by 
a spherical particle may be written as the divergence of the real part of the 
radial component of the complex Poynting vector S' 1 ' as 


V • Re (S*) R = Re / • dA . 

A 


(394) 


The radial components may be written as 


- T [< ff i + *1 x <If i + ff s>] R - T < E *i V E *i fy 


(395) 

E O> s H 0i ) . 


where the subscript i represents the incident wave and s represents the 
scattered wave. 

If the external medium is assumed to be non-conducting, the first 
term in parenthesis on the right is equal to zero. The second term on the 
right measures the outward flow of the scattered energy from a concentric 
spherical surface of radius r which encloses the spherical particle of radius 
rp. The third term on the right is tnus equal to the sum of the absorbed and 
scattered energy which may be written as 


W = W + W = 
t a s 


(396) 



tt 2tt ^ ^ 


r 2 sin0 d0 d<£>. 
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The energies W a and W s may be evaluated by making use of asymptotic 
forms of Bessel functions J m (z) and Hankel functions (z) and with the use 
of certain orthogonal functions as discussed by Stratton [232] . The energy 
terms are thus 



Y (2m + 1) 

m=l 



2 


2 i 

a 

m 

+ 

b 

m 

/ 


and 


w t = 7T 


IT ft) 


Re Y ( 2 m + 
m=l 


1) (a + 
m 


b ) 
m 


(397) 


(398) 


The scattering, absorption, and total (or extinction) cross-sections 
may be defined as the ratio of scattering, absorption, or total energy to the 
energy of the incident wave, Wj or 


a 

s 



W 


°a ” W. 5 
i 


W 



(399) 


Since the mean energy flow of the incident polarized electromagnetic 
wave per unit area (energy density) is 


W. = S* = -J- E 2 (—\ 2 , (400) 

l 2 o \n 2 ) 

the scattering and total cross-sections may be written as the ratio of scattered 
and total energy per second to the energy density of the incident wave, or 


<7 

S 



Re Y (2 m + 1) 
m=l 



and 


a 


t 


OO 



(2 m + 1) (a + b ) . 

m m 


(401) 


(402) 
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These expressions reduce to equations (386) and (393) if the relations 
k 2 - a / rp = 27r/\, Q s = cr s /7rrp and - a^/irr^ are used. 

3. Limiting Case of Rayleigh Scattering. It may be shown (Stratton 
[232] , Born and Wolfe [235]) that a limiting case of the Mie theory of scatter- 
ing occurs when the quantity a « 1 and the Rayleigh theory discussed previously 
may be used to compute the cross-sections and efficiency factors. In this 
manner the Mie coefficients a m and b m may be expanded in powers of a. By 
neglecting all powers above the fifth, and assuming /iq = jiq, the Mie coefficients 
become 


a i 


45 


(n u - 1) cr 


(403) 


2 i / n 2 - 1 
3 \ n 2 + 2 




JL n 4 - 1 

10 n 2 + 2 



(404) 


and 


b 2 


1 n 2 - 1 r, 

-> 0 !° 

15 2 n 2 + 3 


(405) 


If a is so small such that a 5 may be neglected compared to a z , only 
the first-order electric oscillation term needs to be considered, which is 


b i = 


2 i 
3 


n 2 - 1 3 

n 2 + 2 “ 


(406) 


If this expression is used in equations (380) and (381) , it may be 
shown, as in van de Hulst [246] or Stratton [232] , that the resulting field of 
the fundamental mode is that of an electric dipole which has a dipole moment 
equal to 


p - ex - /3E 



(407) 
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This expression was determined in the previous discussion of Rayleigh 
scattering. By using the above procedure, it may be shown (van de Hulst) that 
the first terms oi t We expansion for Q ^ may be written as 


Q 


t 


- Im 



n 2 - 1 
n 2 + 2 



a 


n 2 + 1 \ 2 / n 4 - 27 n 2 + 3s \ 
n 2 + 2 / \ 2 n 2 + 3 / 



+ . . . . 


(408) 


It will be recognized that the first imaginary term is the absorption 
efficiency factor as given by equation (289) and the first real term is the 
scattering efficiency factor as given by equation (288) . 

4. Discussion of the Mie Theory Calculatio n. Although the main pur- 
pose of this entire section is to discuss the use of the Mie theory for A1 2 0 3 
particles immersed in a rocket exhaust gas, it appears necessary to cite the 
work of investigators whose main purpose was to make Mie theory calculations 
for other types of particles besides AI 2 O 3 (other than carbon, discussed in 
part 1 of "Methods of Predicting Carbon Absorption Coefficients and Emissivities 
in the preceding section) for various values of index of refraction, particle 
radius, and wavelengths. Some of the work to be cited involved light-scattering 
experiments and calculations using Mie theory in attempts to obtain particle 
size distributions. These investigations were based on measurements of 
scattered intensity as a function of scattering angle Q and the corresponding 
Mie theory calculations. 

As can be understood from the many equations of the previous sections, 
the calculations of the cross-sections (or efficiency factors) using the Mie 
theory are exceedingly lengthy - especially if the index of refraction of the 
particles is complex (absorbing particles) and a size distribution of the 
particles is used rather than one particular size. Only in the last eight years 
or so when high-speed electronic computers became readily available, have 
the Mie theory calculations become feasible to perform. Throughout the last 
10 to 15 years, however, there have been certain groups in the United States 
that have performed most of the Mie calculations and have published some 
rather lengthy tables for particles of various sizes, indices of refraction, and 
wavelengths. Although all die publications of even these seven major groups 
are too numerous to list completely, some of their most important publications 
are mentioned. 
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This list of groups includes (1) the group at the University of Michigan, 
Gumprecht and Sliepcevich [368-370], Boll, Leacock, Clark, and Churchill 
[371], Churchill, Chu, Leacock, and Chen [48], Churchill, Chu, Evans, Tien, 
and Page [42], Chu and Churchill [372], Boll, Gumprecht, and Sliepcevich 
[373], Boll and Sliepcevich [374], Chin, Sliepcevich, and Tribus [375,376], 
Clark, Chu, and Churchill [377] , Gumprecht, Sung, Chin, and Sliepcevich [ 378] ; 
(2) the group at Rand Corporation, Deirmendjian [379-381] , Deirmendjian and 
Clasen [382], and Deirmendjian, Clasen, and Viezee [383]; (3) the group at 
AVCO, Pendorff [384-392]; (4) the group at Wayne State University, Heller 
and Pangonis [393], Heller and Pugh [394], Heller and Tabibian [395], Heller, 
Nakagaki, and Wallach [396], Heller [397], Heller and Nakagaki [398], 

Nakagaki and Heller [399,400], Pangonis and Heller [401], Stevenson, Heller, 
and Wallach [402], and Tabibian and Heller [403]; (5) the group at Clarkson 
College of Technology, Kerker [404] Kerker and Matejevic [405,406], and 
Kerker and Hampton [407]; (6) the group at the University of Oklahoma, Love 
[29], Love and Wheasler [354], and Love and Beattie [408]; and (7) the group 
at Aeronutronic , Plass [409-411] , Bauer [412] , and Bauer and Carlson [413] . 

Of these seven groups perhaps the only one actively engaged in Mie 
theory calculations for A1 2 0 3 particles is Aeronutronic, although it is known 
that Svaton and Winer [357] at Douglas performed Mie theory calculations for 
A1 2 0 3 particles several years ago (as discussed in the last part of this section) . 
However, the general results of the Mie theory calculations are the same 
whether the particles have the same index of refraction as alumina (A1 2 0 3 ) or 
some other index of refraction. That is, as the value of a increases, the 
value of Q t asymptotically approaches a value of 2.0 (total cross-section 
becomes twice as great as the geometric cross-section, 7 rrp) . This phenomenon 
may be seen from the results of Plass [409] for A1 2 0 3 , reproduced here in 
Figure 36, which also shows the Rayleigh limit. Also, for all values of index 
of refraction and values of a in the Mie region, the intensity functions i t (0) 
and i 2 (0) are comparatively high in the forward scattering region (0 — 0) and 
tend to decrease significantly with increasing 0 until they reach a minimum in 
the neighborhood of 0 = 100° to 110°. The intensity functions then tend to rise 
gradually toward the back scattering region (0 — 180°) , but still have a value 
at 0 = 180° usually an order of magnitude lower than in the forward scattering- 
region. 


To predict cross-sections of alumina particles adequately using the Mie 
theory equations, the index of refraction of alumina (both the real part, n 1} 
and the complex part, n 2 ) must be known as a function of both temperature and 
wavelength. Values of index of refraction at room temperature for alumina in 
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FIGURE 36. TOTAL EFFICIENCY FACTOR FOR ALUMINA 
AS A FUNCTION OF a FOR X - 2 m 

the infrared wavelength region up to X = 5. 5 m are given by McCarthy, Wolfe, 
and Ballard [414] , and Malitson [415,416]. Absorption coefficients and trans- 
missivities (from which n may be determined) for alumina at both room tem- 
peratures and temperatures up to as high as 1000° C have been reported by 
Mergerian [417], McAlister [418], Monroe [419], Openheim and Even [420], 
and Lee and Kingery [421] . Gryvnak and Burch [422,423] at Aeronutronic 
measured the absorption coefficient and index of refraction (n^ and n 2 ) of single 
crystals of A1 2 0 3 (sapphire) which were heated in a furnace at temperatures up 
to 1700° C and over a wavelength range from 0. 546 n to 6. 0 p- Measurements 
of index of refraction and absorption coefficient of molten A1 2 0 3 (at a tempera- 
ture of 2020° C) were also obtained by Gryvnak and Burch by immersing the 
alumina samples into a hydrogen-oxygen flame. Adams and Colucci [424] 
report measurements of the complex part of the index of refraction, n 2 , which 
is equal to (X/47r) k a , of alumina at temperatures from about 2100° C to 2600° C 
based on their own and Carlson's [345] emissivity measurements in hydrogen- 
oxygen flames. The real part of the index of refraction, n t , may be calculated 
from these emissivity measurements by combining equations (64-66) , (267) , 
and the above relation for n 2 . 
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Plass [410] used the values of index of refraction as determined by 
Gryvnak and Burch [422,423] as a function of temperature and wavelength and 
found, by performing Mie theory calculations, that the absorption cross-section 
increased by a factor of about 40 as the temperature increased from 1200® C to 
2020° C. These Mie theory calculations (using Gryvnak' s and Burch's data) 
are believed to be the most valid for A1 2 0 3 at the present time and were used in 
the radiation calculations discussed in the last part of this section. In the next 
section, various experimental data and theoretical predictions for A1 2 0 3 particle 
sizes and distributions are given, as these data are also a very important input 
to the Mie theory for cross-section and, consequently, emissivity predictions. 


Determination of AI2O3 Particle Sizes and Distributions 

The size of A1 2 0 3 particles found in exhausts of aluminized solid propellant 
motors depends on the combustion process of the aluminum in the propellant, 
the A1 2 0 3 condensation process (liquid particles changing to solid particles) , 
and possibly on chamber pressure, chamber residence time, and throat size of 
the engine. Aluminum oxide particle sizes have usually been determined by 
one of three experimental methods: (1) firing small engines into tanks and 
collecting the residue on the walls, from which it is then scraped off and analyzed 
under a microscope or with a Micromerograph; (2) sampling the particle plume 
with probes, microscope slides, and other types of collecting devices and con- 
sequently analyzing the deposits; and (3) using optical (light-scattering) tech- 
niques. Carlson [5] gives an excellent review of sizes and distributions for 
both A1 2 0 3 and BeO particles obtained by these three general methods. In the 
present section Carlson's review is summarized and additional work performed 
in this area which was not mentioned by him is reviewed. Since the combustion 
of aluminum powder to A1 2 0 3 particles is an important process in particle size 
determination, it is briefly treated first in this section. The experimental 
determination of A1 2 0 3 particle sizes and distributions and the theoretical deter- 
mination of A1 2 0 3 particle sizes and distributions are then discussed. 

!• Combustion of Aluminum Powder to Al 2 O a Particles. The size and 
distribution of the A1 2 0 3 particles found in the exhaust depend upon the com- 
bustion process of the aluminum powder (which usually has a size range of 10 
to 25 p) in the propellant. Many investigators have analyzed this phenomenon 
of metal (such as aluminum) combustion in solid propellant motors. However, 
there are some disagreements as to whether the metal burns by a surface or 
by a vapor phase combustion process. These two combustion processes are 
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distinguished, according to Glassman [425], by the fact that, if the boiling 
temperature of the A1 2 0 3 is higher than that of the aluminum powder, the com- 
bustion process takes place in the vapor phase; however, if the aluminum boil- 
ing temperature is greater than that of the A1 2 0 3 , a surface combustion process 
occurs. These two processes are strictly for steady-state combustion, because 
the ignition phenomena could be quite different. 

Some of the many investigators who have analyzed surface and vapor 
phase combustion processes in solid propellant motors include Glassman [425, 
426], Brzustowski and Glassman [427] , Fassell et al. [428-430], Gordon 
[431-433], Wood [434,435], Avery [436], Henderson and Bowen [437], 
Friedman and Macek [438], Markstein [439], Davis [440], and Kurtovich and 
Pinson [441,442] . ' In the particular case of aluminum, experimental data 
such as obtained by Wood [435], Watermeier, Aungust, and Pfaff [443], and 
McCarty [444] tend to favor the vapor phase combustion process involving the 
formation of A1 2 0 3 particles. 

Fassell et al. [428] state that the aluminum combustion process takes 
place in four distinct phases: (1) pre-ignition, (2) ignition, (3) bubble expan- 
sion, and (4) quasi-steady-state combustion. According to Avery [436], the 
oxidizer (such as NH 4 C10 4 ) and the binder are first vaporized, heated to com- 
bustion temperature, and then carry the aluminum particles with them into the 
flame front. The molten aluminum is heated, according to Wood [434] , first 
to a dull red color, which then changes to a yellow-orange color. The actual 
combustion takes place by the establishment of an incandescent reaction zone 
around the particle. Davis [440] states that the diameter of the aluminum 
flame is greater than the diameter of the aluminum particle, thus indicating that 
the reaction zone is not situated on the particle surface, but is actually some 
distance from the particle - similar to the halo of a diffusion flame found 
around a burning liquid droplet. It is known that this zone expands during the 
ignition process and slowly contracts during the remainder of the combustion 
process. The ignition of large particles is delayed longer since the larger 
particles conduct more heat from the surface than the smaller particles. 

As the aluminum particle is heated to ignition, a shell of oxide forms 
on the particle. Because of the intense radiant heating occurring at this time, 
the aluminum is quickly vaporized and, according to Avery [436] , as the oxygen 
immediately surrounding the burning particle is rapidly depleted, the vaporized 
aluminum attempts to obtain oxygen away from its surface. Since this process 


*Also additional information was obtained from D. D. Kirtovich and G. T. 
Pinson in an unpublished Boeing report of January 1960, entitled "Radiation 
Characteristics of Rocket Engine Exhausts. " 
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has proceeded so rapidly, the oxide shell is not able to increase in thickness and 
is consequently vaporized by the high radiant heating. This entire process takes 
place in the vapor phase, with suboxides such as AlO, A1 2 0 2 , and A1 2 0 possibly 
also occurring, according to Gibbons and Siegel [445] . The liquid A1 2 0 3 is then 
formed during the expansion process in the nozzle when it is condensed from 
the vapor phase. 

The surface combustion process occurs as described by Brown and 
McCarty [446] and Kurtovich and Pinson [441,442] . In this process, the 
aluminum metal coalesces to large droplets on the propellant surface which are 
then expelled from the surface and have an oxide coating formed around them. 
After melting, the entire particle assumes a spherical shape because of surface 
tension and vapor pressure. Continued combustion vaporizes the aluminum 
inside the shell as the higher pressure inside the shell causes the aluminum to 
be forced out in a stream of liquid vapor. The A1 2 0 3 particles are thus formed 
by rapid oxidation of the molten aluminum jet coming out of the shell. It is 
possible, according to Kurtovich and Pinson, that the oxide shell, which is 
made rigid by the internal pressure caused by the boiling aluminum, will 
occasionally burst, sending shell fragments out into the exhaust plume. To 
understand this surface mechanism of combustion satisfactorily, it is necessary , 
according to Markstein [439], to analyze the absorption of oxygen on the sur- 
face of the molten oxide and the transport of oxygen or metal through the oxide 
layer. Also, it is necessary to consider the reactions which occur at the oxide- 
gas and at the oxide -aluminum interface as well as those within the oxide layer. 

Because most of the A1 2 0 3 particles recovered from the exhausts of 
aluminized motors are small solid spheres of a-alumina, rather than hollow 
oxide spheres, it is believed that the vapor phase combustion process of alumi- 
num is the type of combustion process actually occurring in this type of motor. 
According to Avery [436] , when hollow spheres are found in exhausts, such as 
mentioned by Kurtovich and Pinson, poor combustion occurred (such as at the 
end of a firing) , the radiation heating in the chamber was lower, and the 
particles were formed by a surface mechanism process. Also, Brown and 
McCarty [446] and Avery [436] state that if the surface mechanism process 
were involved, rather than the vapor phase, the size of the A1 2 0 3 particles 
found would be greater than the size of the aluminum particles (10 p < d p < 25 p) 
placed in the propellant. Tnis is not the case, as will be seen in the next 
section. 


❖ Also additional information was obtained from D. D. Kirtovich and G. T. 
Pinson in an unpublished Boeing report of January 1960, entitled "Radiation 
Characteristics of Rocket Engine Exhausts. " 
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2. Experimental Determination of Al ? Q a Particle Sizes and Distribu- 
tions. The sizes of AI2O3 particles collected from rocket exhausts t>y the major 
ity of investigators are found to have an average particle radius, r p , of between 
1 and 3 ju. Some investigators, such as Sehgal [447] and Allport et al. [448] , 
by firing small rockets into tanks, found a definite size dependence on motor 
chamber pressure. Other investigators, such as Dobbins [449,450] and Dobbins 
and Jizmagian [451,452] using optical techniques, and Brown and McCarty [446] 
by passing microscope slides through the plume, did not find such a chamber 
pressure dependence. 

Carlson [5] presents a figure, reproduced here as Figure 37, which 
shows this mean size (particle diameter) obtained by several investigators as 
a function of chamber pressure. In most of these experimental tests, the 
motors used were of a relatively small size having a throat diameter of one 
inch or less. However, as can be seen from Figure 37, Nack's, Delaney's, 
and Bartlett's data* were for one rather large motor (d t = 7. 5 in. ) and one 
very large (120-in. solid) motor (d t = 37. 8 in.) , and the particle sizes found 
were several times higher than those of the smaller motors. Cheung's and 
Cohen's [453] set of data (for the longer residence time) falls right on Sehgal 's 
curve; however, the short-residence time data shows a much smaller mean 
particle size. Allport's data, curve 3, which were collected from motors 
fired into a tank, agree closely with Sehgal' s data at low values of chamber 
pressure and with Brown and McCarty's data at the higher values of chamber 
pressure. Dobbin's data obtained by light-scattering techniques using Mie 
theory (discussed in the next section) showed the lowest particle size of all 
(dp about 0. 4-0. 6 p) . 

The Air Force Rocket Propulsion Laboratory at Edwards Air Force 
Base, California, has made considerable measurements recently (some of 
which are shown in Fig. 37) of particle size from exhausts of 120-, 156-, and 
260 inch diameter solid-fueled motors. Smith [459] has found that the mean 
particle size can be represented best as a function of throat radius for a large 
range of motor sizes. Figure 38, reproduced from Smith, shows this curve of 
d43 as a function of r fc showing a correlation not seen previously (such as in 
Fig. 37) . In Figure 38 the 260-inch motor had a throat radius of 35 inches; 
the 156-inch motors had throat radii of 30 and 17. 5 inches; and the 120-inch 
motors had throat radii of 18. 9 and 12 inches. The particle collection technique 


* This information is contained in two 1966 Aerospace Corporation unpublished 
reports: "Two Phase Flow Phenomena in Rocket Motors" by T. H. Nack and 
L. J. Delaney, and "Particle Size in Two Phase Nozzle Flow" by R. W. 
Bartlett and L. J. Delaney. Additional information was received through 
private communication with Mr. Delaney and the author. 
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for the vertically-fired 260- and 156-inch motors was rather unique in that an 
RB-57 airplane flew into the cloud of the exhaust about 15 to 20 seconds after 
the end of the firing and collected the particles in a wing-tip pod which had a 
10-inch opening for the particles. 

Since many of the investigators analyzing particle sizes use various 
definitions for an average or a mean size, it is necessary to list the various 
types of sizes (diameters) mentioned in the literature. Carlson [5] presents 
such a list as follows: 
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and 4. Volume-to-surface diameter, d 32 : 
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Some of these values of particle diameter are used in the theoretical 
determination of A1 2 0 3 sizes, as discussed in the next section. 

Other investigators not mentioned by Carlson [5] have obtained A1 2 0 3 
particle sizes similar to those shown in Figure 37. Preckel, Jacobs, and 
Gibson [455] at Allegany Ballistics Laboratory placed microscope slides in 
the exhaust, mounted at distances of 15 to 100 feet behind the motor exit. 
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The mean diameter of the particles collected on the slides in all cases was 
between 2 and 3 p. Optical techniques were also used in which film density 
measurements were made with a Densicron on the photograph, and transmit- 
tances of visible light through the plume at a single wavelength for black-and- 
white photographs were obtained and at three wavelengths for color photographs. 
Using the Beer-Lambert equation (discussed in the next section) , the optical 
depth was found to yield a mean particle diameter of about one micron. 

Kurtovich and Pinson [441,442]'" at Boeing obtained AI2O3 particle 
samples on glass slides from both 1/20^ and full-scale Minuteman firings 
at a chamber pressure of about 760 psia. It was found that two types of parti- 
cles existed on the slides - hollow micronic particles and solid, submicronic 
particles. The hollow micronic particles showed a distribution around dp 4 p, 
while the number of solid submicronic particles appeared to increase as the 
particle size decreased. Measurements at Rohm and Haas* * showed a mean 
A1 2 O s particle size of about 2 to 4 p; however, the results were believed to be 
inconclusive because of the amount of agglomeration and coalescence that 
resulted. Based upon theoretical specific impulse calculations using various 
propellants fired at Rhom and Haas, it is believed that the mean particle 
diameter should have been closer to one micron. Dimmock and Courtney [337] 
at Thiokol, who collected A1 2 0 3 particles in rocket exhausts in which the pro- 
pellant had 8 percent aluminum and where the motor was fired at 750 psia, 
found the mean particle size to be between 0. 1 and 1. 0 p. Povenelli and 
Rosenstein [456] obtained a mean particle diameter of about 1.6 p in aluminized 
(9 percent) propellant strand burning tests. Burns [457] inserted a probe in 
the exhaust of an motor whose propellant had 19 percent aluminum and was 
fired at 630 psia. He obtained the very large mean diameter of 150 p; however, 
the particles collected were of a large porous nature, and it was believed that 
the probe-collecting device was biased in favor of the large particles, because 
many of the small particles followed the gas- stream lines and did not go into 
the probe. 


*Also additional information was obtained from D. D. Kurotvich and G. T. 
Pinson in an unpublished Boeing report of January 1960, entitled "Radiation 
Characteristics of Rocket Engine Exhausts. " 

**W. A. Wood, private communication to the author, August 1966. 
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Measurements of A1 2 0 3 particle size and distribution were obtained at 
Aeronutronic [344,345] from the exhaust of a small hydrogen-oxygen engine 
which was fired with a slurry of AI 2 Q 3 particles dispersed in water. The A1 2 0 3 
particles coming out of the exhaust were collected during the firings by means 
of a May-Welchman cascade particle analyzer and a shutter probe which was 
placed about 15 feet from the motor exit. It was found that the distribution of 
the collected A1 2 0 3 particles (which had a maximum diameter of about 10 p) 
agreed quite well with the A1 2 0 3 particle distribution (analyzed with a Microm- 
erograph) put into the slurry before firing. 

At the time of this writing an investigation is underway by Boeing " 
performed in an altitude cell at MSFC Test Laboratory to determine A1 2 0 3 
particle sizes in small (dt = 3/4 in. ) solid propellant rocket exhausts. The 
chamber pressure is rather high (approximately 1500 psia) ; the expansion 
ratio of the motors is 8:1; and the amount of aluminum is 14 percent and 20 
percent. A particle collector which is pressurized by argon gas to slow the 
particles down (described by McGregor [458]) is mounted at several axial 
distances varying from 27 to 43 inches from the exit and also at several radial 
distances from the axis. The particles impinge on the inside walls of the collec- 
tor (after passing through an orifice) where they are flushed out with a solvent 
and analyzed with an optical microscope. It has been determined that the 
particles are nearly spherical, a result which agrees with most of the investi- 
gations mentioned above. To date (three firings) the average particle size 
(diameter) appears to be about 5 p; however, it is believed that the smallest 
particles are nof being resolved by the optical microscope, so that Micromero- 
graph techniques and possibly electron microscopes will be used in the future 
to obtain more accurate measurements. 

3 . Theoretical Determination of AI 9 O 3 Particle Sizes and Distributions. 
Various theoretical predictions based upon some of the experimental data men- 
tioned above have been suggested for A1 2 0 3 sizes and distributions. Many of 
these have been used in calculations of radiant heating from aluminized rocket 
exhausts. Perhaps the simplest method used has been to assume that all of the 
particles are the same size, such as rp = 2 p as used by Fontenot [459] and 
Bender and Mullin [460,461] . A size distribution method based upon the 
analysis of Kliegel [462] was used by Morizumi and Carpenter [463] , Gulrajani 
[464], Rochelle [465], and Hunt [466] for A1 2 0 3 particle radiation calculations. 
This was based on five sizes of particles, rp = 0. 79 p, 1. 28 p, 1. 76 p, 2. 44 p, 
and 3 . 95 p, each representing 20 percent of the total A1 2 0 3 particles in the 
exhaust. This size distribution was based upon Kliegel's [462] fitting the 
following logarithmic normal particle size distribution based upon Brown and 
McCarty's [446] data and some unpublished data obtained at TRW Systems: 


*M. Baker, private communication to the author, September 1966. 
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where d p = 3. 5 + 1. 0 p, a g = 1. 9 + 1. 0 p, m p is the particle mass flow, and 
mg is the gas mass flow. 

Crowe and Willoughby [467] performed a theoretical investigation of 
A1 2 0 3 particle growth in a solid propellant rocket nozzle by particle collision 
and coalescence in the nozzle and shock waves of the exhaust plume. This 
was done in an attempt to explain certain UTC experimental data [448,468] 
shown as curve 3 in Figure 37 which shows that particle size increases with 
chamber pressure up to a certain point, beyond which there is no increase in 
particle size. Besides chamber pressure effects on particle growth, Crowe 
and Willoughby [467] looked at effects of nozzle scale, aluminum content, and 
initial particle size effects on particle growth. The particular model which 
was used to predict the rate of particle collision per unit time is 
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where 2 is the number of particles of radius r p ^ which collide with parti- 
cles of radius r^ per unit time, v pi and Vp^ are the particle velocities, and 

Np^ is the number of particles of radius ^ per unit volume. 

Nack and Delaney discuss a simplified aerodynamic model to predict 
average particle sizes as a function of propellant composition, nozzle size and 
contour, and motor operating conditions. They also discuss a mathematical 
model which describes the gas-particle flow through a nozzle, taking into 
account both agglomeration and break-up of particles. The four following 
first order partial differential equations are solved in their analysis: (1) 
total momentum, (2) particle motion, (3) particle heat transfer, and (4) 
particle agglomeration. Theoretical calculations for four motors (d^ = 2.5 in. , 
7. 5 in. , 9. 1 in. , and 37. 8 in. ) show agreement within 10 percent of experi- 
mental data. 


Fein [469] presented a theoretical model for predicting A1 2 0 3 particle 
size distribution in which he found the normalized distribution frequency f(V p ) 
to be 

*This information is derived from a i960 Aerospace Corporation unpublished 
report by T. H. Nack and L. J. Delaney, entitled "Two Phase Flow Phenomena 
in Rocket Motors. " 
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where Vp is the volume of a single particle and Vp is the number average 
particle volume. In this analysis Fein assumed the rate of growth of a particle 
dVp/dt to be given by the relation 
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where p, MW, C, and C e are, respectively, the density of the oxide particle 
in the condensed phase, molecular weight of the oxide particle, concentration 
of the oxide particle in the gas phase, and equilibrium concentration of the 
oxide particle in the gas phase at the chamber temperature and pressure of the 
oxide particle. Also, 7 is the rate of growth constant, and Ap is the surface 
area of the particle, assuming that the particle is spherical in shape. 


McGregor [458] assumed an A1 2 0 3 size distribution based upon the 
following normalized error function for the i^ 1 particle (fraction of particles 
of the ith size) : 
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where the mean particle radius 
pressure, as based on Sehgal's 
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[447] chamber-pressure -dependent data as 
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where the chamber pressure P c is in psia and rp is measured in microns. 


The normalizing constant K is given by 
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Dobbins [449-452; 470-472] describes optical measurements using light 
scattering (Mie theory) calculations to obtain the particle size distribution 
(curve 5 in Fig. 37) . The optical transmission of small rocket plumes was 
measured at two wavelengths, 0. 365 and 1. 01 p, two spectral regions in which 
there was no H 2 0 emission and absorption. For a polydispersion of spherical 
particles, the transmission expression may be written as the Beer-Lambert 
equation (9) as 
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where f(dp) is a particle size distribution function defined such that 
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P(dp^ < dp < dp ) is the relative probability of occurrence of particles larger 
than dp and smaller than dp . A mean extinction efficiency factor, Q^, may 


be defined as 
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and, if the particle volume concentration c y may be defined in terms of the 
particle number concentration N p as 


c 
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and if the volume-to-surface mean diameter d 32 is defined as in equation (412) , 
the transmission law, equation (420) , may be written as 
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From this expression, using Mie theory and measuring the ratio I/1 G 
experimentally at the two wavelengths 0.365 and 1. 01 p, the volume-to-surface 
mean diameter d 32 was calculated and plotted in Figure 37. It is known that 
other investigators such as those at United Technology Center' 1 ' have also made 
or are planning to make light-scattering measurements and Mie theory calcu- 
lations to determine A1 2 0 3 particle sizes from rocket exhausts. 


Holland and Draper [473] present the following A1 2 0 3 particle size 
distribution for a polydisperse cloud of particles: 
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where d n is the median volume diameter of the particle, M is the mass of 

f 0 u 

the median volume diameter particle, p is the mass of particles per unit 
volume, and s and b are parameters of the size distribution obtained from 
standard laboratory measurements of particle size. This size distribution 
function was used by Holland and Draper to investigate (theoretically and ex- 
perimentally) the light-scattering effects of talc, carbon black, silicon dioxide, 
as well as A1 2 Q 3 , for the LID effect as discussed by Adams and Holland [474] . 


* W. Lai, private communication with the author, 1965. 
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This LID effect involves the interaction of plume radiation with particles 
brought over the solid propellant missile from the rocket exhaust by plume- 
induced flow separation. 

Peterson et al. [475] fitted the following equation for size distribution, 
n(dp) [number of particles whose diameter lies between dp and dp + d(dp) ] 
based upon Brownian motion theory: 
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This equation, which is similar to Bauer and Carlson's equation discussed in 
the next paragraph, showed very good agreement with the experimental dis- 
tribution (see curve 6 in Fig. 37 for mass average diameter) . 


Bauer and Carlson [413] present the following skew-symmetric equation 
for A1 2 0 3 particle size distribution (probability of finding a particle of radius 
r p to rp + drp) in a rocket exhaust: 
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This type of expression was originally derived for size distribution in 
clouds of water droplets as reported by Bauer [412] , Bartky [476] , and 
Deirmendjian [380,381] , but has been shown to fit A1 2 0 3 particle samples 
collected from rocket exhausts. Based on this distribution function, equation 
(427) , ihe average particle area for a distribution of sizes may be written, 
according to Carlson et al. [477] , as 


oo 



and the average particle volume for a distribution of sizes may be written as 
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Bauer and Carlson [413] curve fitted equation (427) to the data of Brown 
and McCarty [446] and Sehgal [447], obtaining values of a and b, using Mie 
theory, and computed the corresponding scattering and absorption cross-sections 
for several A1 2 0 3 particle distributions. Carlson [477] , using equation (427) , 
calculated the size distribution function for the S-II ullage rocket using a = 1. 895 
and b - 1. 269 as shown in Figure 39. In the analysis of particle radiation heat- 
ing from Saturn ullage and retro motors (See subsection following, entitled 
"Radiation Calculations for Aluminized Solid Propellant Rocket Exhausts. ") 
equation (427) was used to compute the absorption and scattering cross-sections 
and emissivities (discussed in the next section) . 



FIGURE 39. PARTICLE SIZE DISTRIBUTION FOR SATURN 
S-II ULLAGE ROCKET 
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Determination of AI2O3 Particle Cloud Emissivity 


An accurate prediction of the emissivity of the AI 2 O 3 particle cloud in an 
aluminized rocket exhaust depends upon the absorption and scattering cross- 
sections, which depend on the index of refraction of the A 1 2 C >3 particles as dis- 
cussed in the section on Mie theory of scattering, and upon the particular size 
distribution of particles as discussed in the following section. However, this 
cloud emissivity is also a function of optical thickness, t, of the plume and the 
type of particle scattering involved (isotropic or anisotropic, single or multiple) . 
In this section emissivity expressions for optically thick and thin particle clouds 
with various types of scattering are discussed. Three general methods for 
predicting A1 2 0 3 particle cloud emissivity are given: ( 1 ) inverse wavelength 
method, (2) neutron-scattering analogy method, and (3) one-dimensional beam 
approximation method. These methods are described in order of increasing 
complexity (and accuracy) and decreasing value of radiant heating predicted. 

1. Inverse Wavelength Method. This method, which was developed by 
Fontenot [459] at Boeing as a simple approach to the radiation heating problem, 
provides a conservative estimate of particle cloud emissivity. A form of 
equation (298) was used to obtain the particle cloud spectral emissivity as 



where the absorption coefficient k is obtained from Schack [277] as 0. 57/A. 
and V represents the volume fraction of A1 2 0 3 particles in the cloud. The 
total emissivity is obtained as by Sato and Matsumoto [279] from equation (301) 
with the exception that the quantity x is equation (301) may be expressed as 


X = 0.57 Vi T /C 2 , (431) 

Fx 

where Tp^ represents the particle temperature at a distance x from the nozzle 

exit. The quantity V may be expressed in terms of M, the mass of A1 2 0 3 
particles per unit volume, as 
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where pp represents the density of an A1 2 C>3 particle and y represents the 
quantity r - r e where r is the radial distance from the axis and r e is the 
nozzle exit radius. The quantity M can be approximated by the expression 
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where Wp represents the weight fraction of particles in the plume, F is the 
thrust of the motor, I S p is the specific impulse, and A e is the nozzle exit 
area. 


If equations (430), (301) , and (431-433) are combined, the following 
expression for the total emissivity of an AI2O3 particle cloud at a distance x 
from the nozzle exit may be written: 
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Although this prediction for particle cloud emissivity is much easier to 
make than that of the two methods discussed below, it is not believed to be as 
accurate because only one particle size (rp = 2 p) is assumed. Also, the 
effects of particle scattering are neglected, and it is believed that this inverse 
wavelength dependence for spectral emissivity is more valid for carbon particles 
(discussed in the previous section, "Radiation from Carbon Particles") than for 
A1 2 0 3 particles. 

2. Neutron-Scattering Analogy Method. Morizumi and Carpenter [463] 
and Morizumi [478] present a method of predicting the apparent surface 
(hemispherical) emissivity from an A1 2 C>3 particle cloud whose particles are 
radiating isotropically, based upon a neutron-scattering analogy. Subsequently, 
Gulrajani [464] and Hunt [466] produced computer programs based on Morizumi's 
and Carpenter's [463] equations. The particular surface analyzed by these 
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investigators was the surface generated by the limiting particle trajectory of 
the smallest particle (rp = 0.76 p) considered. Five particle sizes were con- 
sidered, corresponding to the size distribution assumed by Kliegel [462] dis- 
cussed in the previous subsection. The cloud emissivity analysis was based 
upon the method of successive generations (multiple scattering) of neutrons as 
discussed by Stuart [479], Stuart and Woodruff [480] and Anthony [481]. In 
this manner, the neutron blackness, /3(t) , which is defined as the probability 
that a neutron incident on a body will be absorbed within it, may be expressed 
as 



(436) 


If photons in a particle cloud are substituted for the neutrons in the 
analysis given [479-481] , the following expression may be written for the 
apparent emissivity of a slab or cylinder composed of A1 2 0 3 particles: 
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where P^ q represents the probability that the radiant heating originating 
isotropically from an A1 2 0 3 particle is reflected by another A1 2 0 3 particle 
within the slab or cylinder after having been reflected (i - 1) times. Accord- 
ing to Morizumi and Carpenter [463], the quantities and g s /ui are 

equivalent to the emissivity and reflectivity, €p and 1 - €p, respectively, of 
an A1 2 0 3 particle (assuming the particle is opaque) . 

Equations for the first, second, and third collision probabilities, P 0) k» 
Pj k, and P 2 k, respectively, have been derived by Morizumi and Carpenter 
[463] . The first collision probability P 0 s k which represents the probability 
that once-scattered photons escape through the surfaces z = 0 and z = l is 
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where E, 

k+1 


is an integro-exponential function which may be written as 
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k+1 p t 
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■(N p (7 t 2)/M ^ 


where p - cos0 and k is a cosine distribution exponent equalling one for iso- 
tropic scattering. The second and third probability functions, Pj ^ and P 2 
respectively, are very complex functions of the optical thickness, which was 
summed over five particle sizes (mentioned in the previous subsection) as 


T 




(439) 


where Np^ is the particle concentration for the n 1 "* 1 size particle obtained from 
a two-phase plume program and at is the total cross-section for the n^ 1 size 
particle. Morizumi and Carpenter [463] used at -- Q^7rrp^ - 27rrp^, using 
the asymptotic value of 2 (see Figure 36) for Qt for large values of a. 


It may be mentioned that an approximation to the particle concentration, 
Np , may be written without going to a two phase plume program. This approx- 
imation, which was suggested by French [482], is 
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where 0^ represents the limiting streamline of the n^ 1 size particle (given 
11 1 

as a function of in French's work) ; is the gas velocity; z is the axial 
distance downstream of the throat; K - v p/ v g (where Vp i s the particle 
velocity) ; in is the particle mass flow; anu Pp is the density of an aluminum 
particle. The quantity 1 - K (velocity lag) was plotted in his work as a function 
of rp(r^) - ^ where r^ is the throat radius. 


To obtain the value of particle emissivity, which Morizumi and 

Carpenter [463] designated as 0 ^/ 0 ^ ~ these investigators extrap- 

olated data of bulk alumina emissivity measurements and plotted their own 


175 



emissivity data obtained from rocket exhausts. The apparent emissivity, e a , 
of the A1 2 0 3 particle plume was calculated from experimental measurements 
obtained during sub-scale Minuteman engines fired at AEDC, described by 
Byrd [483] as 


€ 

a 



(441) 


where q, the heat transfer rate (per area) , was measured with radiometers 
(narrow-view angle) ; radiation calorimeters (wide-view angle) , and spectrom- 
eters (narrow-view angle) . The effective plume temperature, T, was calculated 
as will be discussed in the section immediately following, and the optical thick- 
ness, t, was determined from equation (439) . The particle emissivity €p 
was then obtained from the cylinder curves presented by Morizumi and Carpenter 
[463] in which ? a was plotted as a function of r for various values of €p. It 
was found that a value of €p of about 0.25 was valid for both the experimental 
Minuteman data and also the extrapolated bulk A1 2 0 3 emissivity data. Morizumi 
and Carpenter's apparent emissivity curves, which are presented here as 
Figure 40, can be. represented fairly accurately according to French [464] by 
the following equation: 



where C = 2/3 for a cylinder and 2 for a slab. 

It is mentioned by Carlson [477] that the ratio Q a /Qt should not be used 
for the A1 2 0 3 particle emissivity; instead, merely the absorption efficiency 
factor Q a would be the true particle emissivity as stated by van de Hulst [246] 
on page 452 of his book. If the ratio Q a /Q t were used for a blackbody (e = 1) 
in which the absorption cross-section, o a , would equal the geometric cross- 
section, 7rrp(Q a = 1) , the quantity Q a /Qt would be less than one since Q s 
(and hence Q t ) would be greater than one. Carlson [477] also believes that 
Morizumi and Carpenter's [463] experimental determination of T a is too high 
since it is based upon the Planck blackbody distribution. In the case of solid 
propellant exhausts, a number of sources (particles) will exist at different 
temperatures, and thus a Planckian distribution from such a multi-temperature 
source will not exist. Also, while the effects of one-dimensional multiple 
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Optical Thickness, r = <r % N p X 


FIGURE 40. APPARENT EMISSIVITY FOR A CLOUD OF PARTICLES 
IN THE FORM OF A SLAB OR A CYLINDER 

scattering are included in Morizumi's and Carpenter's analysis, the effects of 
anisotropic scattering (important at large values of r) are not. 

3. One-Dimensional Beam Approximation Method. This method is 
currently being used at Aeronutronic under MSFC contract to analyze particle 
radiation from the exhausts of Saturn ullage and retro motors and from exhausts 
of 120-inch strap-on motors for the Saturn improvement study. This analysis 
involves the study of optically thick plumes (where multiple and anisotropic 
scattering is important) as well as optically thin plumes in which the radiation 
from individual particles is summed over an incremental volume. It is gener- 
ally felt that the plumes of 'die Saturn ullage and retro motors are optically 
intermediate (not thin to scattering, but thin to emission and absorption) , and 
that the 120-inch plumes will possibly be optically thick, especially near the 
exit plane. Carlson [477] specifies the criteria for optical thickness by stating 
that, if 
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the plume is optically thick for emission (anisotropic scattering centers) . If 
the quantity 


i (s) /d OT zr » 1, (444) 


then the plume is optically thin for scattering ( isotropic scattering centers) . 

In these equations the quantity i ( s ) represents the scattering mean path, A p 
is the projected area of a mass-mean particle size, and d p is the particle 
plume diameter. 

In general, for an optically thin particle cloud, the cloud emissivity is 
independent of scattering angle 0, the cloud emissivity for isotropic scattering 
is equal to the cloud emissivity for anisotropic scattering, and the hemispherical 
emissivity ejj is larger than the normal emissivity For an optically thick 

particle cloud, the cloud emissivity is a function of scattering angle 0, the 
cloud emissivity for an anisotropic particle cloud is greater than that for an 
isotropic particle cloud, and the hemispherical emissivity tjj is less than the 
normal emissivity (except for the one -dimensional case) . The above facts 
may be seen in Figure 41 reproduced from Bartky and Bauer [34] , which is 
discussed later in this section. 

The particle emissivity e p , which is equal to Q a , for A1 2 C>3 particles 
was determined experimentally as reported by Carlson and du Puis [345], 
Carlson [344] , and Burch and Gryvnak [422,423] , all from Aeronutronic, and 
was found to be very much lower (especially in the solid phase) than the 0.25 
value used by Morizumi and Carpenter [463] . The Aeronutronic investigators 
found that a large change in order of magnitude of e p occurred when the A1 2 0 3 
sample reached its melting point, as can be seen from Figure 42, which is 
reproduced from Carlson [344] . Additional measurements by Adams and 
Colucci [424] have agreed closely with this set of data above the melting point. 
Carlson [351] suggests several reasons why the A1 2 0 3 particle emissivity 
changes by such a large magnitude upon melting. These reasons include (1) 
the effect of a polycrystalline structure of the A1 2 0 3 , (2) chemical reactions 
occurring at the surface of the hot particles, (3) the "searchlight effect," and 
(4) the effects of impurities on the optical properties of the particles. 

To determine the Al 2 O s particle cloud emissivity (hemispherical 
or normal e^) by the one-dimensional beam approximation, the equation of 
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3D ° 3 Dimensional a * Anisotropic 
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I * Isotropic H * Hemispherical 



FIGURE 41. NORMAL AND HEMISPHERICAL EMISSIVITY OF PLANAR 
PARTICLE CLOUDS OF FINITE THICKNESS 


transfer (72) for the plane -parallel case must be solved. Equation (72) may 
be written in the following form, which is consistent with Aeronutronic's 
notation, as 


dl 




_A ( z ,M,</>) 


dz 


N p( CT a + 0 s* 


a s 2 ” 1 


+ N p4 7 I \<z. < Jp'd*' + N p a a B x (T) , (445) 


if the following substitutions in equation (72) are made: 
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FIGURE 42. SPECTRAL EMISSIVITY OF ALUMINA PARTICLES 

In equation (445) the first term on the right is the attenuation (scattering plus 
absorption) of the radiant beam of intensity 1^, the second term represents 
the energy scattered in the direction 0 , and the third term is the energy 
emitted by the unit volume under consideration. 

In equation (445) the expression 7(0) is the angular distribution func- 
tion sometimes referred to as the differential scattering cross-section for 
which 


/ 


27T 7T 

J 7(0) sin0 d0 d <f> - 4 tt 
0 


and which is defined as 


( 447 ) 
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y(e) =-j [ij(fl) + i 2 ( 0 )] 


where ij(0) and i 2 (0) are the intensity functions defined in equation (387) . For 
an isotropic particle cloud, it may be noted that y(6) - 1. 

To solve equation (445) for 1^ and consequently obtain the particle cloud 
emissivity, various approximations may be considered. For a purely absorbing 
medium (a a » cr s ) , discussed as one of the solutions to the equation of transfer 
in the first section, subtitle "Solutions to the Equation of Radiation Transfer, " 
part 2, the emissivity as a function of 0 may be written as in equation (66) . 

For 0 = 90° , cos0 = p = 1, and an expression may be written for the normal 
emissivity of the cloud (assumed homogeneous) as 


£ N = 1 - e' T = 1 


> (°a.Np z ) 


(449) 


The expression for hemispherical emissivity may be found by integrating 
equation (66) over 0, and may be written according to Bartky and Bauer [34] 
as 


e = 1 - e T + re T - t 2 E^t) = e + re T - t 2 Ej(t) , (450) 

H JN 

where E t (r) is the exponential integral defined as 

oo 

Ej(t) = f (e X /x) dx . (451) 

T 

The one -dimensional beam approximation solution (for anisotropic 
scattering) to equation (72) will now be discussed for the cloud emissivity 
expressions. The details of the solution will be given in the last part of this 
section, and the expressions for emissivity, reflectivity, and transmissivity 
(e,p, and t, respectively) will be given in this section. A further approxi- 
mation can be made for the one -dimensional beam method by assuming the 
intensity scattered in the forward direction (0 = 0) is equal to the intensity 
scattered in the backward direction (0 - 180°) , or 
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(452) 
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7(0) sine d$ =~ = p =-j- 
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f 7(0) sin0 d0 
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where a represents the fraction of radiation scattered forward and /3 = 1 - a 
is the fraction of radiation scattered backward. 

In this manner, the cloud emissivity may be written according to Bartky 
and Bauer [34] and Rossler [484] as 


€( t ) = 7^(1 -e _kT )/(l + Poo e" kT ) 


(453) 


where 
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and 


k = (m) 2 



(455) 


If 2t is inserted for t in equation (453) , the one -dimensional approxi- 
mation for the hemispherical emissivity may be written. 

For the general anisotropic scattering case, the assumption that = a = 
0.5 is not made, and the expressions become much more complex. For this 
general case Bartky and Bauer [34] present the following equations for p( t) , 
t( t) , and e(r) based on the analyses of Mecke [485] and Havard [356]: 

p( T ) = (k' 2 - m 2 ) (e k T - e _k T )/D, t(r) = 4 mk’/D (456) 


and 
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e( T ) = i _ p _ t = 2m I T (k' + m) + e K T (k' - m) - 2k' I /D, (457) 


-k 




where 


o k* T 1 \2 “k'l 

D = ( k' + m) 2 e - ( k - m) 2 e 


(458) 


and 


1 

k' = j m[2(l - m) /3 + m] | 2 . 


(459) 


For an optically thick medium (t~«) , the reflectivity reduces to that 
in equation (454) , and the normal emissivity, for various limiting cases, can 
be written as 
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(460) 


The hemispherical emissivity may be written by replacing r by 2 t. 

For the three-dimensional exact (numerical) case, values of emissivity 
for isotropic scattering have been given by Bartky and Bauer [34] , using values 
of the Chandrasekhar H(/u) function, and for anisotropic scattering have been 
given by Romanova [486] . Bartky and Bauer plotted these three-dimensional 
isotropic (for a = /3 = £) and anisotropic emissivities mentioned above. Both 
hemispherical and normal emissivities were plotted as a function of the optical 
thickness for various values of the parameter m = oJo v Bartky and Bauer's 
set of curves is reproduced here as Figure 41. For an optically thin plume 
( T _*0) , Bartky and Bauer presented the limiting case of cloud emissivity for 
both one- and three-dimensional cases for a semi-infinite medium. These 
limiting values are, for isotropic scattering, 


e iD = 2m 
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= 2.31 m 
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and, for anisotropic scattering, using Romanova's [486] data, as 

i i i 

e iD ~ m2 > € 3D ~ 8 m 2 , and € 3 q = 6. 6 m 2 . (462) 

H N 

It is suggested by Adams and Colucci [424] and Carlson [5] that a 
representative value ofm for an aluminized rocket exhaust is of the order of 0. 005 
Based upon this value, Figure 41 shows that the exhaust of an aluminized rocket is 
not optically thick until t « 100 and then its emissivity is only about 0. 1 - 0. 2. 
Carlson [5] states that BeO has a value of m of about 0. 1; hence, its optically 
thick plume would have a cloud emissivity at least twice as great as that of an 
A1 2 0 3 plume. Figure 41 shows that, when the scattering cross-section ct s 
approaches zero, the value of m approaches one (which is the case for carbon 
particles) , and the optically thick cloud emissivity approaches unity (blackbody 
emission) . 

In the above analysis of one-dimensional anisotropic scattering for a 
particle plume, the effects of multiple scattering were included, but only in a 
one -dimensional sense. A two- or three-dimensional analysis of multiple 
scattering is very complex, even though it is felt by Bauer [412] and Goldstein 
[487] 'that single scattering is far greater than double scattering, except for 
very small scattering angles, 6 . Among the many investigators who have 
analyzed the effects of two- and three-dimensional multiple scattering (although 
not necessarily applied to rocket exhausts) are Chin and Churchill [488], 

Bellman et al. [489,490], Evans, Chu, and Churchill [491], Smart et al. [492], 
Uenoetal. [41], Twersky [493,494], Grosjean [495], Chu et al. [496], Sekera 
[497], Mullikin [498], Bauer [499], de Bary and Bullrich [500], Herman and 
Browning [501], Halpern, Luneberg, and Clark [502], and Richards [503], 

In the next section the effective particle temperature, on which the 
absorption and scattering cross-sections (and consequently the cloud emissivity) 
are dependent, is derived, and finally, in the last part of this section, the 
various radiation programs are discussed, as they are based upon particle cloud 
emissivity and the effective temperature. 


Determination of Particle Cloud Effective Temperature 


Since the particle cloud of an aluminized rocket exhaust is a non- 
isothermal cloud, especially near the exit plane, and has a definite thermal 
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lag between the gas and particles, it is necessary to determine some effective 
temperature at each increment in the plume through which the radiation can be 
analyzed. The particle temperature is a direct function of the radius of the 
particle, with the large particles experiencing a higher temperature than the 
smaller particles whose temperature approaches the gas temperature at the 
same point in the plume. The amount of coupling between the gas and particles 
is an important factor in determining particle temperatures. If the particle 
loadings, <f> = mp/m , are large (<p > 0. 4) , the flow may be strongly coupled 
(gas dynamic flow field influenced by particle flow field and vice-versa) , 
especially in the nozzle for small motors and at low altitudes. A case of weaker 
coupling arises when the particle flow field is influenced by the gas flow field, 
but not vice versa. A third case of completely uncoupled flow arises if both 
flow fields are independent of each other. This is the case in a highly under- 
expanded plume (of low <p) in which the mean free path becomes larger than 
the diameter of the particles and the particles lose more heat by radiation than 
by convection to the surrounding gas. 

Hoglund [504] presents an excellent review of two-phase (gas-particle) 
nozzle flow-fields, surveying references through early 1962, and McGregor 
[458] reviews the effects of two-phase flows, presenting numerous references 
through the end of 1965. Carlson, Lewis, and Bartky [505] discuss the experi- 
mental and theoretical aspects of both near-field (including location of the first 
Mach disc) and far-field gas-particle flows. The early analyses of two-phase 
flow such as by Altman and Carter [506] and Gilbert, Davis, and Altman [507] 
were concerned with effects of the thermal and velocity lags on motor specific 
impulse. In general, these investigators assumed that the thermal lag was 
independent of velocity lag, and found that thermal lag had a smaller effect on 
specific impulse than did velocity lag. Kliegel [462] performed the first one- 
dimensional nozzle two-phase flow analysis that included both thermal and 
velocity lags along with the effects of particle lag on the gas properties. Kliegel 
and Nickerson [508,509] performed the first two-dimensional (axisymmetric) 
coupled two-phase flow analysis in nozzles, and later in exhaust plumes. Since 
the early analyses of two-phase flow, numerous investigators have developed 
computer programs for both coupled and uncoupled solutions. Some of the 
many investigators in this area include Crowe et al. [510-512] , Bailey et al. 
[513], Hoffman [514-516], Rannie [517], Simons [518], Hasson [519], Glauz 
[520], Morganthaler [521], Travis [522], Rudinger [523], Marble [524], 
Torobin and Gauvin [525], Lype [526], Soo [527,528], Price et al. [529], 
Carlson [530,531] , and Carlson and Hoglund [532] . 


185 


In most of the investigations, the basic equations for particle trajectories 
and particle temperatures involve the equation of motion for the particle as 
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and the energy equation for the particle as 


dT 

E 

dt 


7r r 


P C 
P 


M C D K 
g p g H 

Pr 


(T - T 
P g 


) + 


r 4 


B 

P c 
P 


(464) 



In these expressions, pg is the gas viscosity, Vp, 
the gas velocity, Pp the particle density, and C and 


the particle velocity, 

the specific heats 
P g 


of the particle and gas, respectively. The parameters Kq and Kjj are the 
values of the ratio of the actual drag coefficient and Nusselt number, respective- 
ly, to the Stokes values. These quantities, and K H , are inserted in the 
equations to account for Stokes law departures Decause of compressibility, 
inertial, or rarefaction effects. An empirical correlation of Kp is presented 
by Carlson [530] as 
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and the quantity may be expressed as 


K H 2 1 


2 + 0 . 459 Re 
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+ 3. 42 (M/Re) (2 + 0.459 Re 0 ’ 5 *) 


(466) 


where M is the Mach number and Re is the Reynolds number. 
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Kliegel [462] obtained a simplified solution to equations (463) and (464) , 
plus the overall momentum and overall energy equations given by Hoglund [504] , 
by assuming a constant value for the ratio of gas to particle velocities and the 
ratio of gas to particle temperatures. Another simplified solution to equations 
(463) and (464) , given by Carlson [477], may be obtained by assuming that 
there is no lag between the gas and particles, and the coupling effects may be 
accounted for by using an effective isentropic exponent y and an effective 
molecular weight MW as 


1 + (p ( C/ Cp^) 

1 + yep (C/Cp ) 


and 


MW 


MW 

1 - <p 


(467) 


Fontenot [459] made a simplified analysis for the particle temperature 
(assumed constant across the plume at a distance x from the nozzle exit) as 
a function of x by combining the Stefan -Boltzmann equation with the heat flux 
relation for the particle cloud, and obtained the following relation: 
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where nip is the particle mass flow. Solving equation (468) results in the 
relation 
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where tv, is found from equation (434). The only unknown in equation (469) 
Px 

is thus T n which may be obtained from a two-phase nozzle flow program. 

Pe 

French [482] presented a curve showing the particle temperature lag 
relation (T n - T a )/(T - T s ) as a function of r^ based upon data obtained 

O C O P 

by Kliegel [462] and Carlson [532]. This curve is based upon the assumption 
that the temperature lag approaches a constant value in the diverging part of 
the nozzle and remains constant in the rest of the nozzle and in the exhaust 
plume. 
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The two-phase nozzle (and plume) program which is used to predict 
particle temperature, Tp^, and particle concentration, Np^, in Morizumi's and 

Carpenter's [463] radiation program is Klie gel's [462] method, which assumes 
five particle sizes (as discussed in the two preceding sections) . After the 
particle temperature, T D , for a particular size particle is determined from 

the two-phase plume, the average particle plume temperature along a given 
line of sight through the axis is obtained as 

i 
4 


(470) 



where the correction factor r] is found by taking the temperature variation in 
a lateral plane of view and averaging it over the lines of sight contained in that 
plane weighted by the lateral variation of the shape factor. 


Gulrajani [464] defines an effective particle temperature at any point 
in the plume as 


T 4 


P 


eff 



(471) 


and divides the plume into segments perpendicular to its axis. Equations (470) 
and (471) are used by Morizumi and Carpenter [463] and Gulrajani [464], 
respectively, in their radiation analyses, as discussed in the following section. 

The effective plume temperature analysis in the Aeronutronic radiation 
program (one -dimensional beam approximation) includes a pseudo-two- 
dimensional routine for the subsonic flow field to the throat and is combined 
with the Lockheed/Huntsville method-of-characteristics program (at the throat) . 
This program is, in turn, matched with the particle trajectory and thermal 
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history equations (463, 463, and 467) and particle size distribution equation 
(427) . In this manner, the particle momentum and energy equations (463 and 
464) are integrated using a modified Euler technique and a four-point interpo- 
lation scheme to determine the gas properties stored in an orthogonal mesh. 
The gas thermochemical data are obtained from the NASA/Lewis thermochem- 
ical programs, discussed by Svehla [533] , Zeleznik and Gordon [534] , and 
Gordon and Huff [535] . 

For an S-II ullage rocket, the particle temperatures computed by this 
program for 1-p, 2-p, and 5-p radius particles are shown in Figure 43. The 
super-cooling curve is for the case of the heat of fusion never released and 
represents an upper limit to the particle temperature. The zero super-cooling 
curve represents the heat of fusion for liquid particles being released as soon 



FIGURE 43. TEMPERATURE PROFILES FOR SATURN S-II ULLAGE 
ROCKET MOTOR AT 175,000 FEET 


as possible and represents a lower limit for the particle temperature. Figure 
44 shows the limiting particle streamlines (lines outside of which no particles 
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of a particular size exist) also computed by this Aeronutronic two-phase pro- 
gram for 1-p, 2-ju, and 5-p particles superimposed on Mach number contours 
at an altitude of 130, 000 feet. This is the altitude at which the S-II ullage 
rocket was fired in the J-4 test cell at AEDC in which radiation measurements 
were obtained (as discussed in the following section) . 



ALTITUDE SHOWING PARTICLE LIMITING STREAMLINES (AL,0 3 ) 

The effective temperature used for the radiation calculations in the 
Aeronutronic program is the temperature which is assumed in each isothermal; 
homogeneous slab viewed along the line of sight to the particulai object on 
which the radiation is incident, hi this manner, the average radiation emitted 
from a small increment of particle cloud is given by an average absorption 
cross-section, <o a >, multiplied by the Planck blackbody function, B^(T e £j-) 
at the calculated effective temperature, T eff , or 
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In equation (474) , f(r p ) is the size distribution obtained from equation (427) , 
p s is the alumina bulk density, and p g is the gas density. 

The next section discusses the incorporation of the effective tempera- 
ture and particle cloud emissivity relations in the radiation calculations as 
determined by various methods wnicn include the inverse wavelength, neutron 
scattering analogy, and one -dimensional beam approximation methods, as 
discussed in the preceding section. 


Radiation Calculations for Aluminized Solid Propellant 

Rocket Exhausts 


Since the advent of large solid propellant missiles (and smaller ullage 
and retro motors) , predictions of radiation heating from the motor exhausts 
to the surrounding base region (or to other nearby structures and components) 
have been very difficult to make. Before suitable radiation heating programs 
were available, reliance was made on crude blackbody or graybody predictions 
and on scale model and flight data. Some of the early analyses and measure- 
ments of radiant heating to the base of solid propellant missiles include, for 
the Minuteman first stage [536-549], Minuteman second stage [550-554], and 
Minuteman third stage [555-560], for the Polaris first and second stages 
[561-566], and for the 120 -inch strap- ons to the Titan 1I1C [567-573]. 


An early analysis of radiant heating for AI 2 O 3 particles using the Mie 
theory with various mass fractions of AI 2 O 3 for particle sizes ranging from 
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r p = 0. 5 to 1 0 /n was performed by Svaton and Winer [357] . The radiant 
heating equation used was a modification of equation (208) as 

q/A = n Iff °t N p B i/ T ^ exp *"/ °t N dS ^ cos6ds dc0 d ^ (475) 

v co s P v P 


where the quantity N p represented the number of particles per unit volume and 
was given as 
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where is the percentage by weight of the A1 2 0 3 particles in the plume. 

The pressure in Svaton and Winer's [357] analysis was varied from 
1-20 atm for a theoretical isobaric plume of four inches in diameter. Radiant 
intensities ( watts/ cm 2 p) were plotted as a function of wavelength for a tempera- 
ture of 2000° K. The results are not believed to be representative of a practical 
case because the value of index of refraction of the A1 2 0 3 particles used was 
n = 1. 315 - 0. 0143 i, taken from Harris [574] , a value not in agreement with the 
results of Carlson [344] and Gryvnak and Burch [422,423]. Furthermore, 
Svaton and Winder's [357] analysis was made for an arbitrary temperature of 
2000° K, with no suggestions as to the use of an effective temperature or parti- 
cle size distribution throughout the plume. 


Ohrenberger [575] presents another early analysis of particle radiation 
from an aluminized rocket exhaust plume. His final equation for the radiation 
from an isothermal cloud for particle size n to an element of area dA located 
at the vertex of the n dl cone of expansion is 



R. 
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) exp 



» 


(477) 

where p and e are the reflectivity and emissivity, respectively, of the particle 
cloud, ot n is the total cross-section of the n th size particle, and R c is the 

distance from dA to the edge of the cloud. The quantity 
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where dco is the solid angle subtended by dA from a point on the particle. The 
quantity rf is the normal to the surface dA which is located a distance R from 
the particle of size n. Using this method, Ohrenberger was able to show only 
a limited agreement with experiment as he had to make certain geometrical 
approximations, along with the assumption of an optically thin particle cloud, 
in order to make simplified integrations over the volumes considered. 

The inverse wavelength method of predicting radiation discussed by 
Fontenot [459] is based upon the Stefan-Boltzmann law for a particle: 


q /A = F a 
x 


B 
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In predicting radiation heating by this method, the particle total emissivity 
€p x is calculated by equation (434) , and the value of T p ^ is found from equation 

(469) . Calculations have been made using this method by Fontenot [576] and 
Hunt [577] for Minuteman strap-on motors for the S-IC stage and by Mullin and 
Bender [460], 461] for Saturn ullage and retro motors. In general, it appears 
that Hi is method gives higher values of radiative heating than does Morizumi's 
and Carpenter's [463] method or the Aeronutronic method. As stated previously, 
Fontenot's [459] method makes the assumptions that the particles are all of 
the same size, that the particle temperature is constant across radial cuts in 
the plume, and that the effects of scattering may be neglected. 

Morizumi's and Carpenter's [463] neutron- scattering analogy was used 
to predict radiation heating to the base of the Wing VI second stage scale model 
and full scale Minuteman vehicle and also predict radiation heating to the Vela 
Nuclear Detection Satellite [578] from its injection rocket exhaust. The basic 
equation for radiation heating which these investigators used in their numerical 
technique was 


q/A 
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a. i i i 
i 


(480) 


193 


where e a . and are, respectively, the mean value of apparent emissivity and 

effective temperature based on equations (436) and (470) for a particular (the 
ith) isothermal, homogeneous slab through which the radiation is seen. The 
shape factors AFj from the surface to be analyzed to the slabs in the plume 
are determined by methods of Morizumi [579,580] . 

Other investigators have used modifications of Morizumi's and Carpenter's 
[463] theory to predict radiation from solid propellant plumes. Gulrajani [464] 
combined Kliegel's [462] two-phase plume and simplified the radiation analysis 
by segmenting the plume normal to its axis and calculated heating rates from 
the sub scale second stage Minuteman rocket fired at AEDC [483]. He also 
presented a model which was based on an effective length of the plume equal to 
3 d e for purposes of radiation calculation. Hunt [466] programmed Morizumi's 
and Carpenter's [463] method and used the two-phase plume program of Simons 
[518] in an attempt to predict heating rates for the first stage Minuteman motor 
firings discussed by Hunt [581]. Other more simplified analyses of Morizumi's 
and Carpenter’s [463] program were used by Brower [582] to predict heating 
to the base of the S-IB stage from the exhausts of Minuteman strap-ons, and 
by Rochelle [465] to predict radiant heating to the S-IC ordnance disconnect 
caused by the exhaust of the S-II ullage rockets. French [482] also used a 
modification of Morizumi's and Carpenter's program for a generalized particle 
radiation problem based on certain velocity and thermal lag assumptions as 
discussed in the last section. 

It is felt that Morizumi's and Carpenter's method is superior to the other 
methods above; however, it is believed that the particle emissivity, e p, used 
is not equivalent to the quantity m = o- a /crt> and should not be set equal to a 
constant value of 0.25, which is felt to be more than an order of magnitude too 
high. The effects of change of index of refraction with phase state (as discussed 
in the first part of this section) , and consequently change of emissivity with 
phase state, are not considered in their analysis. Also, the assumption of five 
particle sizes with each size represented by 20 percent (by weight) of the total 
A1 2 0 3 weight fraction is not believed to be as valid as an appropriate size distri- 
bution function such as equation (427). Finally, although the effects of one- 
dimensional multiple scattering are included, the effects of anisotropic scatter- 
ing, in which the phase function y(0) is not equal to unity are not considered, 
nor is Mie theory used to calculate the appropriate extinction efficiency factor 
(a value of = 2 was assumed) . 

The particle radiation method developed at Aeronutronic under contract 
to MSFC's Thermal Environment Branch which is based upon the one -dimensional 
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beam approximation technique described in the previous section and the effective 
temperature method mentioned in the preceding part of this section, has been 
described by Kuby [583,584]. This method, which is based upon one-dimensional 
multiple anisotropic scattering with appropriate particle size distribution, is 
believed to be more accurate than any of the methods discussed previously. 

The basis of the one-dimensional beam approximation method as origi- 
nally developed by Schuster [585], modified by Hamaker [586] , Churchill et al. 
[587], and Rossler [484] and further discussed by Klein [588]* is the writing 
of the plane -parallel equation of transfer, equation (445) , in terms of a forward 
component lj (for 0° or forward scattering) and a backward component Ji (for 
180° or backward scattering) . Equation (445) thus reduces to the following 
two equations in one dimension: 

dl. 

l 


and 

dJ. 

l 

dz 

where i represents the i^h slab in the direction of the line of sight which has 
an effective temperature, particle number density, and size distribution, and 
/3 represents the back-scattering fraction obtained from equation (452) . 

A general solution to this set of equations may be found by using the 
optical thickness and ratio of absorbing to total cross-sections, m = <7 a / 
as 
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-"Additional information was obtained from an Aeronutronic unpublished paper of 
February 1966, by C. D. Bartky, entitled "The Reflectance of Homogeneous, 
Plane -Parallel Clouds of Dust and Smoke. " 
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where the quantity k' is defined in equation (459). 

If the radiation passes through n slabs to a surface or detector where 
it is measured, there would thus be 2n equations (I lf J 1# . . . . I n , J n ) with 2n 
constants (A 1} Cj,. . . .A n , C n ) which have to be determined. At the boundaries 
of the slabs, a set of continuity equations would be written as 


I. (t.) - I. , (t.) and J.(r.) = J. , (r.) 
1 1 i+I i i i i+l' i 


(485) 


At the two outer boundaries the intensity in the forward direction 1 (t - 0) 
would equal I 0 , which would be zero if viewed across a diameter of the plume 
(incident intensity at the edge of the plume away from the observer would be 
zero) , and if viewed along the axis of the plume (up the nozzle of the motor) , 
the value of I Q would be 
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The 180° scattered intensity at the outer boundary, J(r - T n ) , if viewed along 
the axis would be input as 
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Equations (483) and (484) would thus be solved to obtain the constants Aj, C l5 
etc. , and the intensity at T n , the position of the surface or detector (if inside 
the plume) on which the radiation is incident, would thus be determined. The 
corresponding heating rate would then be the intensity integrated over wavelength 
and solid angle. 

In considering a simplified solution for the intensity for a three-slab 
case viewed perpendicular to the plume axis near the nozzle exit, Kuby [583] 
writes the following equation for the radiant intensity to a detector mounted 
outside the plume. This equation, which considers the reflections and trans- 
missions for all of the slabs as discussed by Bartky and Bauer [34] and by 
Harrison [310] is 
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I = [i^ + Ii(l - Pi) h + M 1 - Pi) Ml - Pz) t-i] 


(488) 

+ j^l 2 Pi t2 + Ii Pz tj( 1 - P 2 ) taj + £l3( 1 - Pi) ti P 2 P3( i - P 2 ) h +. . .J + £ . . .J . 


The upper (p = 0) and lower (P = 1) bounds for the solution of equations 
(481) and (482) are described by Kuby [583], as well as the "best estimate" 
obtained by using the correct expression for p as a function of 0. Calculations 
of particle radiation heating using this method are currently being performed by 
Aeronutronic under the MSFC contract mentioned previously. These calculations 
will be made to compare with data obtained from both wide- and narrow-view 
angle radiation detectors as shown in Figure 45, which shows the layout of 
instrumentation associated with the S-II ullage rocket tests at AEDC , discussed 
by Rochelle [6, 589] . Comparisons of this theoretical method will be made 
with the data obtained at AEDC from the narrow (2 in. x 2 in. at 20 ft) view 
angle FF-1 radiometer presented in Figure 46, together with the wide veiw angle 
(120°) radiation calorimeter data obtained on the probe shown in Figure 47. 
Figure 48 shows the theoretical spectrum (H 2 0, CO 2 and CO emission) for the 
S-II ullage rocket calculated by the method described in the previous section 
entitled "Gaseous Radiation," subsection "Methods of Predicting Gaseous 
Radiation from Rocket Exhausts," part 3c, as seen by the FF-1 radiometer 
in Figure 45. 

Comparisons will also be made with data from the same FF-1 radiometer 
and 150° view angle radiation calorimeters for the Centaur retro tests described 
by Rochelle [590] . It is expected that comparisons will also be made with data 
obtained from the S-IB retro scale model tests at Cornell discussed by Rochelle 
[591] and Hendershot and Dennis [592-594]. Calculations of radiation heating 
to the stagnation point of copper and Teflon hemispherical-cylindrical probes 
shown in Figure 49, mounted 20 inches from the exit plane during sea level 
static firings of the S-II ullage rocket, are also expected to be made using this 
method. This particular test, described by Rochelle [7] and Datis and Fowler 
[595] resulted in total heating rates at the stagnation point of the order of 3000 
to 3200 Btu/ft? sec. Finally, calculations of heating caused by 120-in. strap-on 
motors are expected to be made to the base of the S-IC and S-IB stages. The 
120-inch strap-on configuration for the Saturn V vehicle is shown in Figure 
50. It is expected that a base heating experimental program will be performed 
for this configuration (first stage only) at Cornell Aeronautical Laboratory in 
the spring of 1967, and a combined base heating/ Launch Umbilical Tower 
impingement heating program will be started at MSFC T s Test Laboratory early 
in 1967. 
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FIGURE 45. SKETCH OF INSIDE OF J-4 TEST CELL 
SHOWING INSTRUMENTATION LOCATIONS 



FIGURE 46. RADIATION HEATING RATES FOR S-II ULLAGE ROCKET 
USING MODEL FF-1 RADIOMETER 6 INCHES FROM EXIT PLANE 
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Total Integrated Heat Transfer Rate 
from X = I.Oa to X = 200a is 3.592 BTU/ft 2 sec 
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FIGURE 47. S-II ULLAGE ROCKET EXHAUST SPECTRAL HEAT 
TRANSFER CALCULATIONS AT 120,000 FEET, AT A POSITION 
6 INCHES FROM NOZZLE EXIT 
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FIGURE 48. HTL COPPER PROBE USED IN AEDC S-II 
ULLAGE ROCKET TESTS 
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FIGURE 49. MOUNTING SUPPORT FOR COPPER PROBE USED 
IN S-II ULLAGE SEA LEVEL TESTS 

Besides the radiation analyses of the three major groups ( Boeing, TRW, 
and Aeronutronic) discussed here, it has been recently discovered by the author 
that another group at Hughes Aircraft has been actively engaged in determining 
radiation from rocket exhausts which scatter as well as emit and absorb. 

Edwards and Bobco [596] discuss solutions of the radiation transport equation 
(51) for an isothermal dispersion which emits, absorbs, and scatters iso- 
tropically. Equation (51) is transformed into a second-order differential 
equation and is solved by using the moment technique to obtain a first-order 
solution and an iterative procedure to obtain a second-order solution. Bobco 
[597] describes a closed form solution to obtain directional emissivities from 
a two-dimensional absorbing and scattering medium. Bobco'' presents an 
analysis of some of Aeronutronic 's alumina measurements and describes cal- 
culations of the parameter m = a /<t for exhausts of the Surveyor main 

at 

*R. P. Bobco, "Extinction, Absorption, and Scattering Parameters for Alumina 
Particles in Exhaust Plumes," Hughes Aircraft Interdepartmental Correspondence, 
November 22, 1966; also "Some Numerical Results on the Radiant Emission From 
Particle Exhaust Plumes," Final Report on IR and D Study of Exhaust Plume 
Radiation, Hughes Aircraft Interdepartmental Correspondence, November 30, 

1966. 
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FIGURE 50. PROPOSED SATURN V 
VEHICLE WITH FOUR 120-INCH 
STRAP-ONS 


retro motors engines. Bobco and Edwards 
[598] describe calculations of directional 
emissivity from non-uniform axisymmet- 
ric dispersions of isotropically scatter- 
ing particles in which the particle num- 
ber density, Np, is assumed to vary 
inversely with the square of the distance 
from the cone vertex. 

This concludes the section on 
radiation from A1 2 0 3 particles. It is felt 
that most of the basic points of radiation 
from aluminized rocket exhausts were 
covered. The details of the Mie theory 
for scattering phenomena, combustion 
of aluminum to A1 2 0 3 , sizes and distri- 
butions of A1 2 0 3 particles, Al 2 O s parti- 
cle cloud emissivities and effective 
temperatures, and particle radiation 
programs were discussed. The advan- 
tages and disadvantages of several pro- 
grams of predicting radiation heating 
from aluminized rocket exhausts were 
given and it is felt that the one -dimen- 
sional beam approximation method of 
Aeronutronic provides probably the 
best approach to the particle radiation 
problem. For additional details and dis- 
cussion of current problem areas involv- 
ing particle radiation and 2 -phase flow, 
it is suggested that the proceedings of 


AFRPL Two- Phase Conference held at 
Norton Air Force Base, California, on March 15-16, 1967, and the final Aero- 
nutronic report on the NASA contract mentioned previously [ 599] be consulted. 


CONCLUSIONS AND RECOMMENDATIONS 


The report has presented a review of the general theory of thermal 
radiation from liquid and solid propellant rocket exhausts. Discussions of 
practical programs for predicting gaseous radiation, radiation from carbon 
particles, and radiation from A1 2 0 3 particles were given. The equation of 
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transfer of radiation for a medium (such as a rocket exhaust) which absorbs, 
emits, and scatters radiation was derived in detail, and various simplified 
solutions (such as for a purely absorbing medium or a purely scattering medium) 
were presented. The gaseous radiation from H 2 O, CO 2 and CO bands in liquid 
propellant rocket exhausts in which the effects of scattering were neglected was 
then discussed in detail. Analyses of radiation from an accelerated charge, 
shape and broadening of spectral lines, band models, and various methods of 
predicting gaseous radiation from rocket exhausts were then treated. It was 
concluded that the generalized method with the modified Curtis-Godson approxi- 
mation, which includes the combined effects of collision and Doppler broadening 
was best suited for the prediction of gaseous radiation from rocket exhausts. 

Use is currently being made of the above generalized radiation method by 
Rocketdyne [266,293] to obtain the radial distribution of temperature, pressure, 
concentration, and emissivity across various plumes at the nozzle exit. This 
procedure, termed "zone radiometry," employs the Abel inversion technique 
and the measured values of transmissivity at the nozzle exit for each zone, 
together with the calculated values of absorption coefficient as a function of 
temperature and pressure, to obtain these radial distributions. 

Radiation from carbon particles, such as found in the exhausts of H-l 
and F-l engines, was next discussed. It was seen that these carbon particles 
radiate at near blackbody level, producing a continuous spectrum (except in 
the regions of atmospheric absorption) . One of the main difficulties in analyzing 
this type of radiation lies in obtaining an accurate mass fraction of carbon pres- 
ent in the exhaust. Recent measurements at Rocketdyne [199,266] have shown 
that this fraction lies between 0. 01 and 0. 05 depending upon the position in the 
plume, the type of hydrocarbon exhaust (gas generator, F-l engine, etc.) 
analyzed, and the O/F ratio. 

In the section on carbon particle radiation, the formation of carbon 
particles in luminous flames was discussed, and various theories were pre- 
sented which suggested possible formation mechanisms for these particles. 

The determination of carbon particle sizes (an important factor in analyses of 
scattering by these particles) was discussed briefly, and it was concluded that 
a mean value for the carbon particle radius, as obtained from luminous flames 
and both scale model and large (F-l) engine exhausts, could be taken as about 
200 A.- The dispersion equation for carbon as it relates to radiation from an 
oscillating dipole (Rayleigh scattering) was discussed, and equations for the 
real and imaginary parts of the index of refraction of carbon were developed. 

The Rayleigh theory of scattering by small particles (such as carbon) in which 
a « 1 was then discussed, various means of predicting carbon absorption 
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coefficients and emissivities were presented, and subsequently methods of 
predicting radiation from luminous flames and rocket exhausts were discussed. 

It was concluded that the experimental methods obtained at General Dynamics/ 
Convair [230,231] were probably the most accurate for predicting carbon 
absorption coefficients (and consequently radiation) , although this method was 
probably not valid at small wavelengths or for large particle sizes. For these 
cases the effects of scattering become important, and methods such as used by 
Stull and Plass [240] and Beheshti [260] should be used to obtain emissivities 
and radiation. Mention was then made of various means of predicting radiation 
to the base of the S-I and S-IB stages from the H-l engine exhausts and to the 
base of the S-IC from the F-l engine exhausts as used by MSFC and its associ- 
ated contractors. 

Radiation from AI 2 O 3 particles, such as those found in the exhausts of 
solid propellant motors, was then discussed. It was seen that the AI 2 O 3 parti- 
cles were larger than the carbon particles and had a lower emissivity, but a 
higher temperature. Since the A1 2 0 3 particles are relatively large, scattering 
effects were seen to be more important than for carbon particles, and the Mie 
theory of scattering, in which o; ~ 1 , was discussed in detail in order to derive 
expressions for total and scattering cross-sections and efficiency factors. The 
determination of A1 2 0 3 particle sizes by various experimental and theoretical 
methods was then discussed, as well as the combustion of aluminum powder in 
the propellant to A1 2 0 3 particles. It was concluded that the A1 2 0 3 particles were 
formed by a vapor-phase combustion process, and that the particles were of 
the general size (radius) of 1 to 3 p, although samples obtained from large 
motors showed sizes larger than this, and particles analyzed by optical tech- 
niques showed sizes smaller than this. The particle size distribution equation 
(427) was concluded to be the best theoretical representation of A1 2 0 3 particle 
sizes for use in radiation calculations. Various methods of predicting particle 
cloud emissivities were presented, and it was believed that the one-dimensional 
beam approximation method developed at Aeronutronic provided the most realis- 
tic analysis. This is because this method considers anisotropic, multiple 
scattering and uses the most accurate data for AI 2 O 3 particle emissivity (and 
consequently index of refraction) . A brief discussion of two-phase flow and 
particle cloud effective temperature was given; it was felt equation (472) should 
be used for the calculations since it did not depend upon the assumption of five 
sizes of particles, as did equation (470) . Discussions of radiant heating cal- 
culated by different methods were then given and various sets of experimental 
data were listed which are currently being compared with the one-dimensional 
beam method of predicting AI 2 O 3 particle radiation. 
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The ultimate objective of all of these radiation programs is a program 
which will combine the best methods mentioned above into one single program 
which can be used with any rocket exhaust. Since a solid propellant exhaust 
would be the most difficult to analyze theoretically, this ideal program would 
have the A1 2 0 3 (and possibly Fe20 3 ) absorption coefficients tabulated as a 
function of wavelength and temperature just as the current MSFC radiation 
programs are being adjusted to have the carbon, CO 2 , CO, H 2 0, and possibly 
HC1 absorption coefficients tabulated. In this combined program the graybody 
continuum emission of A1 2 0 3 (or Fe 2 0 3 ) particles and the near blackbody con- 
tinuum emission of carbon particles would be superimposed on the bands of the 
gaseous species whose radiation would be computed by the generalized method 
with modified Curtis-Godson approximation (with combined collision and Doppler 
broadening) . This is the method currently being used in an attempt to predict 
radiation heating to the copper and Teflon probes which were immersed in an 
S-II ullage rocket exhaust 20 inches from the nozzle exit, as discussed by 
Rochelle [6] and Datis and Fowler [595]. In this particular case, radiation 
from behind the bow shock in front of the probes will possibly contribute a 
sizable amount of the total heating and because of the higher temperatures 
(~ 5500° R) in this region, absorption coefficients from ionized species might 
have to be investigated. 

Most of the radiation processes discussed in this report have generally 
been of the equilibrium type, and the temperatures of all of the gaseous species 
have been assumed to be equal in the same point in the plume. In the after- 
burning mantle of a plume (especially the plume of a hydrocarbon-fueled engine) 
where chemiluminescence occurs, it is known that radiation from molecules 
such as OH, CH, and C 2 may be of the nonequilibrium type (i. e. , radiation is 
emitted when these molecules are in an excited state and the corresponding 
temperature is higher than the Planck blackbody temperature) . Recent studies 
at Rocketdyne [266] have determined that in the afterburning mantle an effective 
rotational temperature for OH at 3064 A was about 8100° R (over 3000° R higher 
than C0 2 gas temperature measured at the same point) and the effective rotational 
temperature for CH at 4312 A was about 5400° R (over 1000° R higher than the 
C0 2 gas temperature) . It is thus felt that these nonequilibrium effects, including 
determination of effective vibrational and electronic temperatures, should be 
investigated further in order to extend the limits of the current radiation pro- 
grams. Since the next generation in liquid propulsion systems will be hydrogen- 
fluorine engines, consideration should also be given to nonequilibrium (and 
equilibrium) radiation from the exhausts of these types of engines. Discussions 
of the rotational fine structure for the HF molecule in the 2. 5 to 2. 8 p vibrational 
band have been given by Simmons and Golden [600] , Golden [ 601] , Kuipers, 

Smith, and Nielson [ 602] , Deeds et al. [ 603] , Benedict et al. [ 604] , Mann [ 605] , 
and Simmons [606, 607] . 
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Further investigations involving solid propellant plumes might include 
a practical three-dimensional multiple anisotropic scattering analysis for 
rocket exhaust radiation. This type of analysis might be applicable for radiant 
heating to the base of the S-IB and S-IC stages caused by the optically thick 
exhausts of the 120-inch and Minuteman strap-ons. The problems of search- 
light emission from the inside of the combustion chamber possibly should be 
investigated further, both for solid propellant motors and for hydrocarbon-fueled 
engines. Since the analysis of radiation from A1 2 C>3 particles is only as accu- 
rate as the particle size distribution assumed, it is felt that additional work 
in particle size analysis is needed, especially with respect to the larger solid 
propellant motors and for higher chamber pressures (up to 2000 psia) . Experimen- 
tal measurements to determine the complex part, n 2 , of the index of refraction 
for liquid alumina should be obtained at wavelengths greater than 2.3 (i, the 
limit of current experimental data. Also, values of n 2 for Fe 2 C >3 should be 
obtained for temperatures up to and beyond its melting point for a wide range 
of wavelengths (visible out to 8 to 10 m) . According to Kuby [608] current 
optical measurements show that, at \ = 1 ju and at room temperature, the value 
of n 2 for Fe 2 C >3 is about 15 times that of A1 2 0 3 at the A1 2 0 3 melting point and 
over seven orders of magnitude higher than the room temperature value of n 2 
for A1 2 0 3 . Hence, if the Fe^ value of n 2 is found to increase with tempera- 
ture as does the A1 2 0 3 value, the Fe 2 0 3 will be found to be a stronger radiation 
emitter than A1 2 0 3 . 

The analysis of scaling effects for radiation measurements obtained 
from plumes of small motors fired in altitude cells, such as those at Cornell, 
should be investigated in detail in order to make better comparison with full 
scale test and flight data. Additional measurements in plumes of velocity, 
pressure, and temperature such as those performed in the UTC 120 motor 
plume by Lai [ 609] and in the Minuteman and larger solid motors by McGee 
[610] and Johnson et al. [611] should be obtained to correlate with theoretical 
flow fields which are prerequisites for radiation calculations. Radiation mea- 
surements using heat transfer probes which have radiation calorimeters mounted 
on the stagnation line or point and which are immersed directly in the rocket 
exhaust, such as those described by Rochelle [6] and [589], should continue 
to be made for additional correlations of experiment with theory. The analysis 
of radiation behind normal and oblique shocks, especially for solid propellant 
plumes, should be investigated further, since this could be a major source of 
radiative heating to an object inside a plume (especially at low altitudes) . 

Also, thorough investigations should be made of radiation effects of intersecting 
plumes of solid propellant motors (including effects of clustered nozzles, as well 
as plumes which intersect some distance downstream) . Finally, flow fields 
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from 120-inch and Minuteman strap-on solid propellant motors which intersect 
H-l and F-l liquid propellant plumes should be analyzed thoroughly so that 
radiation from these composite flow fields can be predicted for the S-IB and 
S-IC base regions. 

In conclusion, it is felt that this memorandum has provided a detailed 
review of the various theories of predicting radiation from liquid and solid 
propellant rocket exhausts, describing what is felt to be the most accurate 
methods for predicting gaseous radiation, radiation from carbon and A1 2 0 3 
particles, and suggesting possible improvements in these general theories 
and methods. 


George C. Marshall Space Flight Center 

National Aeronautical and Space Administration 
Huntsville, Alabama, April 26, 1967 
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